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1 Grand View 
1.1 Preface 

Flavour dynamics and the related origin of quarks and lepton masses and mixings are among 
the least understood topics in the elementary particle physics. While the definite understanding 
of flavour dynamics will most probably come from a fundamental theory at very short distance 
scales, as the GUT scale or the Planck scale, it is commonly accepted that the study of CP- 
violating and rare decay processes plays an important role in the search for this fundamental 
theory. 

In this context an important issue is the question whether the Standard Model (SM) of 
fundamental interactions is capable of describing the violation of CP symmetry observed in 
nature. Actually this question has already been answered through the studies of a dynamical 
generation of the baryon asymmetry in the universe, which is necessary for our existence. It 
turns out that the size of CP violation in the SM is too small to generate a large enough 
matter-antimatter asymmetry observed in the universe today. 

On the other hand it is conceivable that the physics responsible for the baryon asymmetry 
involves only very short distance scales, as the GUT scale or the Planck scale, and the related CP 
violation is unobservable in the experiments performed by humans. Yet even if such an unfor- 
tunate situation is a real possibility, it is unlikely that the SM provides an adequate description 
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of CP violation at scales accessible to experiments peformed on our planet in this millennium. 
On the one hand the Kobayashi-Maskawa (KM) picture of CP violation is so economical that 
it is hard to believe that it will pass future experimental tests. On the other hand almost any 
extention of the SM contains additional sources of CP violating effects. As some kind of new 
physics is required in order to understand the patterns of quark and lepton masses and mixings 
and generally to understand the flavour dynamics, it is very likely that this physics will bring 
new sources of CP violation modifying KM picture considerably. 

Similarly to CP violation, particle-antiparticle mixing and rare decays of hadrons and leptons 
play an important role in the tests of the SM and of its extentions. As particle-antiparticle mixing 
and rare decay branching ratios depend sensitively on the masses and couplings of particles 
involved, these transitions constitute an excellent machinery to study the flavour dynamics of 
quarks, leptons and other particles like sparticles in the supersymmetric extentions of the SM. 

As of January 2001 all existing data on CP violation and rare decays can be described by 
the SM within the theoretical and experimental uncertainties. An important exception are the 
neutrino oscillations, which implying non-vanishing neutrino masses changed the SM picture 
in the lepton sector considerably. It is exciting that in the coming years the new data on CP 
violation and rare decays as well as B® — B® mixing coming from a number of laboratories in 
Europe, USA and Japan may change the SM picture in the quark sector as well. 

These lectures provide a rather non-technical description of this fascinating field. There 
is unavoidably an overlap with our Les Houches Q and Lake Louise lectures [|| and with the 
reviews [||] and 0. On the other hand new developments since the summer 1999 have been taken 
into account, as far as the space allowed for it, and all numerical results have been updated. 
Moreover the discussions of various types of CP violation and of the physics beyond the SM have 
been considerably extended. In particular we discuss in detail the models with minimal flavour 
violation, presenting an improved lower bound on the angle (3 in the unitarity triangle. Finally 
we provide the complete list of references to NLO calculations for weak decays performed until 
the end of 2000. 

The first decade of the new millennium began strictly speaking one month ago. It is a 
common expectation that this decade will bring important, possibly decisive, insights into the 
structure of flavour dynamics that can be most efficiently studied through rare and CP-violating 
decays. We hope that these lecture notes will be helpful in following the new developments. In 
this respect the recent books [||, ^, [5]] , the working group reports || ||| and most recent reviews 



[10] are strongly recommended. 
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1.2 Some Facts about the Standard Model 



In the first eight sections of these lectures we will dominantly work in the context of the SM 
with three generations of quarks and leptons and the interactions described by the gauge group 
SU(3)c <8> SU(2)l (8) U(1)y spontaneously broken to SU(3)c <8> U(1)q. There are excellent text 
books on the dynamics of the SM [|ll]]-[15]. Let us therfore collect here only those ingredients 
of this model which are fundamental for the subject of these lectures. 

• The strong interactions are mediated by eight gluons G a , the electroweak interactions by 
W ± , Z° and 7. 

• Concerning Electroweak Interactions, the left-handed leptons and quarks are put into 
SU(2) L doublets: 





;i.2) 



L \ / L \ J L 

with the corresponding right-handed fields transforming as singlets under SU(2)l- The 

primes in (|1.2j) will be discussed in a moment. 



The charged current processes mediated by are flavour violating with the strength 
of violation given by the gauge coupling g<i and effectively at low energies by the Fermi 
constant 

GF 91 [1.3) 



9l 



y/2 8M? 



w 



and a unitary 3x3 CKM matrix. 



The CKM matrix [16, 17] connects the weak eigenstates (d',s',b') and the corresponding 



mass eigenstates d, s, b through 
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s' 




\ b> J 


\ 



V c d V cs V c b 
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( d\ 
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= VcKM 
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J 









:i-4) 



In the leptonic sector the analogous mixing matrix is a unit matrix due to the masslessness 
of neutrinos in the SM. 

The unitarity of the CKM matrix assures the absence of flavour changing neutral current 
transitions at the tree level. This means that the elementary vertices involving neutral 
gauge bosons (G a , Z°, 7) and the neutral Higgs are flavour conserving. This property is 
known under the name of GIM mechanism lli 
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The fact that the Py's can a priori be complex numbers allows CP violation in the SM 



An important property of the strong interactions described by Quantum Chromodynamics 



(QCD) is the asymptotic freedom [19]. This property implies that at short distance scales 
H > 0(1 GeV) the strong interaction effects in weak decays can be evaluated by means of 
perturbative methods with the expansion parameter a-jjg(fi) |20[ |. The existing analyses 
of high energy processes give a-jjg(Mz) = 0.118 ± 0.003 [21]. 



• At long distances, corresponding to fi < 0(1 GeV), ajjg(fi) becomes large and QCD effects 
in weak decays relevant to these scales can only be evaluated by means of non-perturbative 
methods. 

1.3 CKM Matrix 
1.3.1 General Remarks 

We know from the text books that the CKM matrix can be parametrized by three angles and 
a single complex phase. This phase is necessary to describe CP violation within the framework 
of the SM. 

Many parametrizations of the CKM matrix have been proposed in the literature. The clas- 



sification of different parametrizations can be found in [22|. We will use two parametrizations 



in these lectures: the standard parametrization [23] recommended by the Particle Data Group 



[24] and the Wolfenstein parametrization pa]. In the context of the models for fermion masses 
and mixings a useful parametrization has been proposed by Fritzsch and Xing |26|]. In this 
parametrization, in contrast to the standard and the Wolfenstein parametrization, the complex 
phase recides only in the 2x2 submatrix involving u, d, s and c quarks. 

1.3.2 Standard Parametrization 

With Cjj = cosOij and Sij = sin^jj = 1,2,3), the standard parametrization is given by: 



Vckm 



C12C13 S12C13 s 13 e 



-i5 \ 



-S12C23 - Ci 2 S23Si 3 e lS C12C23 - Si2S23Sl3e i<5 S23C13 
y S12S23 - Ci 2 C23Sl3e l<5 -S23C12 - Si2C 2 3Sl 3 e lS C23C13 j 



;i.5) 



where 5 is the phase necessary for CP violation. cy and can all be chosen to be positive and 
5 may vary in the range < 5 < 2ir. However, the measurements of CP violation in K decays 
force 5 to be in the range < 5 < n. 
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From phenomenological applications we know that S13 and S23 are small numbers: O(10~ 3 ) 
and 0(1O~ 2 ), respectively. Consequently to an excellent accuracy C13 = C23 = 1 and the four 
independent parameters are given as 



S12 



IK. 



Sl3 = \Vub\, S 2 3 = \V cb \, 5. 



;i.6) 



The first three can be extracted from tree level decays mediated by the transitions s — ► u, 
b —* ti and 6 — > c respectively. The phase (5 can be extracted from CP violating transitions or 
loop processes sensitive to |Ptd|. The latter fact is based on the observation that for < 5 < ir, as 
required by the analysis of CP violation in the K system, there is a one-to-one correspondence 
between 5 and \ Vtd\ given by 



\V td \ = Va 2 + b 2 -2abcos5, 



\V c dVcb\, 



\VudVub\ 



1.7) 



The main phenomenological advantages of (|T1^) over other parametrizations proposed in the 
literature are basically these two: 

• S12, S13 and S23 being related in a very simple way to \V US \, \V u b\ and \V c b\ respectively, 
can be measured independently in three decays. 

• The CP violating phase is always multiplied by the very small S13. This shows clearly the 
suppression of CP violation independently of the actual size of 5. 

For numerical evaluations the use of the standard parametrization is strongly recommended. 
However once the four parameters in ( |1.6| ) have been determined it is often useful to make a 
change of basic parameters in order to expose the structure of the results more transparently. 



This brings us to the Wolfenstein parametrization [25] and its generalization given in [27] 



1.3.3 Wolfenstein Parameterization 

The Wolfenstein parametrization is an approximate parametrization of the CKM matrix in which 
each element is expanded as a power series in the small parameter A = \ V US \ = 0.22, 



V 



1 



2 



irj) 



A 

1 2 
-AX 2 



AX 3 (g - ir/) ' 



AX 2 
1 



+ 0(X 4 



li 



y ax 3 {i - e 

and the set ( |1.6|) is replaced by 

A, A, q, V . (1.9) 

Because of the smallness of A and the fact that for each element the expansion parameter is 
actually A 2 , it is sufficient to keep only the first few terms in this expansion. 
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The Wolfenstein parameterization is certainly more transparent than the standard parametriza- 
tion. However, if one requires sufficient level of accuracy, the higher order terms in A have to 
be included in phenomenological applications. This can be done in many ways. The point is 



that since ( |1.8| ) is only an approximation the exact definiton of the parameters in (1.9) is not 
unique by terms of the neglected order 0(A 4 ). This situation is familiar from any perturbative 
expansion, where different definitions of expansion parameters (coupling constants) are possible. 
This is also the reason why in different papers in the literature different 0(A 4 ) terms in (|1.8|) 
can be found. They simply correspond to different definitions of the parameters in Ql.9| ). Since 
the physics does not depend on a particular definition, it is useful to make a choice for which 
the transparency of the original Wolfenstein parametrization is not lost. Here we present one 
way of achieving this. 

1.3.4 Wolfenstein Parametrization Beyond LO 

An efficient and systematic way of finding higher order terms in A is to go back to the standard 



parametrization ( |ljj ) and to define the parameters (X,A,g,rj) through [27, 28] 

si2 = A, s 23 = ^A 2 , s 13 e- iS = A\ z { e -irj) (1.10) 
to all orders in A. It follows that 

coso, 7] = smo. (1-11) 



S12S23 S12S23 
( 1.10| ) and ( |l.ll| ) represent simply the change of variables from fll.6| ) to ( |l.9[) . Making this 



change of variables in the standard parametrization ( |1.5| ) we find the CKM matrix as a function 
of (A, A, g, rj) which satisfies unitarity exactly. Expanding next each element in powers of A we 
recover the matrix in (|1.8| ) and in addition find explicit corrections of 0(A 4 ) and higher order 
terms: 

V ud = 1 - \x 2 - \x A + 0(A 6 ) (1.12) 
V us = X + 0(X 7 ), V ub = AX 3 (e-ir]) (1.13) 
Vcd = -A + ^A 2 X 5 [1 - 2(g + i v )] + 0(X 7 ) (1.14) 

V cs = l- \x 2 - \x\l + AA 2 ) + 0(A 6 ) (1.15) 

2, o 

V cb = AX 2 + 0(X 8 ), V tb = l-^A 2 X 4 + 0(X G ) (1.16) 



V U = AX 



l-{Q + iri)(l-\>?) 



+ 0(X 7 ) (1.17) 
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V ts = -AX 2 + ±A(1 - 2g)X* - i v AX 4 + 0(A 6 ) . (1.18) 

We note that by definition V ub remains unchanged and the corrections to V us and V cb appear 
only at 0(X 7 ) and 0(A 8 ), respectively. Consequently to an an excellent accuracy we have: 

V us = A, V cb = AX 2 , (1.19) 

V ub = AX 3 (g-i V ), v td = AX 3 (l-Q-ifj) (1.20) 

with |27| 

e=e(l-y), V = r,(l-^-). (1.21) 

The advantage of this generalization of the Wolfenstein parametrization over other generaliza- 
tions found in the literature is the absence of relevant corrections to V us , V cb and V ub and an 
elegant change in Vtd which allows a simple generalization of the so-called unitarity triangle 
beyond LO. For these reasons this generalization of the Wolfenstein parametrization has been 
adopted by most authors in the literature. 

Finally let us collect useful approximate analytic expressions for Aj = Vi d V* s with i = c,t: 

lmXt = -lmX c = r]A 2 X 5 =\ V ub \ \V cb \ sin 5 , (1.22) 



A 2 

- 

A_ 2 

2 

Expressions (|i~2l) and flOf ) represent to an accuracy of 0.2% the exact formulae obtained 



ReA c = -A(l - — ) , (1.23) 
RcA, - ! J -y)i 2 A 5 (l-p). (1.24) 



using ( |1.5| ). The expression ( 1.24 ) deviates by at most 2% from the exact formula in the full 
range of parameters considered. For g close to zero this deviation is below 1%. After inserting 
the expressions ( |l.22| )-( [l.24| ) in the exact formulae for quantities of interest, a further expansion 
in A should not be made. 

1.3.5 Unitarity Triangle 

The unitarity of the CKM-matrix implies various relations between its elements. In particular, 
we have 

V ud V: b + V cd V; b + V td V t * b = 0. (1.25) 

Phenomenologically this relation is very interesting as it involves simultaneously the elements 
V ub , V cb and V td which are under extensive discussion at present. 

The relation ( pT25| ) can be represented as a "unitarity" triangle in the complex (g, rf) plane. 



The invariance of ( 1.25 ) under any phase-transformations implies that the corresponding triangle 
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is rotated in the (g, fj) plane under such transformations. Since the angles and the sides (given 
by the moduli of the elements of the mixing matrix) in this triangle remain unchanged, they 
are phase convention independent and are physical observables. Consequently they can be 
measured directly in suitable experiments. One can construct additional five unitarity triangles 



corresponding to other orthogonality relations, like the one in ( |1.25 ). They are discussed in [p9|| . 
Some of them should be useful when LHC-B and BTeV experiments will provide data. The 



areas of all unitarity triangles are equal and related to the measure of CP violation Jcp [3C, pip : 

I Jcp 1= 2 • A A , (1.26) 

where denotes the area of the unitarity triangle. 

The construction of the unitarity triangle proceeds as follows: 



• We note first that 

V cd V: b = -AX 3 + 0(X 7 ). 
Thus to an excellent accuracy V cd V* b is real with IV^V^I = AX 3 . 



(1.27) 



Keeping 0(A 5 ) corrections and rescaling all terms in ( |1.25 ) by AX 3 we find 



1 



v ud v ub = Q + *ri, 



AX 3 ud 

with g and fj defined in ( p.. 21 ). 



AX 3 



V td V t t = l-(g + ifj) 



(1.28) 



Thus we can represent ( 1.25| ) as the unitarity triangle in the complex (g, rf) plane as shown 
fig- i 



m 



A=(p,r|) 



p+iri 




C=(0,0) 



1-p-iT] 



B=(1,0) 



Figure 1: Unitarity Triangle. 



Let us collect useful formulae related to this triangle: 
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Using simple trigonometry one can express sin(20j), = a, /?, 7, in terms of (^>, 77) as 
follows: 

27? (if + g 2 - g) 



sin(2a) 



sin(2/3) 
sin(27) = 



{g 2 +f) 2 ){{l-g) 2 +f) 2 )' 

2r?(l - g) 



' (l_^)2 + ^2' 

2^f/ 2gn 



g 2 +T] 2 g 2 + r/ 2 



(1.29) 
(1.30) 
(1.31) 



The lengths CA and i?^4 in the rescaled triangle to be denoted by R b and R t , respectively, 
are given by 

IT/ V*\ 1 X 2 t y 

(1.32) 



At = 



K d v; b \ ' 

1^1 



g2 + f) 2 = {l __ ) _ 



1 — £>) 2 + T] 2 



v td 



V, 



cb 



(1.33) 



The angles /? and 7 of the unitarity triangle are related directly to the complex phases of 
the CKM-elements Vtd and V ub , respectively, through 



V, 



td 



\v t d\<r i0 , v ub = \v ub \e-^. 



(1.34) 



• The angle a can be obtained through the relation 

a + /J + 7 = 180° (1.35) 
expressing the unitarity of the CKM-matrix. 

The triangle depicted in fig. [I], \V US \ and \V cb \ give the full description of the CKM matrix. 
Looking at the expressions for R b and Rt, we observe that within the SM the measurements of 
four CP conserving decays sensitive to | V us | , | V ub | , | V cb | and | V td | can tell us whether CP 
violation (77 7^ 0) is predicted in the SM. This fact is often used to determine the angles of the 
unitarity triangle without the study of CP violating quantities. 

Indeed, measuring the ratio jV^/V^I in tree-level B decays and \ Vtd\ through B^ — mixing 
allows to determine R b and Rt respectively. If so determined R b and Rt satisfy 

1 - R b < Rt < 1 + R b (1.36) 

then f\ is predicted to be non-zero on the basis of CP conserving transitions in the B-system 
alone without any reference to CP violation discovered in Kl — > it + it~ in 1964 p^ ]. Moreover 
one finds 

fj = ±^Rj^i 2 , - e= l±3_R. (1.37) 
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1.4 Grand Picture 



What do we know about the CKM matrix and the unitarity triangle on the basis of tree level 
decays? A detailed answer to this question can be found in the reports of the Particle Data 
Group |24j as well as other reviews to be mentioned in Section 4, where references to the relevant 
experiments and related theoretical work can be found. Using the information given there we 
find in particular 



IK, 



\Vub\ 

\V cb \ 



A = 0.2205 ±0.0018 

0.085 ±0.018, |y u6 | 



| Vcb | = 0.041 ±0.002, 
= (3.49 ±0.76) • 10~ 3 . 



Using ( |j~T9| ) and (|0|) we find then (A = 0.22) 

A = 0.847 ± 0.041, R b = 0.38 ± 0.08 . 



(1.38) 
(1.39) 

(1.40) 



This tells us only that the apex A of the unitarity triangle lies in the band shown in fig. ||. In 
order to answer the question where the apex A lies on this "unitarity clock" we have to look 
at different decays. Most promising in this respect are the so-called "loop induced" decays and 
transitions and CP-violating B decays which will be discussed in these lectures. 



0.5 



T| 



-0.5 



-0.5 0.5 



P 



Figure 2: "Unitarity Clock". 

These two different routes for explorations of the CKM matrix and of the related unitarity 
triangle may answer the important question, whether the Kobayashi-Maskawa picture of CP 
violation is correct and more generally whether the Standard Model offers a correct description 
of weak decays of hadrons. Indeed, in order to answer these important questions it is essential to 
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calculate as many branching ratios as possible, measure them experimentally and check if they all 
can be described by the same set of the parameters (A, A, g, rj). In the language of the unitarity 
triangle the question is whether the various curves in the (g, fj) plane extracted from different 
decays and transitions will cross each other at a single point as shown in fig. 3 and whether the 
angles 7) in the resulting triangle will agree with those extracted from CP-asymmetries in 
B decays and CP-conserving B decays. It is truely exciting that during the present decade we 
should be able to answer all these questions and in the case of the inconsistencies in the (g, fj) 
plane get some hints about the physics beyond the SM. One obvious inconsistency would be the 
violation of the constraint ( 1.36| ). 



1.5 

T| 1 
0.5 



K^7I°VV 
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ub cb 



■1 




B° 


-B° 




B 




X d vv 
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X d ll 


B d 
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K + -* 71 vv 



Figure 3: The ideal Unitarity Triangle. 



Clearly the plot in fig. || is highly idealized because in order to extract such nice curves 
from various decays one needs perfect experiments and perfect theory. One of the goals of these 
lectures is to identify those decays for which at least the theory is under control. For such decays, 
if they can be measured with a sufficient precision, the curves in fig. [| are not fully unrealistic. 
Let us then briefly discuss the theoretical framework for weak decays. 



2 Theoretical Framework 

2.1 OPE and Renormalization Group 

The basis for any serious phenomenology of weak decays of hadrons is the Operator Product 
Expansion (OPE) p3, 34 1, which allows to write the effective weak Hamiltonian simply as 
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follows 

Gp 



(2.1) 



Here is the Fermi constant and Qi are the relevant local operators which govern the decays 
in question. They are built out of quark and lepton fields. The Cabibbo-Kobayashi-Maskawa 



factors V^ KM [16, 17] and the Wilson coefficients Cj(/t) J33| describe the strength with which 
a given operator enters the Hamiltonian. The latter coefficients can be considered as scale 
dependent "couplings" related to "vertices" Qi and as discussed below can be calculated using 
perturbative methods as long as // is not too small. 

An amplitude for a decay of a given meson M = K,B, .. into a final state F = ttvu, irn, DK 
is then simply given by 

A (M - F) = (F\H eff \M) = ^ E V&KMCifamQiMlM), (2.2) 

where (F\Qi(fj,)\M) are the matrix elements of Qi between M and F, evaluated at the renormal- 
ization scale //. 

The essential virtue of OPE is this one. It allows to separate the problem of calculating 
the amplitude A(M — » F) into two distinct parts: the short distance (perturbative) calculation 
of the coefficients Cj(/i) and the long-distance (generally non-perturbative) calculation of the 
matrix elements (Qj(/i)). The scale // separates, roughly speaking, the physics contributions into 
short distance contributions contained in Ci(/i) and the long distance contributions contained 
in (Qi (//)}. Thus Cj include the top quark contributions and contributions from other heavy 
particles such as W-, Z-bosons and charged Higgs particles or supersymmetric particles in the 
supersymmetric extensions of the SM. Consequently Cj(/t) depend generally on mt and also 
on the masses of new particles if extensions of the SM are considered. This dependence can 
be found by evaluating so-called box and penguin diagrams with full W-, Z-, top- and new 
particles exchanges and properly including short distance QCD effects. The latter govern the 
//-dependence of Ci(//). 

The value of // can be chosen arbitrarily but the final result must be //-independent. There- 
fore the //-dependence of Cj(//) has to cancel the //-dependence of (Qj (//)). In other words it is 
a matter of choice what exactly belongs to Cj(/t) and what to (Qi (/*))• This cancellation of the 
//-dependence involves generally several terms in the expansion in ( ^.2[) . The coefficients Ci(fj) 
depend also on the renormalization scheme. This scheme dependence must also be canceled 
by the one of (Qj(/x)) so that the physical amplitudes are renormalization scheme indepen- 
dent. Again, as in the case of the //-dependence, the cancellation of the renormalization scheme 
dependence involves generally several terms in the expansion (|2.2|). 



12 



Although n is in principle arbitrary, it is customary to choose [/, to be of the order of the 
mass of the decaying hadron. This is 0{m\ ) ) and 0{m c ) for B decays and D decays respectively. 
In the case of K decays the typical choice is fi = 0(1 — 2 GeV) instead of 0(i71k), which is 
much too low for any perturbative calculation of the couplings Cj. Now due to the fact that 
[i <C Myyz, mt, large logarithms lnM\v/^i compensate in the evaluation of Ci(fi) the smallness 
of the QCD coupling constant a s and terms a"(lnMw//i) n , a™(lnM\v/V) n_1 etc. have to be 
resummed to all orders in a s before a reliable result for Cj can be obtained. This can be done 
very efficiently by means of the renormalization group methods. The resulting renormalization 
group improved perturbative expansion for Ci(fi) in terms of the effective coupling constant 
at s (fJ>) does not involve large logarithms and is more reliable. The related technical issues are 
discussed in detail in |IJ and [|| . 

All this looks rather formal but in fact should be familiar. Indeed, in the simplest case of 
the /3-decay, T~t e ff takes the familiar form 

Hf} f = cos0 c [u 7m (1 - lb )d ® erf {I - j 5 )u e ] , (2.3) 

where V u d has been expressed in terms of the Cabibbo angle. In this particular case the Wilson 
coefficient is equal unity and the local operator, the object between the square brackets, is given 
by a product of two V — A currents. Equation (|2.3|) represents the Fermi theory for /^-decays as 



formulated by Sudarshan and Marshak [35] and Feynman and Gell-Mann |}6| more than forty 
years ago, except that in ( |2.3| ) the quark language has been used and following Cabibbo a small 
departure of V u d from unity has been incorporated. In this context the basic formula (|2.l|) can 
be regarded as a generalization of the Fermi Theory to include all known quarks and leptons as 
well as their strong and electroweak interactions as summarized by the SM. 

Due to the interplay of electroweak and strong interactions the structure of the local operators 
is much richer than in the case of the /3-decay They can be classified with respect to the Dirac 
structure, colour structure and the type of quarks and leptons relevant for a given decay. Of 
particular interest are the operators involving quarks only. In the case of the AS = 1 transitions 
the relevant set of operators is given as follows: 

Current Current : 

Qi = {s a up)v-A (upd a ) V -A Q2 = (su) V -A (ud) V - A (2.4) 
QCD Penguins : 

Q3 = (sd)v-A (QQ)v-A Qi = (s a dp)v„ A {QpQa)v-A (2.5) 

q=u,d,s q=u,d,s 

Q 5 = (s<£)v-a (.9q)v+a Qg = (s a d/3) V -A {qpqa)v+A (2.6) 

q=u,d,s q=u,d,s 
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Electroweak Penguins : 



= T, (sd) V -A Y e i (9Q)v+a Qs = X {s a d(3) V -A Y e g(Q/3Qa)v+A (2.7) 




Q 9 = o 



(sd)y_A XI e g(^)v-A QlO = T, (Sad/3)v-A Y e Q (Wa)v-A ■ (2-8) 




Here, a, (3 denote colours and e q denotes the electric quark charges reflecting the electroweak 
origin of Q7, . . . , Q w . Finally, (su)v-A = ^7^(1 - ^>)u a . 

Clearly, in order to calculate the amplitude A{M — * F) the matrix elements (Qi(fM)) have 
to be evaluated. Since they involve long distance contributions one is forced in this case to use 
non-perturbative methods such as lattice calculations, the 1/N expansion (N is the number of 
colours), QCD sum rules, hadronic sum rules, chiral perturbation theory and so on. In the case of 
certain B-meson decays, the Heavy Quark Effective Theory (HQET) also turns out to be a useful 
tool. Needless to say, all these non-perturbative methods have some limitations. Consequently 
the dominant theoretical uncertainties in the decay amplitudes reside in the matrix elements 



The fact that in many cases the matrix elements (Qi(/J,)) cannot be reliably calculated at 
present, is very unfortunate. One of the main goals of the experimental studies of weak decays 
is the determination of the CKM factors Vckm and the search for the physics beyond the SM. 
Without a reliable estimate of (Qi(fi)) this goal cannot be achieved unless these matrix elements 
can be determined experimentally or removed from the final measurable quantities by taking 
suitable ratios and combinations of decay amplitudes or branching ratios. Flavour symmetries 
like SU(2)p and SU(3)f relating various matrix elements can be useful in this respect, provided 
flavour breaking effects can be reliably calculated. However, this can be achieved rarely and 
often one has to face directly the calculation of (Qi(fi)). We will discuss these problems later 
on. 

One of the outstanding issues in the calculation of (Qi(n)} is the compatibility ("matching") 
of (Qi(ff)) with Ci(ff). (Qi(ff)) have to carry the correct \i and renormalization scheme depen- 
dence in order to ensure the \x and scheme independence of physical amplitudes. Most of the 
non-perturbative methods struggle still with this problem. Moreover, it has been emphasised 
recently in (37| that the presence of higher dimensional operators can in the case of low matching 
scales complicate further this issue. It appears to me that in the future lattice methods have 
the best chance to get the matching in question under control. On the other hand, analytic 
solutions would certainly be preferable. 



(Q<(m)>- 
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2.2 Wilson Coefficients at NLO 



In order to achieve sufficient precision for the theoretical predictions it is desirable to have 
accurate values of Ci(/x) . Indeed it has been realized at the end of the 1980's that the lead- 
ing term (LO) in the renormalization group improved perturbation theory, in which the terms 
a™(lnM\y//u)"' are summed, is generally insufficient and the inclusion of next-to-leading correc- 
tions (NLO) corresponding to summing the terms a"(lnMw//u) n_1 is necessary. In particular, 
the proper matching of Cj(/i) and (Qi(fj)) discussed above can only be done meaningfully after 
NLO corrections have been taken into account. One finds then that unphysical fi- and renor- 
malization scheme dependences in the decay amplitudes and branching ratios resulting from the 
truncation of the perturbative series are considerably reduced by including NLO corrections. It 
is then instructive to discuss briefly the general formulae for Cj(/i) at the NLO level. Detailed 
exposition can be found in Q and |ij . 

The general expression for Ci(jj) is given by: 

C( f i) = U(fi,M w )C(M w ) (2.9) 

where C is a column vector built out of Cj's. C(Mjy) are the initial conditions which depend 
on the short distance physics at high energy scales. In particular they depend on m t . We set 
the high energy scale at M\v, but other choices are clearly possible. U(n, Mjy), the evolution 
matrix, is given as follows: 

U(H,M W ) = T 9 exp 

with g denoting the QCD effective coupling constant and T g an ordering operation defined in 
|l]]. (3(g) governs the evolution of g and 7 is the anomalous dimension matrix of the opera- 
tors involved. The structure of this equation makes it clear that the renormalization group 
approach goes beyond the usual perturbation theory. Indeed U(fi,M]y) sums automatically 
large logarithms log Mw j [i which appear for /j, <C My/. In the so-called leading logarithmic 
approximation (LO) terms (g 2 log Mw/fj,) n are summed. The next-to- leading logarithmic cor- 
rection (NLO) to this result involves summation of terms (g 2 ) n (log M\y / ^) n_1 and so on. This 
hierarchic structure gives the renormalization group improved perturbation theory. 

As an example let us consider only QCD effects and the case of a single operator so that 
(|2.9|) reduces to 

C(fi) = U(fi, M W )C(M W ) (2.11) 
with C((i) denoting the coefficient of the operator in question. Keeping the first two terms in 
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the expansions of 7(5) and (3(g) in powers of g: 

2 



7(9) 



7 (0)^ + 7 (D 

47T 



4tt 



13(g) = -/% 



16vr 2 



(16vr 



2\2 



(2.12) 



and inserting these expansions into (2.10) gives: 
U(fi,M w ) 



47T 



a s (M, 



a,(M, 



w , 



Air 



J 



where 



P 



7 



(0) 



2/?o 



J 



7 



(i) 



2/?o 



(2.13) 



(2.14) 



General formulae for U(ii,M\y) in the case of operator mixing and valid also for electroweak 
effects can be found in ^ . The leading logarithmic approximation corresponds to setting J = 
in Q. 

At NLO, C(M W ) is given by 



C(M W ) = C + ^l Cl 

47T 



(2.15) 



where Cq and C\ depend generally on m t , M\y and the masses of the new particles in the 
extentions of the SM. It should be stressed that the renormalization scheme dependence of C\ 
is canceled by the one of J in the last square bracket in Q2.13] ). The scheme dependence of J in 
the first square bracket in ( [2.13 ) is canceled by the scheme dependence of (Q(fi)). The power P 
is scheme independent. The methods for the calculation of U(fj,, Mw) and the discussion of the 
cancellation of the fi- and scheme dependence are presented in detail in O. 

As an example consider the case of the operator (bd)v-A(bd)v—A relevant for B® — B® mixing. 
In this case using the so-called NDR renormalization scheme one has []39| , |52[| 

6/25 



U(m h ,M w ) 



1 + 1.63 



a s (m h ) 
Air 



a s (M w ) 
a s (m h ) 



1 - 1.63 



« S (M W ) 

47T 



0.86 . 



(2.16) 



where we have used / = 5, a s (Mz) = 0.118 and a s (m^) = 0.222. The departure of C/(mb, Mw) 
from unity is rather small in this example, due to the small value of the power P. In the case 
of (V — A) ® (V + A) operators the renormalization group effects are larger. 



2.3 Status of the NLO Calculations 
2.3.1 General Comments 

During the last decade the NLO corrections to Ci(n) have been calculated within the SM for 
the most important and interesting decays. They will be taken into account in these lectures. 
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In table [T] we give references to all NLO calculations within the SM done until the end of 
2000. While these calculations improved considerably the precision of theoretical predictions in 
weak decays and can be considered as an important progress in this field, the pioneering LO 
calculations for current-current operators [77], penguin operators [78| and AS = 2 operators 



[79] should not be forgotten. 



2.3.2 NNLO Calculations 

In the case of the CP violating ratio e'/e and the rare decays B — > X s l + l~ and Kl — > 7r°e + e _ , 
the NLO matching conditions for electroweak operators do not involve QCD corrections to box 
and penguin diagrams and consequently the renormalization scale dependence in the top quark 
mass in these processes is not negligible. In order to reduce this unphysical dependence, QCD 
corrections to the relevant box and penguin diagrams have to be computed. In the renormal- 
ization group improved perturbation theory these corrections are a part of next-next-to-leading 



(NNLO) corrections. In [80] and [81] such corrections have been computed for e'/e and the rare 



decays in question, respectively. 

2.3.3 Two-Loop Anomalous Dimensions Beyond the SM 

In the extentions of the SM new operators are present. The two loop anomalous dimensions for 



all AF = 2 four-quark dimension-six operators have been computed in [82, pa\ . In |83|] also the 
corresponding results for AF = 1 can be found. The applications of these results to (AMx,£k) 
and AB = 1 decays in the MSSM can be found in [|84| and ||85[] , respectively. 

2.3.4 Two-Loop Electroweak Corrections 

In order to reduce scheme and scale dependences related to the definition of electroweak pa- 
rameters like sin 2 6w, and uqed, two- loop electroweak contributions to rare decays have to be 
computed. For K® —* ir Q vv, B —* l + l~ and B — > X & vv they can be found in |pq[ , for B° — B° 



mixing in p7j and for B -» X s -/ in fl8& 39, 9C, M. 



2.3.5 NLO Calculations Beyond the SM 

There exist also a number of partial or complete NLO QCD calculations within the Two-Higgs- 
Doublet Model and the MSSM. In the case of the Two-Higgs-Doublet Model such calculations 



for B — B mixing and B — > X s ^ can be found in pi] and [70, 71, 92] respectively. The 



corresponding calculations for B — > X s ^ in the MSSM can be found in |93| and 173]. The latter 
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paper gives also the results for B — ► X s gluon. Finally gluino-mediated NLO-QCD corrections 
to B° — B° mixing in the MSSM have been considered in [94|. 



Table 1: References to NLO Calculations within the SM 



Decay 



Reference 



AF = 1 Decays 



current-current operators 
QCD penguin operators 
electroweak penguin operators 
magnetic penguin operators 
Br(B) SL 

inclusive AS = 1 decays 



, [42], @§ |5| 
, M, pip] 



M 

, pi, m 



Particle-Antiparticle Mixing 



m 



5§ 



Rare K- and U-Meson Decays 



K° -> vr ^, £ l+r, B 



X s vv 



7T + /X/i 



MM 



#L 

B - 
B - 
B - 

Ar Bs 

inclusive B — > Charmonium 

— > D7T, — > 7T7T 



TT°e + e~ 



, 56, 57, 5 



a-0 



H Ell 



2.4 QCD Factorization for Exclusive Non-Leptonic B-Meson Decays 

A simple method for the evaluation of the hadronic matrix elements of four quark operators 
relevant for B decays is the factorization approach in which the matrix elements are expressed 
in terms of products of meson decay constants and formfactors |95| . In its naive formulation, this 
approach gives the /^-independent hadronic matrix elements and consequently //-dependent decay 
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amplitudes. Moreover final state interactions are not taken into account. Various generalizations 
of this method have been proposed in the literature |)(| with the hope to include non-factorizable 
contributions and to remove the //-dependence from the decay amplitudes. Critical reviews of 
these attempts can be found in |)7], |76| . Parallel to these efforts general parametrizations of decay 



amplitudes by means of flavour flow diagrams |9g] and Wick contractions [99, |100(| supplemented 
by dynamical assumptions have been proposed. These parametrizations may turn out to be 
useful when more data will be available. 

Recently factorization for a large class of non-leptonic two-body B-meson decays has been 
shown by Beneke, Buchalla, Neubert and Sachrajda [[Tfj] to follow from QCD in the heavy- 
quark limit. The resulting factorization formula incorporates elements of the naive factorization 
approach but allows to compute systematically non-factorizable corrections. In this approach 
the /i-dependence of hadronic matrix elements is under control. Moreover spectator quark effects 
are taken into account and final state interaction phases can be computed perturbatively. While, 
in my opinion, an important progress in evaluating non-leptonic amplitudes has been made in 



[76], the usefulness of this approach at the quantitative level has still to be demonstrated when 



the data improve. In particular the role of the l/m^ corrections has to be considerably better 



understood. Recent lectures on this approach can be found in [101]. 

There is an alternative perturbative QCD approach to non-leptonic decays fll02fl which has 
been developed earlier from the QCD hard-scattering approach. Some elements of this approach 
are present in the QCD factorization formula of |76|1 . The main difference between these two 
approaches is the treatment of soft spectator contributions which are assumed to be negligible 
in the perturbative QCD approach. While the QCD factorization approach is more general and 
systematic, the perturbative QCD approach is an interesting possibility. Competition is always 
healthy and only time will show which of these two frameworks is more successful and whether 
they have to be replaced by still more powerful approches in the future. 

Finally a new method to calculate the B — > mr hadronic matrix elements from QCD light- 
cone sum rules has been proposed very recently by Khodjamirian pi. This work may shed light 



on the importance of 1/Wb and soft-gluon effects in the QCD factorization approach. Reviews 
of QCD light-cone sum rules can be found in \104]. 



2.5 Inclusive Decays 

So far we have discussed only exclusive decays. It turns out that in the case of inclusive decays 
of heavy mesons, like B-mesons, things turn out to be easier. In an inclusive decay one sums 
over all or over a special class of accessible final states. A well known example is the decay 
B — > Xsj, where X s includes all accessible final states with the net strange quantum number 
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5 = 1. 

At first sight things look as complicated as in the case of exclusive decays. It turns out, 
however, that the resulting branching ratio can be calculated in the expansion in inverse powers 
of TOb with the leading term described by the spectator model in which the B-meson decay is 
modelled by the decay of the 6-quark: 

Br(S -► X) = Br(6 -► q) + £>(— ) . (2.17) 



This formula is known under the name of the Heavy Quark Expansion (HQE) |105| , 106 1. Since 
the leading term in this expansion represents the decay of the quark, it can be calculated in per- 
turbation theory or more correctly in the renormalization group improved perturbation theory. 
It should be emphasized that also here the basic starting point is the effective Hamiltonian Q2.1| ) 
and that the knowledge of Cj(//) is essential for the evaluation of the leading term in ( |2.17| ). 
But there is an important difference relative to the exclusive case: the matrix elements of the 
operators Qi can be "effectively" evaluated in perturbation theory. This means, in particular, 
that their /i and renormalization scheme dependences can be evaluated and the cancellation of 
these dependences by those present in Cj(/i) can be explicitly investigated. 

Clearly in order to complete the evaluation of Br(B — ► X) also the remaining terms in 
( |2,17j ) have to be considered. These terms are of a non-perturbative origin, but fortunately 
they are suppressed by at least two powers of m&. They have been studied by several authors 
in the literature |106| with the result that they affect various branching ratios by less then 
10% and often by only a few percent. Consequently the inclusive decays give generally more 
precise theoretical predictions at present than the exclusive decays. On the other hand their 
measurements are harder. There are of course some important theoretical issues related to the 
validity of HQE in fl2T7| ) which appear in the literature under the name of quark-hadron duality. 



Since these matters are rather involved I will not discuss them here. 
2.6 Penguin— Box Expansion 

The rare and CP violating decays of K and B mesons are governed by various penguin and box 
diagrams with internal top quark and charm quark exchanges. Some examples are shown in 
fig. H. Evaluating these diagrams one finds a set of basic universal (process independent) m t - 



dependent functions F r (xt) [107] where xt = mjr /M§j. Explicit expressions for these functions 
will be given below. 



It is useful to express the OPE formula (2.2) directly in terms of the functions F r {xt) |10 



A(M -> F) = P (M ^F) + J2 P r( M F ) F r(xt), (2-18) 
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u,c,t 



u,c,t 




Figure 4: Typical Penguin and Box Diagrams. 



where the sum runs over all possible functions contributing to a given amplitude. Pq summarizes 
contributions stemming from internal quarks other than the top, in particular the charm quark. 
In the OPE formula ( |2.2| ), the functions F r {xt) are hidden in the initial conditions for Cj(/x) 
represented by C(M\y) in (fOj). 

The coefficients Pq and P r are process dependent and include QCD corrections contained 
in the evolution matrix ?7(/i,M\y). They depend also on hadronic matrix elements of local 
operators and the relevant CKM factors. An efficient and straightforward method for finding 
the coefficients P r is presented in [108]. I would like to call ( |2.18| ) Penguin-Box Expansion 
(PBE). We will encounter many examples of PBE in the course of these lectures. 

Originally PBE was designed to expose the m t -dependence of FCNC processes []108f1 . After 
the top quark mass has been measured precisely this role of PBE is less important. On the other 
hand, PBE is very well suited for the study of the extentions of the SM in which new particles 
are exchanged in the loops. If there are no new local operators the mere change is to modify 
the functions F r (xt) which now acquire the dependence on the masses of new particles such as 
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charged Higgs particles and supersymmetric particles. The process dependent coefficients Pq 
and P r remain unchanged. The effects of new physics can be then transparently seen. However, 
if new effective operators with different Dirac and colour structures are present, new functions 
multiplied by new coefficients P r (M — > F) contribute to ( [2.18 ). 

In the rest of this section we present the functions F r (xt) within the SM. To this end, let us 
denote by Bq, Co and Dq the functions F r (xt) resulting from AF = 1 (F stands for flavour) box 
diagram, Z°-penguin and 7-penguin diagram respectively. These diagrams are gauge dependent 



and it is useful to introduce gauge independent combinations [108 



Xq = Cq — ABq, Y = Co — B , 



Co + -D 



(2.19) 



Then the set of gauge independent basic functions which govern the FCNC processes in the SM 
is given to a very good approximation as follows (x{ = mf/M^): 



So(x t ) = 2.46 
S (x c ,x t 



m t 



Xo{x t ) 
Zo{x t ) = 



1.57 



0.71 



mt 



170 GeV 



In 2 

x c 

1.15 



1.52 



3x t 



170 GeV 



mt 



D' {x t ) = 0.38 



170 GeV 

/ mt 



1.86 



0.60 



4(1 -x*) 
Yo(x t ) 

Eo(xt) -- 



So(x. 
3xt In xt 

= 1.02 

0.26 



Xc 



mt 



E' {x t ) = 0.19 



170 GeV 

mt 
170 GeV 

mt 



1.56 



-1.02 



0.38 



(2.20) 

(2.21) 

(2.22) 
(2.23) 
(2.24) 



V170GeVy ' uv ~ w "" " Vl70GeV, 
The first three functions correspond to AF = 2 box diagrams with (t, t), (c, c) and (t, c) ex- 
changes. Eq results from QCD penguin diagram with off-shell gluon, D' and E' from 7 and 
QCD penguins with on-shell photons and gluons respectively. The subscript "0" indicates that 
these functions do not include QCD corrections to the relevant penguin and box diagrams. 

In the range 150 GeV < mt < 200 GeV these approximations reproduce the exact expres- 
sions to an accuracy better than 1%. These formulae will allow us to exhibit elegantly the mt 
dependence of various branching ratios in the phenomenological sections of these lectures. Exact 
expressions for all functions can be found in JD], 

Generally, several basic functions contribute to a given decay, although decays exist which 
depend only on a single function. We have the following correspondence between the most 
interesting FCNC processes and the basic functions: 
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The super symmetric contributions to the functions ^o, Xq, Yq, Zq and Eq within the MSSM 
with minimal flavour violation (see Section 9) have been recently compiled in |109|. See also 



|llC]-[ |TT3j ]. QCD corrections to these functions can be extracted from papers in table [l] and 



from the section on NLO calculations beyond the SM. In the SM it is convenient in most cases 
to include these corrections in the coefficients P r . Beyond the SM it is better to retain them in 
F r as these corrections depend on the new parameters present in the extentions of the SM. 

3 Particle- Antiparticle Mixing and Various Types 
of CP Violation 

3.1 Preliminaries 

Let us next discuss the formalism of particle-antiparticle mixing and CP violation. Much more 
elaborate discussion can be found in two recent books || [jj. We will concentrate here on K° — K° 
mixing, B® s — B® s mixings and CP violation in K-meson and B-meson decays. As this section 
is rather long it is useful to specify our goals. These are: 

• Presentation of basic concepts of particle-antiparticle mixing and CP violation. 

• Introduction of the CP violating parameters e and e' that describe the so-called indirect 
a direct CP violation in Kl — * irir, respectively. 

• Presentation of a different and more useful classification of different types of CP violation 
that distinguishes between: CP violation in mixing, CP violation in decay and CP violation 
in the interference between mixing and decay. 

• Derivation of a number of formulae that will turn out to be useful in subsequent more 
phenomenological sections. 

It is important to emphasize at this moment that particle-antiparticle mixing and CP vio- 
lation have been of fundamental importance for the construction and testing of the SM. They 
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have also proven often to be undefeatable challenges for suggested extensions of this model. 
In this context the seminal papers of Glashow, Iliopoulos and Maiani [|l^] and of Kobayashi 



and Maskawa [17] should be mentioned. From the calculation of the — K$ mass difference, 
Gaillard and Lee [114] were able to estimate the value of the charm quark mass before charm 



discovery. On the other hand B® — B® mixing 1 1 15 1 gave the first indication of a large top quark 
mass. Next CP violation in the K° - K° mixing offers within the SM a plausible description of 
CP violation in Kl — ► vr7r discovered in 1964 |3^]. Finally the very small values of the measured 
i^L — Ks mass difference and of the CP violating parameter e put severe restrictions on the 
flavour structure and the pattern of complex phases in the extentions of the SM. This is in 
particular the case of general supersymmetric extentions of the SM, in which the mismatch in 
the alignment of the quark mass matrices and the squark mass matrices is very restricted by 
the i^L — Ks mass difference and e. 

It is important to stress that in the SM the phenomena discussed in this section appear 
first at the one-loop level and as such they are sensitive measures of the top quark couplings 
Vti(i = d,s, b) and in particular of the phase (5 = 7. They allow then to construct the unitarity 
triangle as explicitly demonstrated in Section 4. 



Let us next enter some details. The following subsection borrows a lot from [116, 117|. The 
discussion of different types of CP violation benefited from several very nice lectures by Nir 



[115], although the presentation of this topic below differs occassionally from his. 
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Figure 5: Box diagrams contributing to K° — K° mixing in the SM. 



3.2 Express Review of K° — K° Mixing 

K° = (sd) and K° = (sd) are flavour eigenstates which in the SM may mix via weak interactions 
through the box diagrams in fig. [|. We will choose the phase conventions so that 

CP\K°) = -\K°), CP\K°) = -\K°). (3.1) 
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In the absence of mixing the time evolution of \K°(t)) is given by 

\K°(t)) = \K°(Q)) exp(-ifft) , H = M - i- , (3.2) 

where M is the mass and F the width of K . Similar formula exists for K°. 

On the other hand, in the presence of flavour mixing the time evolution of the K — K° 
system is described by 



where 



.f / Mn -i^ M12 -i^ 



H = M-i-=\ ' ' (3.4) 

- \ M 21 -i\ M 22 -i^f 

with M and T being hermitian matrices having positive (real) eigenvalues in analogy with M 
and r. Mij and Tij are the transition matrix elements from virtual and physical intermediate 
states respectively. Using 

M 21 = M{ 2 , r 2i = T* 12 , (hermiticity) (3.5) 

Mn = M 22 = M , Tn = T 22 = T , (CPT) (3.6) 



we have 



M -q Ma -j 



Diagonalizing ( |3.3| ) we find: 
Eigenstates: 



(1 + e)K° ± (1 - e)K° 
= V2(1+ | £ - p) ^ 3 - 8 ) 



where e is a small complex parameter given by 

1 - £ 



M 12 ~ 4 F 12 _ AM ~ 4 Ar _ 2M 12 - ir !2 



rexp(m) . (3-9) 



1 + e \ Mi 2 - i|ri 2 2Mi 2 - tTia AM - iiAT 

with Ar and AM given below. 
Eigenvalues: 

M Li5 = M±ReQ r Li5 = r^2ImQ (3.10) 

where 

Q = \j{M 12 -i l -T 12 ){Ml 2 -i l -T\ 2 ). (3.11) 

Consequently we have 

AM = Ml - Ms = 2ReQ , Af = T L - T s = -4ImQ. (3.12) 
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It should be noted that the mass eigenstates K$ and Kl differ from CP eigenstates 

Ki = -^(K° - K°), CP\Ki) = \K X ) , (3.13) 

K 2 = -±=(K° + K°), CP\K 2 ) = -\K 2 ) , (3.14) 

by a small admixture of the other CP eigenstate: 

K x + eK 2 K 2 + eK x 

Ks = vmw Kh = vmw {3A5) 

Since e is O(10 -3 ), one has to a very good approximation: 

AM K = 2ReM 12 , Ar^ = 2R e r 12 , (3.16) 

where we have introduced the subscript K to stress that these formulae apply only to the K°—K° 
system. 



The i^L — Ks mass difference is experimentally measured to be [24 



AM K = M(K L ) - M(K S ) = (3.489 ± 0.008) • 10~ 15 GeV . (3.17) 

In the SM roughly 70% of the measured AMk is described by the real parts of the box diagrams 
with charm quark and top quark exchanges, whereby the contribution of the charm exchanges is 
by far dominant. This is related to the smallness of the real parts of the CKM top quark couplings 
compared with the corresponding charm quark couplings. Some non-negligible contribution 
comes from the box diagrams with simultaneous charm and top exchanges. The remaining 20% 
of the measured AMr is attributed to long distance contributions which are difficult to estimate 
[HE]- Further information with the relevant references can be found in pip. 



The situation with Ar^ is rather different. It is fully dominated by long distance effects. 
Experimentally one has Ar^ ~ — 2AM^. We will use this relation in what follows. 

Generally to observe CP violation one needs an interference between various amplitudes that 
carry complex phases. As these phases are obviously convention dependent, the CP-violating 
effects depend only on the differences of these phases. In this context it should be stressed that 
the small parameter e depends on the phase convention chosen for K° and K°. Therefore it 
may not be taken as a physical measure of CP violation. On the other hand Re e and r, defined 
in ( p.Sj ) are independent of phase conventions. In particular the departure of r from 1 measures 
CP violation in the K° - K° mixing: 

— — -.mi l 1 ~ I + mi l l 

4|JV/ 12 |2 + |r 12 |2 Vr 12 ; \r u J 



! + JEhLm (Mr) „ j + , m (Mr) . (3M) 
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This type of CP violation can be best isolated in semi-leptonic decays of the Kl meson. The 
non-vanishing asymmetry 

OSL(Kl) _ ! e 3> - (3.i9) 



F(K L -> Tr-e+iy e ) + T(K L -► ^+ e -p e ) V AW A' 
signals this type of CP violation. Equivalently 

asi(K L ) = ^2 = 2Ree (3.20) 

1 + T z 

Note that asx (Kl) is determined purely by the quantities related to K°—K° mixing. Specifically, 
it measures the difference between the phases of Ti2 and M\2- 

That a non-vanishing osl(-Kl) is indeed a signal of CP violation can also be understood in 
the following manner. Kl, that should be a CP eigenstate K2 in the case of CP conservation, 
decays into CP conjugate final states with different rates. As Ree > 0, Kl prefers slightly to 
decay into Tr~e + v e than 7r + e~v e . This would not be possible in a CP-conserving world. 

3.3 The First Look at e and e' 

Since a two pion final state is CP even while a three pion final state is CP odd, K$ and Ki 
preferably decay to 2n and 37r, respectively via the following CP-conserving decay modes |120[|: 



K L -» 3vr (via K 2 ), K s -► 2vr (via K x ). (3.21) 

This difference is responsible for the large disparity in their life-times. A factor of 579. However, 
Ki, and K§ are not CP eigenstates and may decay with small branching fractions as follows: 

Ki 2vr (via Ki), K s -» 3vr (via K 2 ). (3.22) 

This violation of CP is called indirect as it proceeds not via explicit breaking of the CP symmetry 
in the decay itself but via the admixture of the CP state with opposite CP parity to the dominant 
one. The measure for this indirect CP violation is defined as (I=isospin) 

A(K L -> (7T7r) 7=0 ) 



A(K S -> (7T7T) /= 



(3.23) 



Following the derivation in [116] one finds 



-^^W 12 + ^M 12 >, f=^. ,3,4, 

The phase convention dependence of the term involving ^ cancells the convention dependence 
of e so that e is free from this dependence. The isospin amplitude Aq is defined below. 
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The important point in the definition ( |3.23| ) is that only the transition to (vr7r)/ = o enters. 
The transition to (tttt)j = 2 is absent. This allows to remove a certain type of CP violation that 
originates in decays only. Yet as e ^ e and only Ree = Ree, it is clear that e includes a type 
of CP violation represented by Ime which is absent in the semileptonic asymmetry ( |3.19| ). We 
will identify this type of CP violation in Section 3.7, where a more systematic classification of 
different types of CP violation will be given. 
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Figure 6: Indirect versus direct CP violation in Kl —* mr. 



While indirect CP violation reflects the fact that the mass eigenstates are not CP eigenstates, 
so-called direct CP violation is realized via a direct transition of a CP odd to a CP even state 
or vice versa (see fig. H). A measure of such a direct CP violation in Kl — ► tttt is characterized 
by a complex parameter e' defined as 



1 A 



2.L 



A2,S A),L 



(3.25) 



y/2 \A ,s A 0tS A 0i sJ 

where short hand notation Aj l = A{Kl — ► (inr)i) and Ai s = A(K$ — * (irir)i) has been used. 

This time the transitions to (7T7t)/ = o and (7T7t)/ = 2 are included which allows to study CP 
violation in the decay itself. We will discuss this issue in general terms in Section 3.7. For the 



time being it is useful to cast ( 3.25 ) into a more transparent formula 



i=Im(^) ( ,M)(,<I> 



7T 



So, 



where the isospin amplitudes Aj in K — » 7T7t decays are introduced through 



A(K+ -» TT+TT ) 



-A 2 e i52 



A(K° -► vt+tt-) = \ -A e i5 ° + 



A(K C 



7T 7T 



A n e iSo 



2\i —j±o 
3 ' 



A 2 e iS * 



A 2 e iS * 



(3.26) 

(3.27) 

(3.28) 
(3.29) 



Here the subscript I = 0, 2 denotes states with isospin 0, 2 equivalent to AI =1/2 and AI = 3/2 
transitions, respectively, and 6q 2 are the corresponding strong phases. The weak CKM phases 
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are contained in Aq and A 2 . The isospin amplitudes Aj are complex quantities which depend 
on phase conventions. On the other hand, e' measures the difference between the phases of A2 
and Aq an d is a physical quantity. The strong phases Jo, 2 can be extracted from tttt scattering. 
Then $ £ , tt/4. 

Experimentally e and e' can be found by measuring the ratios 

7T 7r°) 



??00 



A(K L 



A(K S -> vr°7r ) ' V+ ~ A(K S 
Indeed, assuming e and e' to be small numbers one finds 

2s' 



A(K L 



TT 1 7T 



7M7T 



Voo 



1- V2< 



uJ 



V+- 



e + 



(3.30) 



(3.31) 



where u = ReA 2 /ReA = 0.045. 

In the absence of direct CP violation 7700 
through 

Re(s'/s] 



The ratio s'/s can then be measured 



??oo 



_ (3 32) 

6(l + o;/V2) 1/ ??+- y 

To my knowledge J121| ] the experimental groups in giving their results for Re(e'/e) omitt the 

term in (|3.32j ). Yet in order to be consistent with the definitions ( 3.23 ) and ( |3.25| ) used by 

theorists this term should be kept. Consequently the existing experimental results for Re(e'/e) 

quoted below should be rescaled down by 3.2%. Clearly at present this rescaling is academic as 

the experimental error in Re(e'/e) is roughly ±20% and the theoretical one at least ±50%. We 

will therefore omitt this rescaling in what follows. 



3.4 Basic Formula for e 

With all this information at hand let us derive a formula for e which can be efficiently used 
in pheneomenological applications. The off-diagonal element M12 in the neutral ET-meson mass 
matrix representing K°-K° mixing is given by 

2m K M{ 2 = (K°\n e{i (AS = 2)\K°) , (3.33) 

where H. e f[(AS = 2) is the effective Hamiltonian for the AS = 2 transitions. That M* 2 and 
not M12 stands on the l.h.s of this formula, is evident from fl3.7| ). The factor 2rriK reflects our 
normalization of external states. 

To lowest order in electroweak interactions AS = 2 transitions are induced through the box 
diagrams of fig. |5[ Including leading and next-to-leading QCD corrections in the renormalization 
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group improved perturbation theory one has for fi < \i c = 0(m c 



n_,AS=2 



167T 2 



F M 2 



X c r]iSo(x c ) + \ t r} 2 So(x t ) + 2\ c \ t r]3So(x c , x t ) 



-2/9 



Q(AS = 2) + h.c. 



(3.34) 



(3) 

where Aj = V* s V id , a s is the strong coupling constant in an effective three flavour theory and 
J 3 = 1.895 in the NDR scheme |§. In fl3.34| ), the relevant operator 



Q(A5 = 2) = (fd) v _ A (ld) v _ A , 



(3.35) 



is multiplied by the corresponding Wilson coefficient function. This function is decomposed into 
a charm-, a top- and a mixed charm-top contribution. The functions So are given in ( 2.20| ) and 

(H). 

Short-distance QCD effects are described through the correction factors rji, 7] 2 , r/3 and the 
explicitly (^-dependent terms in (3.34). The NLO values of rji are given as follows [51, 52. l54j: 



r/i = 1.38 ± 0.20, 



r/ 2 = 0.57 ± 0.01, % = 0.47 ± 0.04 . 



(3.36) 



The quoted errors reflect the remaining theoretical uncertainties due to leftover //-dependences 
at Ofa 2 .) and A-jjg, the scale in the QCD running coupling. 

Defining the renormalization group invariant parameter Bk by 

ai 3) (/i) 



B K = B K (fi) afV) 



-2/9 



1 + 



47T 



-J 3 



(K°\(sd)v-A(sd) V - A \K ) = -B K (n)F K m 2 K 



and using ( 3.34| ) one finds 

<2 



M 



G 2 

12 = -^^F K B K m K M^j X*cT]iSq{x c ) + \fr) 2 So(xt) + 2A*A t *r/ 3 S , (x c , 



(3.37) 



(3.38) 



(3.39) 



where = 160 MeV is the iT-meson decay constant and mx the X-meson mass. It should 
be mentioned that in the normalization of external states in which the factor 2rriK in ( |3.33 ) is 
absent, the r.h.s of ( 3.38j ) is divided by 2m k so that ( 3.39j ) remains unchanged. 



To proceed further we neglect the last term in ( 3.24 ) as it constitutes at most a 2 % correction 
to e. This is justified in view of other uncertainties, in particular those connected with Bk- 
Inserting flOg ) into (^2§) we find 



e = C e B K lm.\ t {ReA c [r/i5 (x c ) - rj 3 S (x c , x t )] - ReA i r/ 2 5o(x t )} exp(i7r/4) 



(3.40) 
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where we have used the unitarity relation ImA* = ImAt and have neglected ReAj/ReA c = 0{\ 
in evaluating Im(A*A|). The numerical constant C £ is given by 

i2 t?2 „, ji#2 



= G F F K m J ,M w = 4 



To this end we have used the experimental value of AMk in ( 3.17j ) and M\y = 80.4 GeV. 



Using the standard parametrization of (1.5) to evaluate ImAj and ReAj, setting the values 



for S12, sis, S23 and mt in accordance with experiment and taking a value for Bk (see below), 
one can determine the phase S by comparing (|3.40|) with the experimental value for e 



e exp = (2.280 ± 0.013) • 10~ 3 exp i$ £ , $ e = ^. (3.42) 



Once S has been determined in this manner one can find the apex (g, if) of the unitarity 
triangle in fig. [l] by using 

3 cos 5, i] = Sl3 sin 5 (3.43) 



•S12S23 S12S23 
and 2 2 

e = e(l-y), ^ = »?(l-y)- (3-44) 

For a given set (s±2, S13, S23, mt, Bjr) there are two solutions for 5 and consequently two 
solutions for (g,fj). This will be evident from the analysis of the unitarity triangle discussed in 
detail below. 

Finally we have to say a few words about the non-perturbative parameter Bk, the main 
uncertainty in this analysis. References to older estimates can be found in @. In our numerical 
analysis presented below we will use 

5^ = 0.85 ±0.15 (3.45) 



which is very close to most recent lattice estimates, reviewed recently in [122, |123j], slightly 



higher that the large-N estimates [124, 125, 126] and somewhat lower than chiral quark model 
estimates [1127 



Recently an interesting large-N calculation of Bk in the chiral limit and including next-to- 



leading 1/N corrections has been performed by Peris and de Rafael [125]. The nice feature of 
this calculation is the explicit cancellation of the ^/-dependence and the renormalization scheme 
dependence between the Wilson coefficient and the matrix element of the operator Q(AS = 2). 
The resulting Bk = 0.41 ± 0.09 is by a factor of two lower than the value in ( |3.45| ) and the 
lattice results. However, before a meaningful comparision with the lattice values can be made, 



higher order corrections in the chiral expansion have to be added to the result in [ 128 
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3.5 Express Review of B ds -B ds Mixing 
The flavour eigenstates in this case are 

B° d = (bd), B° d = (bd), B° s = (bs), B° s = (bs). (3.46) 

They mix via the box diagrams in fig. || with s replaced by b in the case of B d -B d mixing. In 
the case B®-B® mixing also d has to be replaced by s. 

Dropping the subscripts (d, s) for a moment, it is customary to denote the mass eigenstates 

by 

B H = pB° + qB°, B L = pB° - qB° (3.47) 

where 

p= vwrwm' "'vwrnrr (3 ' 48) 

with sb corresponding to e in the K° — K° system. Here "H" and "L" denote Heavy and Light 
respectively. As in the B° — B° system one has Ar <C AM, it is more suitable to distinguish 
the mass eigenstates by their masses than the corresponding life-times. 

The strength of the B® — B® mixings is described by the mass differences 

AM dtS = M^ s - M^ s . (3.49) 

In contrast to AMx , in this case the long distance contributions are estimated to be very small 
and AM(j )S is very well approximated by the relevant box diagrams. Moreover, due m U)C <C mt 
only the top sector can contribute significantly to B ds — B ds mixings. The charm sector and 
the mixed top-charm contributions are entirely negligible. 

AM,j 8 can be expressed in terms of the off-diagonal element in the neutral B-meson mass 
matrix by using the formulae developed previously for the K-meson system. One finds 

AM q = 2\M[f\, q = d,s. (3.50) 

This formula differs from AMk = 2ReMi2 because in the B-system Ti2 <C M±2- 
We also have 

AB = T(B H ) - T(B L ) = 2 Re( ^ 2) (3.51) 



and 



q _ 2M X * 2 - iT\ 2 _ M 



12 



1 Im 



2 \M 



12 



'12 



P AM-i\AT \M 12 \ 
where higher order terms in the small quantity T12/M12 have been neglected. 
The smallness of Im(ri2/Mi2) < O(10 -3 ) has two important consequences: 



(3.52) 
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• The semileptonic asymmetry osl(-B) discussed a few pages below is even smaller than 
osl(-Kl)- Typically 0(1CT 4 ). These are bad news. 

• The ratio q/p is a pure phase to an excellent approximation. These are very good news as 
we will see below. 

Inspecting the relevant box diagrams we find 



(M* 2 )d oc (ytdVtty , {M* l2 ) s (x (VtsVtt 



*\2 



(3.53) 



Now, from Section 1 we know that 



Vtd = Pule-*, 



Vt 



Is 



(3.54) 



with j3 s = O(10 ). Consequently to an excellent approximation 



\P/d,s 



with (fyj^ given entirely by the weak phases in the CKM matrix. 



(3.55) 



3.6 Basic Formulae for AMd }S 

Let us next find AM^ S . The off-diagonal term M\i in the neutral U-meson mass matrix is given 
by a formula analogous to fljT33D 



2m Bq \M[f\ = \(B° q \H eS (AB = 2)\B° q ) 



(3.56) 



where in the case of B^ — B^ mixing 



n _rAB=2 



16vr 2 ' 



-6/23 



1 + 



Q(AB = 2)+h.c 



Here /i& = 0(mb), J5 = 1.627, 



(3.57) 



Q(AB = 2) = (bd) V - A (bd) V -A 



(3.58) 
(3.59) 



and H, [53 1 

r7B =0.55 ±0.01. 

In the case of B® — B® mixing one should simply replace d — > s in (3.57) and (3.58) with all 
other quantities and numerical values unchanged. 
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Defining the renormalization group invariant parameters B q in analogy to ( 3.37| ) and ( |3.3£ 
one finds using (|3,57| ) 



AM, = ^VBm Bq (B Bq F^)M4,S (x t )\V tq \ 



where Fb is the fL-meson decay constant. This implies two useful formulae 



and 



AM d = 0.50/ps • 



AM S = 15.1/ps 



B B d FB d 



200 MeV 



m t {m t ) 
170 GeV 



1.52 



td 



8.8 • 10" 



VB 



0.55 



B B F Ba 



240 MeV 



-/„ \ n 1.52 



m t (m t 



170 GeV 



\v t 



n2 r 



0.040 



riB 



0.55 



(3.60) 



(3.61) 



(3.62) 



There is a vast literature on the calculations of F B . and B^s- The most recent lattice 



results are summarized in [129, |123f| . They are compatible with the results obtained with the 
help of QCD sum rules fll30ll . Guided by [l2| we will use in our numerical analysis the 



value for Fs d y Bb a given in table [|. The experimental situation on AM^ S is also given there. 



3.7 Classification of CP Violation 
3.7.1 Preliminaries 

We have mentioned in Section 2 that due to the presence of hadronic matrix elements, various 
decay amplitudes contain large theoretical uncertainties. It is of interest to investigate which 
measurements of CP-violating effects do not suffer from hadronic uncertainties. To this end it is 
useful to make a classification of CP-violating effects that is more transparent than the division 
into the indirect and direct CP violation considered so far. To my knowledge this classification 
has been developed first for B decays but it can also be useful for K decays. A nice detailed 
presentation can be found in [ 131 j . 

Generally complex phases may enter particle-antiparticle mixing and the decay process itself. 
It is then natural to consider three types of CP violation: 

• CP Violation in Mixing 

• CP Violation in Decay 

• CP Violation in the Interference of Mixing and Decay 

As the phases in mixing and decay are convention dependent, the CP-violating effects depend 
only on the differences of these phases. This is clearly seen in the classification given below. 
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3.7.2 CP Violation in Mixing 

This type of CP violation can be best isolated in semi-leptonic decays of neutral B and K mesons. 
We have discussed the asymmetry clsl{Kl) before. In the case of B decays the non-vanishing 
asymmetry (we suppress the indices (d,s)) 

, m r(B (t)^i + vX)-r(B (t)^rpx) 1 - \g/ P \ 4 ( r 12 \ 

signals this type of CP violation. Here B°(0) = B°, B°(0) = B°. For t / the formulae 
analogous to Q3.3| ) should be used. Note that the final states in ( |3.63| ) contain "wrong charge" 
leptons and can only be reached in the presence of B — B° mixing. That is one studies effectively 
the difference between the rates for B° —* B° — ► l + vX and B° — ► B° —* l~vX. As the phases in 
the transitions B° — ► B° and B° — > B° differ from each other, a non-vanishing CP asymmetry 
follows. Specifically ash(B) measures the difference between the phases of T12 and M\i- 

As Mi 2 and in particular Yyi suffer from large hadronic uncertainties no precise extraction 
of CP-violating phases from this type of CP violation can be expected. Moreover as q/p is 
almost a pure phase, see (|3.52 ) and ( 3.55| ), the asymmetry is very small and outside the reach 



of experiments performed in the coming years. 
3.7.3 CP Violation in Decay 

This type of CP violation is best isolated in charged B and charged K decays as mixing effects 
do not enter here. However, it can also be measured in neutral B and K decays. The relevant 
asymmetry is given by 

decay _ T(B+ -> /+) - T(B~ -> f~) _ 1 - \A f - /A f+ \* 

a /± -r(B+-/+) + r(B--/-) " i + \A f -/A f+ \i 

where 

A f+ = {f + \H wcak \B + ), A f - = {f-\7i wcak \B-) . (3.65) 

For this asymmetry to be non-zero one needs at least two different contributions with different 
weak (<pi) and strong (Si) phases. These could be for instance two tree diagrams, two penguin 
diagrams or one tree and one penguin. Indeed writing the decay amplitude Aj+ and its CP 
conjugate Aj- as 

A f+ = =A ie i(s ^\ A f - = 2 =A i e i ^~^\ (3.66) 

i=l,2 i=l,2 

with Ai being real, one finds 



decay _ -2A X A 2 sm(5i - 5 2 ) sin(^i - 2 ) 
lf± " A\ + A 2 2 + 2A 1 A 2 cos(<5i - S 2 ) cos(^>i - fa) 



(3.67) 
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The sign of strong phases <5j is the same for Aj+ and Aj- because CP is conserved by strong 
interactions. The corresponding weak phases have opposite sign. 

The presence of hadronic uncertainties in Ai and the presence of strong phases Si complicates 
the extraction of the weak phases <fn from data. An example of this type of CP violation in K 
decays is e'. We will demonstrate this below. 



3.7.4 CP Violation in the Interference of Mixing and Decay 

This type of CP violation is only possible in neutral B and K decays. We will use B decays 
for illustration suppressing the subscripts d and s. Moreover, we set Ar = 0. Formulae with 
Ar / can be found in @j. 

Most interesting are the decays into final states which are CP-eigenstates. Then a time 
dependent asymmetry defined by 

acp{t ' f) - rW^PFT) ( } 

is given by 

a CP (t, f) = Acp y (B -» /) cos(AMt) + jtig P {B -> /) sin(AMt) (3.69) 

where we have separated the decay CP-violating contributions from those describing CP violation 
in the interference of mixing and decay: 

A^(B - /) = A { Sp(B - /) = 2ImC/ 2 . (3.70) 



The later type of CP violation is sometimes called the mixing-induced CP violation |132(| . The 
quantity £ f containing essentially all the information needed to evaluate the asymmetries ( |3.70| ) 
is given by 

= p A(BO^f) = WW") A{B o ^ f) ( 3 - 71 ) 
with 4>m, introduced in ( |3.55| ), denoting the weak phase in the B° — B° mixing. A(B° — > /) and 
A(B° — * f) are decay amplitudes. The time dependence of acp(t,f) allows to extract A^p 3,7 
and Aq P as coefficients of cos(AMt) and sin(AMt) respectively. 

Generally several decay mechanisms with different weak and strong phases can contribute to 
A(B° — » /). These are tree diagram (current-current) contributions, QCD penguin contributions 
and electroweak penguin contributions. If they contribute with similar strength to a given decay 
amplitude the resulting CP asymmetries suffer from hadronic uncertainies related to matrix 
elements of the relevant operators Qi. The situation is then analogous to the class just discussed. 
Indeed 

"^4 T e*( 5 T-<^T) _|_ A pe i(6p-ct>p)~ 



A(B° -> /) 
A(B° -+ f) Vf 



A T e i ^ T+< t >T ^> + Ape^p^p^ 



(3.72) 
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with rjf = ±1 being the CP-parity of the final state, depends on strong phases 5t,p and hadronic 
matrix elements present in At. p. Thus the measurement of the asymmetry does not allow a 
clean determination of the weak phases 0t,p- The minus sign in ( 3.72] ) follows from our CP 



phase convention CP\B } = — \B ), that has also been used in writing the phase factor in ( |3.71| ). 
Only £ is phase convention independent. Explicit derivation can be found in section 8.4.1 of 0]. 

An interesting case arises when a single mechanism dominates the decay amplitude or the 
contributing mechanisms have the same weak phases. Then the hadronic matrix elements and 
strong phases drop out and 

is a pure phase with tp£> being the weak phase in the decay amplitude. Consequently 

£/ = -?7/exp(i20 M )exp(-i20£)), | | 2 = 1 . (3.74) 

In this particular case A < ^^ >,J (B — ► /) vanishes and the CP asymmetry is given entirely in terms 
of the weak phases m and 4>d- 

a C p(t, f) = Im£/ sin(AMt) Im£/ = rjf sin(20 D - 20 M ) • (3.75) 

Thus the corresponding measurement of weak phases is free from hadronic uncertainties. A well 
known example is the decay Bd — > tpKs- Here 0m = — and 0£> = 0. As in this case rjf = — 1, 
we find 

a C p(t, f) = - sin(2/3) sin(AMt) (3.76) 

which allows a very clean measurement of the angle /3 in the unitarity triangle. We will return 
to this decay and other decays in which this type of CP violation can be tested. 

We observe that the asymmetry acp(t, f) measures directly the difference between the phases 
of B° — l? -mixing (24>m) and of the decay amplitude (24>d)- This tells us immediately that 
we are dealing with the interference of mixing and decay. As 4>m and 4>d are obviously phase 
convention dependent quantities, only their difference is physical, it is impossible to state on 
the basis of a single asymmetry whether CP violation takes place in the decay or in the mixing. 
To this end at least two asymmetries for B^(B^) decays to different final states fa have to be 
measured. As 0m does not depend on the final state, Im^ / Im£/ 2 is a signal of CP violation 
in the decay. The same applies to B®(B®) decays. 

In the case of K decays, this type of CP violation can be cleanly measured in the rare decay 
Kl — ► tt^vv. Here the difference between the weak phase in the K° — K° mixing and in the 
decay s — ► dvv matters. 

We can now compare the two classifications of different types of CP violation. CP violation 
in mixing is a manifestation of indirect CP violation. CP violation in decay is a manifestation 



37 



of direct CP violation. CP violation in interference of mixing and decay contains elements of 
both the indirect and direct CP violation. 

It is clear from this discussion that only in the case of the third type of CP violation there are 
possibilities to measure weak phases without hadronic uncertainties. This takes place provided 
a single mechanism (diagram) is responsible for the decay or the contributing decay mechanisms 
have the same weak phases. 

3.7.5 Another Look at e and e' 

Let us finally investigate what type of CP violation is represented by e and e' . Here instead of 
different mechanism it is sufficient to talk about different isospin amplitudes. 

In the case of e, CP violation in decay is not possible as only the isospin amplitude Aq is 



involved. See ( |3.23| ). We know also that only Ree = Ree is related to CP violation in mixing. 
Consequently: 

• Ree represents CP violation in mixing, 

• Ime represents CP violation in the interfernce of mixing and decay. 



In order to analyze the case of e' we use the formula ( 3.26| ) to find 

A 2 



Ree' = — -= 
V2 



Ime' = — = 
V2 



A 
A 2 



A 



sm(4> 2 - 4> ) sin(<5 2 - <5o) (3.77) 
sin (02 — 4>o) cos(S 2 — So) ■ (3.78) 



Consequently: 



• Re e' represents CP violation in decay as it is only non zero provided simultaneously 
<fi 2 ¥= 4>o and S 2 / S - 

• Im e' exists even for S 2 = So but as it requires 4> 2 ^ (fio it represents CP violation in decay 
as well. 

Experimentally S 2 ^ So. Within the SM, <p 2 and </>o are connected with electroweak penguins 
and QCD penguins respectively. Do these phases differ from each other so that a nonvanishing 
e' is obtained? We will return to this question in Section 5. 



4 Standard Analysis of the Unitarity Triangle 
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4.1 Basic Procedure 



With all these formulae at hand we can now summarize the standard analysis of the unitarity 
triangle in fig. |l[ It proceeds in five steps. 
Step 1: 

From b — > c transition in inclusive and exclusive leading B-meson decays one finds \V cb \ and 
consequently the scale of the unitarity triangle: 



\V cb \ 



X\V cb \ = X 3 A 



(4.79) 



Step 2: 

From b — > u transition in inclusive and exclusive B meson decays one finds IK^/V^I and 



consequently using ( 1.32 ) the side CA = R b of the unitarity triangle: 



Vub 



V cb 



Ri, 



+ n 2 = 4.44 • 



v ub 



v cb 



(4.80) 



Step 3: 

From the experimental value of e ( |3.42| ) and the formula ( |3.40| ) one derives, using the ap- 

the constraint 

"(1 - g)A 2 m S (x t ) + P c (e)] A 2 B K = 0.226, (4.81) 



proximations (1.22)— (1.24), the constraint 



where 



P c (e) = [r) 3 S (x c ,x t ) 



x t 



m t 



(4.82) 



P c (e) = 0.31 ± 0.05 summarizes the contributions of box diagrams with two charm quark ex- 
changes and the mixed charm-top exchanges. The main uncertainties in the constraint ( |4.81| ) 
reside in Bk and to some extent in A^ which multiplies the leading term. Equation ( |4,81| ) 
specifies a hyperbola in the (g, fj) plane. This hyperbola intersects the circle found in step 2 in 
two points which correspond to the two solutions for 5 mentioned earlier. This is illustrated in 
fig. ^. The position of the hyperbola ( f4.81| ) in the (g, fj) plane depends on m t , \V cb \ = AX 2 and 
Bk- With decreasing m t , \ V cb \ and Bk the e-hyperbola moves away from the origin of the (g, f)) 
plane. 

Step 4: 

From the observed B^ — B^ mixing parametrized by AM^ the side BA = Rt of the unitarity 
triangle can be determined: 



Rt = x- 



l Wt 



td\ 



X \V r 



1.0 



cb\ 



" \Vtd\ ' 
8.8- 10-3 



"0.040" 



W, 



cb\ 



(4.83) 
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Figure 7: Schematic determination of the Unitarity Triangle. 



with 



|^| = 8.8-10- 



200 MeV 



B B d F B , 



170 GeV 



-I 0.76 r 



m t (m t ) 



AMd 



-|0.5 



0.50/ps 



/0.55 
f]B 



(4.84) 



Since m t , AMd and rjs are already rather precisely known, the main uncertainty in the 
determination of | Vtd | from — B® mixing comes from Fg d J B Bd . Note that Rt suffers from 
additional uncertainty in |V^,|, which is absent in the determination of \Vtd\ this way. The 
constraint in the (g, fj) plane coming from this step is illustrated in fig. |7[ 

Step 5: 

The measurement of B® — B® mixing parametrized by AM S together with AMd allows to 
determine Rt in a different manner. Using ( |3.60| ) one finds 



\Vu\ 



\V t 



ts 




Fb s \Bb. 



(4.85) 



Now to an excellent accuracy [27]: 
\V td \ = \Vcb\XRt, 



\v t 



ts 



\Vcb\(l 



F Bd JB Bt 



-A 2 + ^A 2 ) 



(4.86) 



We note next that through the unitarity of the CKM matrix, the present experimental upper 
bound on AMd/ AM S (see table ||) and the value of |V^b/V c fc| one has < g < 0.5, where 
^ = 1.15 ± 0.06 fl2$ , |12S|1 has been used. Consequently \Vt s \ deviates from \V c b\ by at most 2%. 
This means that to a very good accuracy we can set \Vt s \ = \V c b\. Consequently (4.85) and the 
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first formula in (4.86) imply 



Rt = 0.84 £ 



cff1 



AM d 



Using next AM d m 



0.50/ps' 

0.50/ps one finds a useful formula 



£cff = 



15.0/ps 



(4.87) 



(A, 



t J max 



0.84 C\ 



15.0/ps 



(4. 



(AM s ) min ' 

If necessary the C(A 2 ) corrections in ( 4.86j) can be incorporated in ( 4.87 ). This will be only 
required when the error on £ will be decreased below 2%, which is clearly a very difficult task. 

One should note that mt and |V^,| dependences have been eliminated this way and that 
£ should in principle contain much smaller theoretical uncertainties than the hadronic matrix 



elements in AM^ and AM S separately. The most recent values relevant for ( 4.88j ) are summarized 
in table 0. 



Table 2: The ranges of the input parameters. 



Quantity 


Central 


Error 


\Vcb\ 


0.041 


±0.002 


\V u b/V c b\ 


0.085 


±0.018 


\v ub \ 


0.00349 


±0.00076 


B K 


0.85 


±0.15 


\[B~ d F Bd 


230 MeV 


±40 MeV 


m t 


166 GeV 


±5 GeV 


(AM) d 


0.487/ps 


±0.014/ps 


(AM) S 


> 15.0/ps 






1.15 


±0.06 



4.2 Numerical Results 
4.2.1 Input Parameters 

The input parameters needed to perform the standard analysis of the unitarity triangle are 
given in table ||. In constructing this table I was guided to a large extent by the reviews 
122, 122, 129, 133, 134]. I am aware of the fact that other authors would possibly use slightly 



different ranges for input parameters. Still table || is representative for the present situation. 
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Please note, however, that the error on \V u b\ in table g is by a factor of two larger than in |]123]| . 
I do not think that our present understanding of theoretical uncertainties in the determination 
of \V u b\ is sufficiently good that an error of ±10% on this element can be defended. 

The great progress during the last year has been the improved lower limit on AM S from 



LEP and SLD as reviewed in [133]. The value of mt refers to the running current top quark 
mass defined at [i = m^ ole . It corresponds to mf ole = 174.3 ± 5.1 GeV measured by CDF and 
DO (135 1 . 




Figure 8: The Unitarity Triangle as of January 2001. 



4.2.2 Various Error Analyses 



Having set the input parameters and their uncertainties there is the question how to treat the 
theoretical uncertainties in a quantitative analysis. There is a hot discussion on this issue, that 
can be traced by reading [p], 122, 132, 134, |136|, 137, |138||. Basically three different approaches 



can be found in the literature: 

• Gaussian Method: The experimentally measured numbers and the theoretical input pa- 



rameters are used with Gaussian errors. Examples are the analyses in [123, 137, |139| 1. 
There are some small differences in the actual treatment of errors in these papers. We 



refer to [123] for a detailed presentation. 



• BaBar 95% C.L. Scanning Method: This method has been developed in [136] and is the 
official method of the BaBar collaboration ||. In this method one sets the theoretical 
input parameters at some fixed values and finds the allowed (95% C.L.) region for (g,fj) 
by using gaussian errors for the experimental input parameters. Repeating this procedure 
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for different sets of the values of the theoretical input parameters one obtains an envelope 
of 95% C.L. regions. The latest application of this method can be found in |138[| . 



• Simple Scanning: Both the experimentally measured numbers and the theoretical input 
parameters are scanned independently within ranges, given for instance in table This 
is the method used for instance by Rosner [|i]] and by myself below. 

In my opinion the use of Gaussian errors for theoretical input parameters is questionable 
but I do not want to enter this discussion here. A recent attempt to justify this method can 



be found in [122, 133 1. Yet, when the lattice calculations improve dramatically I could imagine 
that one could defend this approach. On the other hand the simple scanning method appears 
to be too conservative. It should be stressed that only the gaussian method pretends to give 
standard deviations for the output quantities. The scanning methods can only give ranges for 
the output quantities. The BaBar scanning gives generally the ranges for quantities of interest 
which are comparable to the ones found by the more naive scanning used here. The 95%C.L. 
ranges from the gaussian method are not so different from the ones obtained by the scanning 
methods and consequently global pictures of the unitarity triangle obtained by these methods 
are compatible with each other. In order to get the full story the interested reader should have 



a look at [123, 137, |13§| [. In particular the first paper contains very useful material. 

In this section I will present the results of my simple scanning analysis and of a Gaussian 
analysis by Stefan Schael who used the same input parameters. For the rest of these lectures I 
will only use simple scanning, except for e 1 je where also the results of the Gaussian method will 
be presented. The quoted results of the scanning method for a quantity Q given below should 
be understood as follows: 

{Q = A±B} = {A-B<Q<A + B} (4.89) 

This means that the central value A does not generally correspond to central values of the input 
parameters. 



4.2.3 Output of the Standard Analysis 

Using simultaneously the five steps discussed above one finds the allowed region of (g,fj)- In 
fig. H we show the result of an analysis by Stefan Schael which uses the input parameters of 
table |2[ Only the dark region is allowed. From this figure one extracts 

a = 93° ±11°, (3 = 23.6 tts°, 7 = 64° ± 11°, (4.90) 

sin 2(3 = 0.73±°;°I, \V td \ = (8.0 ± 0.7) • lO" 3 . (4.91) 
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In this analysis Gaussian errors for all input parameters have been used. My own, more conser- 
vative analysis that uses scanning for all input parameters gives 

78.8° < a < 120°, 15.1° </?< 28.6°, 37.9° < 7 < 76.5° (4.92) 

sin 2/3 = 0.67 ±0.17, \V td \ = (8.0 ± 1.3) • 10~ 3 . (4.93) 

The "true" errors are probably between these two estimates. The results of both analyses are 
summarized in table ||. 




P 

Figure 9: Conservative Unitarity Triangle as of January 2001. 

The allowed region for (g,fj) resulting from the scanning method is presented in fig. ||[ It 
is the shaded area on the right hand side of the circle representing the lower bound for (AM) S , 
that is (AM) 8 > 15/ps. The hyperbolas in fig. || give the constraint from e and the two circles 
centered at (0,0) the constraint from |Kt&/^cft|- The circle on the right comes from B® — B® 
mixing and excludes the region to its right. We observe that B^ — B^ mixing is almost ineffective 
for the chosen ranges of the input parameters within the SM and the allowed region is governed 
by |Ktb/Vrf,|, AM S and e. We also observe that the region g < is practically excluded by the 
lower bound on AM S . It is clear from this figure that (AM) S is a very important ingredient 
in this analysis and that the measurement of (AM) S giving also lower bound on Rt will have 
a large impact on the plots in figs. || and [9|. Finally we find that whereas the angle [5 is rather 
constrained, the uncertainties in a and 7 are substantially larger. 
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Table 3: Output of the Standard Analysis. X t = V t *Vtd- 



Quantity 


Scanning 


Gaussian 


1 t/ . 1 /i n — 3 


P. 7 Q 3 


O.U lu. 1 


I Vts/Vcb | 


0.979 - 0.993 


0.984 ± 0.004 


1 Vtd/Vts I 


0.16 - 0.22 


0.20 ±0.02 


sin(2/3) 


0.50 - 0.84 


7S+ - 07 


sin(2a) 


-0.87 - 0.38 


-0.1 ±0.3 


sin(7) 


0.61 - 0.97 


0.9 ±0.1 


ImAt/10" 4 


0.94-1.60 


1.2 ±0.2 


fj 


0.22 - 0.46 


0.36 ± 0.07 


Q 


0.06 - 0.34 


0.18 ±0.09 



4.2.4 An Upper Bound on AM S 

In view of the expected measurement of AM S in 2001 at Tevatron it is of interest to find its 
upper bound within the SM. The most straightforward manner to obtain this bound is to set 



\Vts\ max — 0.043, TR\,{rn%) m ax — 

171 GeV and {\J B Bs F Bs ) max = 300 MeV fl2|, |12| in (p2j) . 

We find 

AM S < 27.5/ps (4.94) 
with 16.0/ps obtained for central values of the input parameters. 

4.3 First Conclusions 

In this section we have completed the determination of the CKM matrix. It is given by the values 
of \ V U s\, \Vcb\ an d \ V u b\ in ( |1.38j ) and ( |1.39| ), the results in table ^ and the unitarity triangle shown 
in figs. H and ^. We should stress, that soon this analysis will be improved through the new 
information on the angle (3 coming from B factories and Tevatron and on AM S from Tevatron. 
In particular the latter measurement should have an important impact on the allowed area in 
the (g, fj) plane. 

We conclude that the SM is capable of describing the observed indirect CP violation in Kl 
decays, taking into account the data on B^ s — B^ s mixings, \V c b\ and |lub/K*|- We also observe 



that Rb and Rt from {V^/V^l and B^ s — B^ s mixings alone satisfy the condition ( 1.36 ). Taking 
{Rb)min = 0.30, this condition reads 0.70 < Rt < 1.30. From fl4.88|) and the lower bound on 
AM S one has Rt < 1.01. On the other hand |Vtd| m j n is governed by B® — B^ mixing. From 
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and ( p3|) one has then R t > 0.71. Consequently CP violation in B-decays is predicted 
on the basis of | / V^t, | and B® s — B^ s mixings alone, even if our conservative analysis does 
not show this so clearly. Indeed the first result for the CP asymmetry in — > ip^s with 
sin 2/3 = 0.79 ±0.42 from CDF [p0|] , presented in 1999, established CP violation in B decays at 
93% C.L. and moreover was consistent with ( 4.91 ) and fl4.93j ). Yet the reality could turn out to 
be different as we will see in a moment. 



4.4 First Results on sin 2/3 from B Factories 

Since the summer 2000 we have the first results on the time dependent CP asymmetry, a^x 



from BaBar [ 141 1 and Belle [142]. These results indicate that the value of the angle [3 in the 
unitarity triangle could turn out to be substantially smaller than expected on the basis of the 
standard analysis of the unitarity triangle within the SM presented above. Indeed the three 
measurements of this asymmetry 

{0.79 ± 0.42 (CDF) [140] 

0.12 ±0.37 ±0.09 (BaBar) ||| (4.95) 
0.45 ±0.44 ±0.08 (Belle) [|l42| 

imply the grand average 

(sin 2(3)^ Ks = 0.42 ± 0.24 . (4.96) 

This should be compared with the results of the standard analysis of the unitarity triangle within 
the SM presented above which gave sin 2/3 = 0.73±° : ?I and sin 2/3 = 0.67 ± 0.17 for the gaussian 
and scanning analyses respectively. Similar values can be found in the literature 



(sin 2/3) 



SM 



0.69 ±0.07 
0.70 ±0.24 
0.63 ±0.12 



[123] 



(4.97) 



where the last two results represent 95% C.L. ranges. Clearly in view of the large spread of 
experimental results and large statistical errors in ( |4,95| ), the SM estimates in ( [4.91| ), ( 4.93| ) 
and ( 4.97|) are compatible with the experimental value of (sm2/3)^K s m ( 4.96] ). Yet the small 
values of sin 2/3 found by BaBar and Belle may give some hints for new physics contributions to 



B® — B^ and K° — K° mixings. In particular as discussed recently in several papers [143]-[14C] 



new CP violating phases in B® — B® and K° — K° mixing could be responsible for small values 
of sin 2/3 in ( |4.96 ). In this context an absolut lower bound on sin 2/3 in models with minimal 



flavour violation has been derived in [147]. We will discuss this bound and its update in detail 
in Section 9. Finally a general model independent analysis of new physics in B — » ipK^ and 



Bj — > ipKs has been presented very recently in |14 
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On the other hand as stressed in [143] the SM estimates of sin 2/3 are sensitive to the assumed 
ranges for the parameters 

\Vcbl \V ub /V cb \, B K , \fs~ d F Bd , £, (4.98) 

that enter the standard analysis of the unitarity triangle. While for "reasonable ranges" (see 
table |2|) of these parameters, values of sin 2/3 < 0.5 are essentially excluded, such low values 



within the SM could still be possible if some of the parameters in ( 4.98j ) were chosen outside 
these ranges. In particular for [Ktfe/^cbl < 0.06 or Bk > 1-3 or £ > 1.4 the value for sin 2/3 lower 
than 0.5 could be obtained within the SM. We agree with these findings. 

In what follows we will assume the " reasonable ranges" for the parameters as given in table || 
and we will use the values of the CKM parameters determined in this section to predict various 
branching ratios for rare and CP-violating decays. This means we will not take into account 
in our numerical analyses the recent results from BaBar and Belle. We will return to them in 
Section 9 where some aspects of the physics beyond the SM will be discussed. Thus the next four 
sections apply only to the SM and could be affected considerably if the improved measurements 
of the asymmetry a^,K s by BaBar, Belle, CDF and DO will confirm with higher precision the 
present low values of sin 2(3 from BaBar and Belle. This would be truely exciting! 

5 e'/e in the Standard Model 
5.1 Preliminaries 

Direct CP violation remains one of the important targets of contemporary particle physics |149f| . 
In the case of K — > tttt, a non- vanishing value of the ratio Re(e'/e) would give the first signal for 



direct CP violation ruling out superweak models [15C] in which e'/e = 0. Until February 1999 



the experimental situation on e'/e was rather unclear: 



Re(e'/e) = { (NA31) [151] , 

1 (7.4 ±5.9) • 10- 4 (E731) pi. 



While the result of the NA31 collaboration at CERN [151] clearly indicated direct CP violation, 



the value of E731 at Fermilab [152], was compatible with superweak theories. This controversy 
is now settled in favour of NA31. The most recent experimental results for the ratio e'/e from 
Fermilab and CERN presented during 1999 and 2000 read 

H*7.) = < f 28 - *")- 1 "- 4 (KTeV) [153], 
} [14.0 ± 4.3) • 10~ 4 (NA48) fll5j . 
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Together with the NA31 measurement these data confidently establish direct CP violation in 



nature and taking also the E731 result into account one finds the grand average [154 

Re(e'/e) = (19.2 ± 2.4) • 10~ 4 , (5.3) 

close to the NA31 result. 

While the experimentalists should be congratulated for these results, the situation is clearly 
unsatisfactory. The KTeV result is by a factor of four larger than the E731 result. In addition 
the substantial difference between KTeV and NA48 results is disturbing. Let us hope that these 
issues will be clarified soon. In this context an independent measurement of e'/e by KLOE 
at Frascati, which uses a different experimental technique than KTeV and NA48 will be very 
important. 

There is a long history of calculations of e'/e in the SM beginning with the pioneering 



calculations in 155 ] and |156f| and the first extensive phenomenological analyses in [157]. Over 
the 1980's these calculations were refined through the inclusion of QED penguin effects for 



m t <C M\v JL58| , 159, 160 1 , the inclusion of isospin breaking in the quark masses J159| , |160| , |161| j, 
through improved estimates of hadronic matrix elements in the framework of the 1 /N approach 
[162] and the inclusion of QCD penguins, electroweak penguins (7 and Z° penguins) and the 
relevant box diagrams for arbitrary top quark mass [162, 164]. The strong cancellation between 
QCD penguins and electroweak penguins for mt > 150 GeV found in these two papers was 



confirmed by other authors [165j. 

During the 1990's considerable progress has been made by calculating complete NLO correc- 



tions to e E0]-[H3]. Together with the NLO corrections to e and B — B mixing pfi, 52 



this allowed a complete NLO analysis of e'/e including constraints from the observed indirect 
CP violation (e) and B° ds - B° ds mixings (AM^ jS ). The improved determination of the V u b and 
V c b elements of the CKM matrix and in particular the determination of the top quark mass mt 
had of course also an important impact on e'/e. 



While until 1999 only a few groups were involved in detailed analyses of e'/e |42|,[163]- |169[| , 
the new results from KTeV and NA48 prompted a number of theorists to join the efforts to 
calculate this important ratio. By the end of 2000 the theoretical e'/e-club has more than 30 
members. 

The developments of the last two years include new estimates of hadronic matrix elements 
using lattice simulations |122j ], various versions of the large-N approach, chiral quark model, 
QCD sum rules and other approches mentioned below. Simultaneously extensive studies of final 
state interactions (FSI), isospin breaking effects and electromagnetic effects in e'/e have been 
pursued. Finally the impact of new physics on e'/e has been analyzed in various extentions of 
the SM. 
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We will briefly review these topics below. To this end we have to recall a number of useful 
formulae for e'/e. Other reviews of e'/e can be found in [169]- [172] 



5.2 Basic Formulae 



The parameter e' is given in terms of the isospin amplitudes Aj in ( 3.26 ). Applying OPE to 



these amplitudes one finds [42, 167 1 



— = IniA t • F e >, 

£ 



where 



with 



F £ > = 

p(l/2) 
p(3/2) 



p(l/2) _ p(3/2) 



exp(i$), 



E%(Oi>o(l-filB) 

-Y\vi{Qi}2 • 

1 , 1 ' ■ 



r 
r 



Here 



Since 



G-puj 



2 Id Re A, 



(Qi)/ = ((vr^/IQilif) 



Co' 



Re^2 
Re^o 



$ £ ~0, 



(5.4) 
(5.5) 

(5.6) 
(5.7) 

(5.8) 
(5.9) 



F £ i and e'/e are real to an excellent approximation. The operators Qi have been given already 
in (|2.4|)-( |2~8| ), The Wilson coefficient functions yi{p) were calculated including the complete 
next-to-leading order (NLO) corrections in ]4Q| — p4]| . The details of these calculations can be 



found there and in the review ||. Their numerical values for Ai^ corresponding to ajijL(Mz 



™ corresponding to cJ^o 

Mo M o 

0.119 ± 0.003 and two renormalization schemes (NDR and HV) are given in table || |173f |. 

It is customary in phenomenological applications to take ReAo and u from experiment, i.e. 

Re^ 2 



ReA = 3.33 • 10"' GeV, 



Re^n 



0.045, 



(5.10) 



where the last relation reflects the so-called AI = 1/2 rule. This strategy avoids to a large 
extent the hadronic uncertainties in the real parts of the isospin amplitudes Aj. 

The sum in (5.6) and ( |5.7| ) runs over all contributing operators. p( 3 / 2 ) is fully dominated 
by electroweak penguin contributions, i-^ 1 / 2 ) on the other hand is governed by QCD penguin 
contributions which are suppressed by isospin breaking in the quark masses (m u ^ nid)- The 
latter effect is described by 

1 (Irm4 2 )iB 



Im^n 



(5.11) 
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Table 4: AS = 1 Wilson coefficients at fj, = m c = 1.3 GeV for nit = 165 GeV and / = 3 effective 
flavours. y\ = 3/2 = 0- 





MS 


290 MeV 


MS 


340 MeV 


MS 


390 MeV 


Scheme 


NDR 


HV 


NDR 


HV 


NDR 


HV 


2/3 


0.027 


0.030 


0.030 


0.034 


0.033 


0.038 


2/4 


-0.054 


-0.056 


-0.059 


-0.061 


-0.064 


-0.067 


1/5 


0.006 


0.015 


0.005 


0.016 


0.003 


0.017 


2/6 


-0.082 


-0.074 


-0.092 


-0.083 


-0.105 


-0.093 


2/7/a 


-0.038 


-0.037 


-0.037 


-0.036 


-0.037 


-0.034 


W« 


0.118 


0.127 


0.134 


0.143 


0.152 


0.161 


2/9 A* 


-1.410 


-1.410 


-1.437 


-1.437 


-1.466 


-1.466 


yio/a 


0.496 


0.502 


0.539 


0.546 


0.585 


0.593 



The main source of uncertainty in the calculation of e'/e are the hadronic matrix elements 
{Qi)i and r^iB as we will see later on. {Qi)i generally depend on the renormalization scale \x 
and on the scheme used to renormalize the operators Q{. These two dependences are canceled 
by those present in the Wilson coefficients yi(n) so that the resulting physical e'/e does not (in 
principle) depend on [i and on the renormalization scheme for the operators. Unfortunately, the 
accuracy of the present non-perturbative methods used to evalutate (Qi)i is not sufficient to have 
the \i and scheme dependences of (Qi)i fully under control. We believe that this situation will 
change once the lattice calculations improve. A brief review of the existing methods including 
most recent developments will be given below. 



5.3 An Analytic Formula for e'/e 

The basic formulae for e'/e presented above are not very transparent and we would like to 
improve on this. To this end let us first determine as many matrix elements (Qi)i as possible 
from the leading CP conserving K — > -kit decays, for which the experimental data is summarized 
in ( |5.10 ), This determination is particularly easy if one sets fj, = m c j42|. The details of this 
approach will not be discussed here. It sufficies to say that this method allows to determine 
only the matrix elements of the (V — A) (g) (V — A) operators. Explicit formulae can be found in 
| J42j | . For the central value of ImXt these operators give a negative contribution to e'/e of about 
—2.5 • 10~ 4 . This shows that they are only relevant if e'/e is below 1 • 10~ 3 . 

Unfortunately the method in [j42j | does not provide the matrix elements of the dominant 
( V — A) (g) (V + ^4) operators and one has to use non-perturbative methods to estimate them. In 
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this context the calculations of Bertolini and collaborators [|169|1 within the chiral quark model 
should be mentioned. These authors have determined the parameters of this model by taking 
the constraints ( |5.10| ) into account. In this manner they were able to determine not only the 
matrix elements of the (V — A) ® (V — A) operators, as in [42], but also the matrix elements of 



the (V — A) <8> (V + A) operators. This is clearly interesting. On the other hand, it is not clear 
how well the chiral quark model approximates QCD and consequently whether the extracted 
values of the relevant matrix elements have anything to do with the true QCD values. For this 
reason we will not use them here. 

In order to exhibit the matrix elements of the dominant (V — A) <8> (V + A) operators in an 
analytic formula for e'/e, we first express them in terms of non-perturbative parameters B^ 2 ^ 
and B^ 2 ^ as follows: 

(Qi)o ee B? IV (Q^ , (Qih ee B? /2) (Q t )t c) • (5-12) 

The label "vac" stands for the vacuum insertion estimate of the hadronic matrix elements in 

question for which B^ 2 ^ = B^ 2 ^ = 1. 

The Wilson coefficients y$, yj and ys are much smaller than ?/q. On the other hand the 

contributions of (Qr,8)2 are enhanced by the factor 1/lo ps 22 as seen in (f>.7\). But yj is 

substantially smaller than y$ and ((^7)0,2 are colour suppressed. It is then not surprising [Q 

that e'/e is only weakly sensitive to the values of the parameters B^ 2 \ B^ 2 \ B^ 2 ^ and 
(3/2) 

B-f as long as their absolute values are not substantially larger than 1. On the basis of 
existing non-perturbative approaches we can assume that this is indeed the case. 
Following p2| we set then 



d(1/2)/ \ -1 d(1/2)/ \ d(1/2)/ x _ D (l/2) d(3/2)/ x d(3/2)/ n _ D (3/2) 

B 7,8 ( m c) = 1, £5 (m c ) = Bq >(m c ) = By By >{m c ) = B% (m c ) = By ' . 

(5.13) 

This strategy |12| allows then to express e'/e or equivalently F e > in terms of of the following 
parameters 

The appearence of m s originates in 

(Q 6 ) { y c) ~ l/m 2 M, (Q 8 )t C) ~ 1MV) • (5-15) 
A^. and m-t enter through the coefficients yi. 



With all this information at hand, it is possible |167| , 174] to cast the formal expressions for 
e'/e in (|5.4|)-( |5"7r] ) into an analytic formula which exhibits the dependence on the parameters in 
( 5.14| ). We set Qib = 0.16 for the moment. We will return to it below. The analytic formula 
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given below, while being rather accurate, exhibits various features which are not transparent in 
a pure numerical analysis. It can be used in phenomenological applications if one is satisfied 
with a few percent accuracy. Needless to say, in the numerical analysis |173| , |109| 1 presented 
below we have used exact expressions. 

In this formulation the function F £ > is given simply as follows: 



F £ , =P + P X X (x t ) + P Y Y (x t ) + P z Zo(x t ) + P E E (x t ) 



(5.16) 



with the mt-dependent functions given in subsection 2.6. 

(1/2) 



The coefficients Pi are given in terms of B G 

Pi 

where 



Ug 3//2 ' ) and m s (m c 



as follows: 



Rk = B, 



(1/2) 



137 MeV 



i 2 



m s (m c ) + m d (m c ) 



i?« = B. 



(3/2) 



137 MeV 



m s (m c ) + m d (m c ) 



(5.17) 



(5.18) 



The Pi are renormalization scale and scheme independent. They depend, however, on A^. In 



table | we give the numerical values of r 



(0) (6) 



and r • for different values of at fi = m c in 



(i) 



MS 



the NDR renormalization scheme [10£]. Actually at NLO only ro coefficients are renormalization 
scheme dependent. The last row gives them in the HV scheme. The inspection of table || shows 
that the terms involving and dominate the ratio e'/s. Moreover, the function Zo(xt) 
representing a gauge invariant combination of Z°- and 7-penguins grows rapidly with m t and 
due to < these contributions suppress e'/e strongly for large m t |163| , |164f| • 



Table 5: Coefficients in the formula ( |5.17l ) for various A^ in the NDR scheme. The last row 
gives the ro coefficients in the HV scheme. 





AgL = 290MeV 




1 = 340 MeV 


M 


1= 390 MeV 




-(0) 


-(«) 


r (8) 


JO) 


r (6) 


r (8) 


JO) 


J6) 


J8) 


i 


i 


i 


i 


i 


i 


i 


i 


i 


i 





-3.122 


10.905 


1.423 


-3.167 


12.409 


1.262 


-3.210 


14.152 


1.076 


X 


0.556 


0.019 





0.540 


0.023 





0.526 


0.027 





Y 


0.404 


0.080 





0.387 


0.088 





0.371 


0.097 





Z 


0.412 


-0.015 


-9.363 


0.474 


-0.017 


-10.186 


0.542 


-0.019 


-11.115 


E 


0.204 


-1.276 


0.409 


0.188 


-1.399 


0.459 


0.172 


-1.533 


0.515 





-3.097 


9.586 


1.045 


-3.141 


10.748 


0.867 


-3.183 


12.058 


0.666 
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Finally by investigating numerically the formula ( |5.16| ) it is possible to find a crude approx- 
imation for F £ r. 



F £ , « 13 • 



110 MeV 



1 2 



m s (2 GeV) 



6 V y 8 V 165 GeV 



MS 



340 MeV 



. (5.19) 



This formula while exhibiting very transparently the dependence of e' je on the parameters ( [5.14] ) 
should not be used for any serious numerical analysis. 
We observe the known features: 

• F £ i increases with increasing A^L. In fact the renormalization group effects play an im- 
portant role here. If one did not include them e' je would be typically 0(1O~ 5 ). On the 
other hand the formula ( |5.19| ) is only valid for A^L > 200 MeV. F £ / does not vanish for 

• F e i increases with increasing B^ 2 ^ that represents QCD penguins. It decreases with 

(3/2) 

increasing B$ that represents electroweak penguins. The latter contribution increases 



with increasing m t . 

(1/2) (3/2) 

• The partial cancellation between QCD penguin (Bq ) and electroweak penguin (B% ) 

(1/2) (3/2) 

contributions requires accurate values of Bg and B% for an acceptable estimate of 
e'/e. Because of the accurate value mt(mt) = 166 ± 5 GeV, the uncertainty in e'/e due 
to the top quark mass amounts only to a few percent. 

• The 1/m, 2 dependence is an artifact of the decomposition of the matrix elements into B{- 
factors and the vacuum insertion matrix elements as given in ( |5.12| ), We will return to 
this point below. 

5.4 The Status of m s , B^ /2 \ Bf ,2 \ VL m and 
5.4.1 m s 

The most recent values for m s (2 GeV) extracted from lattice calculations, QCD sum rules and r— 
decays are summarized in table H The lattice result is the most recent "world average" of Lubicz 



|175| ], See also | |180| ]. It is higher than the older result from Gupta [181] but could decrease a 
bit when the calculations with dynamical fermions are completed. Significant progress has 
been made in the last two years by the introduction of improved actions and non-perturbative 
renormalization techniques. Quenched calculations give m s (2GeV) = (110 ± 15) MeV and the 
increased error in table ^ gives the estimate of the quenching error based on the caculations with 
two flavours. Details with references to calculations of different groups can be found in |17 
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The recent QCD sum rule values are considerably smaller than the older ones [182]. The same 
comment applies to m s from r-decays. 

While the situation regarding m s improved considerably during the last two years and the 
values obtained using different methods are coming closer to each other, the error on m s is still 
large. In our numerical analysis of e'/e we will use 

, , f (110 ±20) MeV / u = 2GeV 
mJu) = < v ' 5.20 

[ (130 ± 25) MeV \i = m c 

which is in the ball park of the results given in table ||. Here m c = 1.3 GeV. 

Five years ago values like m s (2GeV) = 130 — 150 MeV could be found in the literature. 
Beginning with the work of Gupta [ |181| ] the values of m s decreased considerably, which makes 
the values of e'/e larger. On the other hand QCD sum rules allow to derive lower bounds on the 
strange quark mass. It is found that generally m s (2GeV) > 100 MeV fl83| . We observe that 
lattice results are very close to this bound. 



Table 6: m s (2GeV) obtained using various methods. 



Method 


m s (2GeV) [MeV] 


Reference 


Lattice 


110 ±25 


Lubicz 


175 


1 


QCDS 


119 ±12 


Narison 


176| 


QCDS 


116 ±22 


Jamin 


177 




QCDS 


115 ±8 


Maltman 17 


8| 


r-Decays 


114 ± 23 


Pich-Prades ] 


79 1 



The large sensitivity of e'/e to m s has been known since the analyses in the 1980's. In the 



context of the KTeV result this issue has been analyzed in [184]. It has been found that provided 

(1/2) (3/2) 

2B K & ' ' - B£< ' < 2 the consistency of the SM with the KTeV result requires the 2a bound 
m s (2GeV) < 110 MeV. Our analysis below is compatible with these findings. 

On the other hand, it has been emphasized by Guido Martinelli, Eduardo de Rafael and 
other researchers in the field that the study of e'/e as a function of m s , independently of B^ 2 ^ 

(3/2) 

and .Bg , could be misleading. Indeed a correlation could exist between m s and the latter 
parameters. In the large-N approach of |162, 185] such correlation is not observed. However, 
I can imagine that in an approach, like lattice simulations, that calculates the hadronic matrix 
elements directly, without parametrizing them in terms of m s and i?j-factors, some correlations 
in values of these parameters could certainly be present. It will be interesting to study this issue 
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once the lattice calculations improve. 

Prom my point of view, the strange quark mass must be somehow related to the matrix 
elements (Qe)o and (Qs)2- Afterall Qe and Qs are density-density operators and in the large-N 
limit the anomalous dimensions of these operators are equal to — 27 m , where 7 m is the anomalous 
dimension of the mass operator |160fl . Consequently the //-dependence of (Qe)o and (Qs)2 is 
governed to a very good approximation by the //-dependence of m s (//) as shown in ( |5.15| ). This 
has been verified numerically beyond the large-N limit in |^ ]. This implies that B^ 2 ^ and 
-Bg 3//2 ^ are practically //-independent and from the point of view of the //-dependence m s is 

(1/2) (3/2) 

indeed uncorrelated with B§ and B$ . On the other hand I can imagine that at fixed /t, 

(1/2) (3/2) 

there is some correlation between the actual values of m s , Bg and B% . 



In this context it should be recalled [160] that the //-dependence of (Qq)q and (Qs)2, just 



discussed, is to a very good approximation canceled by the /i-dependence of ye(fJ-) and ys(^) 
respectively. This means that the unphysical /i-dependence in e'/e is very small already in 
the present calculations. This is a nice feature, which can be seen explicitly within the lare-N 
approach. 

5.4.2 2#/2) and B (3/ 2) 

The values for B^ 2 ^ and B^ 2 ^ obtained in various approaches are collected in table [?]. The 
lattice results have been obtained at // = 2 GeV. The results in the large-N approach and 
the chiral quark model correspond to scales below 1 GeV. However, as explained above, B^ 2 ^ 

(3/2) 

and Bg are only weakly dependent on //. Consequently the comparison of these parameters 
obtained in different approaches at different // is meaningful as long as // > 0.8 GeV. 

Next, the values coming from lattice and chiral quark model are given in the NDR renormal- 
ization scheme. The corresponding values in the HV scheme can be found using approximate 
relations 1173|1 

(4 1/2) )hv « 1.2(4 1/2) ) N dr, (^ 3/2) )hv « 1.2(i4 3/2) ) N DR. (5.21) 
The scheme dependence of these parameters in the large-N approach are to my knowledge not 



yet under full control but some progress is being made 128 , 190 , 192 1 



Concerning the lattice results for B^ 2 \ the old results read B^q 2 \2 GeV) = 1.0 ± 0.2 



1 193 , 194]. However, over the last years it has been realized beginning with [|195[] that lattice 



calculations of B^ 2 ^ are very uncertain. One has to conclude that there are no solid predictions 
for B^ 2 ^ from the lattice at present. There is a hope that soon results using domain wall 

(1/2) (3/2) 

fermions not only for B§ 1 ' and Bg ' ' but for all hadronic matrix elements entering e /e will be 
available. A very informative and critical review of the present situation has been given recently 



by Laurent Lellouch [ 122 ], 
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The result for B$ from the dispersive approach [ |189[ | corresponds to m s {2 GeV) ~ 100 MeV. 



In this approach B^ 2 ^ increases with m s . The results for B^ 1 *' from the 1/N approach of de 
Rafael and collaborators are not yet available. Their analysis of the electroweak operator Qj 



(3/2) 



can be found in [ 19C ] 



(1/2) (3/2) 

Table 7: Results for Bf' J and B%'*' obtained in various approaches. 



Method 


R (l/2) 


R (3/2) 


Lattice[ 186, 187, 188] 
Large-N|185, 126] 
ChQMjl69| 
Dispersive Approach |189j] 
Large-N (NJL)]l90|] 
Sum Rules PH 


0.72 - 1.10 
1.07-1.58 

2.9 ±0.7 
2.8 ±0.5 


0.69 - 1.06 
0.42 - 0.64 
0.75 - 0.79 
1.11 ±0.16 ±0.23 

1.3 ±0.2 

1.7±0.4 



Table demonstrates very clearly that there is no consensus on the values of B^ 2 ^ and 



B 



(3/2) 



In particular the last two groups find both parameters to be substantially higher than 



obtained by the other groups. It should be remarked that the Dortmund group |185, 126] 



generalized the previous leading order large-N calculations |124j , |160 , 162] to include higher 
order 1/N corrections. They find that in the chiral limit 0(p 2 /N) correction not included in 
the result in table [7| enhance B^ 2 ^ to roughly 1.5. As the chiral logs have not been taken into 
account it is probably premature to use this result in phenomenological applications. 

Biased to some extent by the results from the large-N approach and lattice calculations, we 

(1/2) (3/2) 

will use in our numerical analysis below B§ and Be in the ranges: 



B, 



(1/2) 



1.0 ±0.3, 



B. 



(3/2) 



0.8 ±0.2 



(5.22) 



keeping always B, 



(1/2) 



> B 



(3/2) 



5.4.3 $7tr and Ai£ 



MS 



The older estimates of ^ib in the 1/N approach [16C] and in chiral perturbation theory [15£, 161] 



gave the value 0.25 ± 0.10. The most recent refined calculation in [197] gives 

n m = 0.16 ±0.03 . 



(5.23) 



We will adopt this value here. On the other hand the story of ^ib is not finished yet. We will 
discuss it below. 
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The dependence of e'/e on S7ib can be studied numerically by using the formula (5.6) or 



incorporated approximately into the analytic formula ( |5.16 ) by simply replacing B^ 2 ^ with an 



effective parameter 



r R (i/2h 



B, 



(i /2) (1 - 0.9 gig) 
0.856 



(5.24) 



A numerical analysis shows that using (1 — ^ib) overestimates the role of JIib- 

In table ^ we summarize the input parameters used in the numerical analysis of e'/e below. 
The range for A^L in table |8| corresponds roughly to a s (Mz) = 0.119±0.003, which is very close 
to the recent determination in pi]] . 

(1/2) (3/2) 

Table 8: Collection of input parameters. We impose 



Quantity 



Central 



Error 



Reference 



MS 

m s (m c ) 



B. 



(1/2) 
6 

(3/2) 



ImA t 
ImA t 



340 MeV 
130 MeV 
1.0 

0.8 

1.33 

1.33 



±50 MeV 
±25 MeV 

±0.3 

±0.2 
±0.14 (G) 
±0.30 (G) 



See Text 
See Text 
See Text 



[173] 



mi 



5.5 Numerical Results for e'/e 



We list the results from various groups in table |9[ The labels (G) and (S) in the second column 
stand for two error estimates: "Gaussian" and "Scanning" respectively. The result from the 



Munich group given here is an update of the analysis in |173| 1 done in 109 1. It uses the input 
parameters of table [8] where the value of ImAf from [10S, 173] is slightly higher than the one 
found in the previous section. Similarly the result of the Rome group [199] is the most recent 



estimate given in [171]. 



In [173, 109, 199, 171] e'/e has been found to be typically by a factor of 2-3 below the data 
and the KTeV result in (|5.2| ) could only be accomodated if all relevant parameters were chosen 
simultaneously close to their extreme values. On the other hand the NA48 result is essentially 



compatible with [173, 109, |199|, 171] within experimental and theoretical errors. Higher values 



of e'/e than in |73|, Ell HHi IIZ1> in the ballpark of (|J), have been found in g69], [L9|, [191 
190, |20j, |20^]. The result in (200| corresponds to B^ /2) 



B, 



(3/2) 



1 as the Dubna-DESY 



group has no estimate of these non-perturbative parameters. Recent reviews can be found in 
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|16Sl- ||172]1 . Furthermore it has also been suggested that the final state interactions (FSI) could 
enhance e' je by a factor of two [169, 203 , 204]. We say a few words about it below. 



Reference 


p i Ip rin -4 ! 

c / c [1U 


Munich 
Munich 


172 
172 


,s 


Q n+6.8 (fl\ 


Rome 
Rome 


L99 
L99 


171 
171 




°- J --9.S IW 

_13 n _» 37 n (S) 


Trieste 1 169 
Trieste 169 


1 
1 


22 ± 8 (G) 
9 -> 48 (S) 


Dortmund 


[198] 


6.8 -> 63.9 (S) 


Montpellier [191 




24.2 ±8.0 


Granada-Lund 


19CJ 


34 ± 18 


Dubna-DESY | 


20C| 


-3.2 -» 3.3 (S) 


Taipei [2 


01 




7 -> 16 


Barcelona- Valencia [|203| 


17 ±6 


Beijing [2 


02] 


16.2 39.3 



Table 9: Results for e'/e in the SM in units of 10 -4 . 



5.6 Renormalization Scheme Dependence 

All the results presented here apply to the NDR scheme. They are lower by roughly 30% in the 

(1/2) (3/2) 

HV scheme if the same values for (Bq ' , Bg ' ) are used. On the other hand, if simultaneously 
(B^ 2 \ B^ 2 ^) are transfered to the HV scheme by means of (5.21), the scheme dependence is 
reduced. However, the game of changing renormalization schemes is only meaningful if a given 
non-perturbative method has renormalization scheme dependence of the matrix elements fully 
under control. This is not really the case in the present approaches. Future lattice calculations 
have the best chance to make progress here, but it is very desirable to have also analytic solutions 
to this problem. 

In this context it should be emphasized that the cancellation of the renormalization scheme 
dependence in the electroweak penguin sector requires to go beyond the NLO approximation 



for yi(/J,) [173|. The reason is that the dominant effect of the electroweak penguins related to 
the functions Xq, Yq and Zq in ( |5.16j ) enters first at the NLO level in the renormalization group 
improved perturbation theory. Consequently the issue of the scheme dependence in this sector is 
shifted to the NNLO level. In particular QCD corrections to the one-loop functions Xq, Yq and 
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Zq have to be calculated. This calculation is now available |80|| . However, in order to complete 
the NNLO analysis, corresponding corrections to gluon penguin diagrams and the three loop 
anomalous dimension matrices must be calculated. A very difficult task. The general structure 
for Wilson coefficients at NNLO can be found in 



5.7 Final State Interactions 

What is the role of final state interactions in e' /el This question has been addressed recently 



by several authors [169, 198, |203| , 204]. Actually the issue of FSI in K —* 7T7t decays is not 



new. It has been known from the analyses in |162|, 205] that the pion loops representing FSI 



provide a sizable enhancement of A J = 1/2 transitions and suppression of AI = 3/2 transitions 
helping in the explanation of the AI = 1/2 rule. The analyses addressing this rule are given in 



p06| , |202j , 207], where further references can be found. 

If this pattern of enhancements and suppressions also applies to the matrix elements of the 
penguin operators relevant for e'/e, this would mean that FSI effects enhance the QCD penguin 
contributions and suppress the electroweak penguin contributions making e'/e larger than in 
the case of neglecting these effects. In this context an interesting proposal has been made by 
Pallante and Pich [ |203| ], who following an older work by Truong [ |208| ], suggested that the FSI 
effects in question can be unambiguosly resummed to all orders in chiral perturbation thory 
using the so-called Omnes factor. Assuming in addition that this factor should be the dominant 
FSI effect relative to large-N estimates of hadronic matrix elements, they found that the result 



for e'/e from the Munich group [173, 109] can be enhanced by roughly a factor of two. 

A critical analysis of this suggestion has been presented in [ |209| ] . A nice summary of the 
points made in this paper can be found in |210|]. Personally, I think that the sign of the effect 



calculated by Pallante and Pich could be correct but in view of the critics made in |209| , the 
actual size of FSI in the evaluation of e'/e given in [ |203| ] cannot be really trusted. Most probably 
a satisfactory solution to the problem of FSI can only be achieved in a complete calculation of 
the hadronic matrix elements within a self-consistent framework. In this context, an interesting 
recent proposal in ] ]211| |, if realized, could put the inclusion of FSI in the lattice calculations of 
hadronic matrix elements in principle under control. 

5.8 Isospin Breaking Effects 

Isospin breaking effects are generated by the m u — difference and by electromagnetic correc- 
tions. While they are generally small in the kaon system (roughly 1%) and can be neglected in 



view of other uncertainties, they may have an impact on e'/e. Indeed as seen in ( 5.11 ), fiie is 



enhanced by 1/u ~ 22, implying that f2iB in the ball park of 0.20 could certainly be expected. 
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Let us first concentrate on the isospin breaking corrections due to m u ^ m^. In the lowest 
order of chiral perturbation theory ((m u — md)p ) only the 7r° — rj mixing in the strong La- 
grangian contributes and one finds Ojb = 0.13 |161| , |197f |. The older estimates of ((m u — m^)p 2 ) 
contributions to S7ib , which include also 7r° — rj — rj' mixing, within the large-N approach |160|] 



and in chiral perturbation theory |15S| , 161 ] gave the total value 0.25 ± 0.10. The most recent 
refined calculation in [ |197| ] gives somewhat smaller value J7ib = 0.16 ±0.03 that we have adopted 
in our numerical analysis. 



However, as stressed in [212, 213], at this level also 0{p ) weak counterterms have to be 



taken into account. At present these counterterms can only be estimated by making some 



specific dynamical assumptions. As demonstrated in |212| , |213| 1 these additional contributions 
can compete with the ones included sofar. Interestingly, they have the tendency of decreasing 
fiiB and even to reverse its sign making s' je larger. Unfortunately a quantitative estimate of 
these contributions is plaqued by very large uncertainties. Consequently, if one wants to be 
conservative, an error of 100% should be assigned to f^iB with a central value arround 0.1. 

Concerning the electromagnetic corrections, a recent analysis indicates that they have only 
a small impact on Ojb [214, |21S[|. On the other hand, as pointed out in [215, 216], the electro- 



magnatic corrections modify the amplitude decompositions in fl3.27| )-( |3.29 ). In particular the 



AI = 5/2 contributions to all three decays have to be included. These contributions could have, 



in principle, an impact on e' /e, but their size is not understood at present [215|-pT 



5.9 Summary 

As of the beginning of the third millennium, we know confidently that there is a direct CP 
violation in Kl — > tttt decays. The ratio e'/e is measured to be around 2 • 10 -3 but values as 
low as 1 • 10 -3 are in view of the most recent NA48 results and the older E731 results not yet 
excluded. We should know the truth in the next few years. 

While the present theoretical status of e'/e is rather unsatisfactory, with various estimates 
ranging from 5 • 10" 4 to 4 • 10~ 3 , e'/e within the SM agrees within theoretical and experimental 
uncertainties with the data. The sign and the order of magnitude have been correctly predicted 
by theorists. Unfortunately, in view of very large hadronic and substantial parametric uncer- 
tainties, it is impossible at present to see what is precisely the impact of the e'/e-data on the 
CKM matrix. 

The short distance contributions to e'/e within the SM are fully under control. On the other 
hand considerable progress has to be made in the evaluation of the matrix elements {Qi)o,2- 
The calculations should not be confined to {Qq)o and {Qs)2 but should include matrix elements 
of all contributing operators. Personally, I believe that lattice calculations have eventually the 
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best chance to achieve this goal, but this may still take several years of extensive work [122]. On 
the other, it is very important to continue analytic efforts in order to confront the future lattice 
results. Finally isospin breaking corrections should also be taken into account. Here I expect 
analytic tools to be more powerful than lattice simulations for some time. 

In any case in view of very large theoretical uncertainties and sizable experimental errors 
there is still a lot of room for non-standard contributions to e' je. Indeed results from NA31, 
KTeV and NA48 prompted several analyses of e'/e within various extensions of the SM like 
general supersymmetric models [184, 219, 220.(221, 222], models with anomalous gauge couplings 
[223], four-generation models [224], left-right symmetric models [225] and models with additional 
fermions and gauge bosons 226 ]. Unfortunately several of these extensions have many free 
parameters and are not very conclusive at present. One should also remember that the hadronic 
uncertainties in the SM are also present in its extensions. The situation may change in the 
future when the calculations of hadronic matrix elements will improve and new data from high 
energy colliders will restrict the possible ranges of parameters involved. For instance a recent 
analysis of the bounds on anomalous gauge couplings from LEP2 data indicates that substantial 



enhancements of e'/e from this sector are very unlikely [227]. 

On the other hand from the point of view of Munich and Rome groups the e'/e data puts 
models in which there are new positive contributions to e and negative contibutions to e' in 
difficulties. In particular as analyzed in |173[| the two Higgs Doublet Model II can either be 
ruled out with improved hadronic matrix elements or a powerful lower bound on tan (5 can be 
obtained from e'/e. In the Minimal Supersymmetric SM (MSSM), in addition to charged Higgs 
exchanges in loop diagrams, also charginos contribute. For suitable choice of the supersymmetric 
parameters, the chargino contribution can enhance e'/e with respect to the SM expectations 
|111| [. Yet, as found in |111| , |109j [, the most conspicuous effect of minimal supersymmetry is a 
depletion of e'/e. We will quantify this in section 9. The situation can be different in more general 
models in which there are more parameters than in the two Higgs doublet model II and in the 
MSSM, in particular new CP violating phases. As an example, in general supersymmetric models 
e'/e can be considerably enhanced through the contributions of the chromomagnetic penguins 
[229|, |184|, gjl, 122(1], Z°-penguins with the opposite sign to the one in the SM pO, [228L pO, pi 



and isospin breaking effects in the squark sector [221]. 

The future of e'/e in the SM and in its extensions depends on the progress in the reduction of 
parametric and hadronic uncertainties. In any case e'/e already played a decisive role in estab- 
lishing direct CP violation in nature and its quite large value gives additional strong motivation 
for searching for this phenomenon in cleaner K decays like Ki — > ifivv and Ki — > 7r°e + e _ , in 
B decays, in D decays and elsewhere. We now turn to discuss some of these topics. 
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6 The Decays K + — ► ir + vv and — > n Q vv 
6.1 General Remarks 

We will now move to discuss the semileptonic rare FCNC transitions K + — > 7r + zvi/ and I^l — ► 
ir°vv. Within the SM these decays are loop-induced semileptonic FCNC processes determined 
only by Z°-penguin and box diagrams and are governed by the single function Xo(xt) given in 

(PI )- 



A particular and very important virtue of K —* 7ri^ is their clean theoretical character. This 
is related to the fact that the low energy hadronic matrix elements required are just the matrix 
elements of quark currents between hadron states, which can be extracted from the leading 
(non-rare) semileptonic decays. Other long-distance contributions are negligibly small p32[| . 
The contributions of higher dimensional operators are found to be negligibly small in the case of 

— > ix^vv [233] and below 5% of the charm contribution in the case of K + — > ■k+vv [234]. As 



a consequence of these features, the scale ambiguities, inherent to perturbative QCD, constitute 
the dominant theoretical uncertainties present in the analysis of these decays. These theoretical 
uncertainties have been considerably reduced through the inclusion of the next-to-leading QCD 
corrections |55|]-[|59|1. 

The investigation of these low energy rare decay processes in conjunction with their the- 
oretical cleanliness, allows to probe, albeit indirectly, high energy scales of the theory and in 
particular to measure Vtd and lm\ t = ImV^Vtd from K+ — ► tt + uv and Kj_, — ► iftvv respectively. 
Moreover, the combination of these two decays offers one of the cleanest measurements of sin 2/3 
p35| ]. However, the very fact that these processes are based on higher order electroweak ef- 
fects implies that their branching ratios are expected to be very small and not easy to access 
exp er iment ally. 

6.2 The Decay K + -> -k+vv 
6.2.1 The effective Hamiltonian 

The effective Hamiltonian for K + — > n + vv can be written as 

n& = % n " £ (v; s v cd x l NL + v t :v td x(x t ))(sd)v^A(^i)v-A. (e.i) 



y/2 2vr sin 2 9 W l= 



The index l=e, fj,, r denotes the lepton flavour. The dependence on the charged lepton mass 
resulting from the box-graph is negligible for the top contribution. In the charm sector this is 
the case only for the electron and the muon but not for the T-lepton. 
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The function X(xt) relevant for the top part is given by 



f-V.s 



X(x t ) = X (x t ) + —Xi(x t ) = rjx ■ X (x t ), 

47T 



j]x = 0.994, 



with the QCD correction [55]-|5 



X^xt) = Xi(st) + 8x-, 



dX {x t ) 
dx t 



Inx 



n ■ 



(6.2) 



(6.3) 



Here x M = fif/Myy with [i t = 0(m t ) and X\(x t ) is a complicated function given in |55]-|5q|. 
The //^-dependence of the last term in ( |6.3| ) cancels to the considered order the //t-dependence 
of the leading term Xo(xt(/i)). The leftover ^-dependence in X(xt) is below 1%. The factor r/x 
summarizes the NLO corrections represented by the second term in (|6,2j). With nit = mt(mt) 
the QCD factor rjx is practically independent of nit an d an d is very close to unity. 

The expression corresponding to X(xt) in the charm sector is the function X^ L . It results 



from the NLO calculation |5^] and is given explicitly in [58]. The inclusion of NLO corrections 
reduced considerably the large ^ c dependence (with [i c = 0{m c )) present in the leading order 
expressions for the charm contribution [|236j 1 . Varying /x c in the range 1 GeV < /x c < 3 GeV 
changes Xnl by roughly 24% after the inclusion of NLO corrections to be compared with 56% 
in the leading order. Further details can be found in |||, 5£]. The impact of the fi c uncertainties 
on the resulting branching ratio Br(K + — > 7r + z^) is discussed below. 

The numerical values for X^ L for /i = m c and several values of and m c (m c ) can be found 



in |58[| . The net effect of QCD corrections is to suppress the charm contribution by roughly 30%. 
For our purposes we need only 



Po(X) 



1 

T 1 



where the error results from the variation of and m c (m c ). The contribution of dimension 

MS v ' 



0.42 ± 0.06 



(6.4) 



eight operators, not included here, is estimated to be below the error in (6.4) |23 



6.2.2 Deriving the Branching Ratio 

The relevant hadronic matrix element of the weak current (sd)y-A in ( |6.1[ ) can be extracted 
with the help of isospin symmetry from the leading decay K + — > Consequently the 

resulting theoretical expression for the branching fraction Br{K + — > ix^vv) can be related to 
the experimentally well known quantity Br{K + — > 7r°e + ^). Let us demonstrate this. 
The effective Hamiltonian for the tree level decay K + — > 7r°e + ^ is given by 



7r°e + u) 



G F 

—^Vus( su )v-A(Vee)v-A 



(6.5) 
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Using isospin symmetry we have 

{ir + \{sd) V - A \K+) = V2(ir°\(su) v _ A \K + ). 



(6.6) 



Consequently neglecting differences in the phase space of these two decays, due to m n + ^ m n o 
and m e ^ 0, we find 



a 



Br(K+ -» w+uu) 
Br{K+ -»• Tr°e+v) ~ \V us \ 2 2tt 2 sin 4 6 W 



6.2.3 Basic Phenomenology 



]T \v; s v cd x^ L + v t * s v td x(x t ) 



l=e,fi,T 



(6.7) 



Using (6.7) and including isospin breaking corrections one finds 

2 



Br(K + tt+i/P) = k + • 



/ImA+ „. .\ z /ReA c „,^ ReAt „„. , 
(^^X(x t )J + ^^-£p (X) + ~^X(x t , 

3a 2 Br(K+ -► vr°e H 



«+ = o2 . 4n = 411 • 10_11 

Z7r z sin Bw 



where we have used 



o 



129' 



sin 2 9 W = 0.23, Br(K + 7r°e + i/) = 4.82 • 10 



-2 



(6.8) 
(6.9) 

(6.10) 



Here Aj = V^*Vi<i with A c being real to a very high accuracy. r K + = 0.901 summarizes isospin 
breaking corrections in relating K + — ► 7r + ^P to — > 7r°e + ^. These isospin breaking corrections 
are due to quark mass effects and electroweak radiative corrections and have been calculated in 
[ 237) 1. Finally P (X) is given in Q. 

Using the improved Wolfenstein parametrization and the approximate formulae ( 1.22j ) - 
( 1.24| ) we can next put ( |6.8| ) into a more transparent form [27]: 

1 



Br(K + -> Tr+z/i/) = 4.11 • 10~ u A i X 2 (x t )- {afjf + (go - qY 

a 



where 



a 



1- ^ 

. 2 < 



(6.11) 



(6.12) 



The measured value of Br(K + — ► 7r + z^) then determines an ellipse in the (g, 77) plane 
centered at (go>0) with 



go = 1 + 



and having the squared axes 



-2 2 
0i = r o> 



^ 2 X(x t ) 



-2 / r o 



(6.13) 
(6.14) 
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where 



1 



a ■ Br(K + — * tt + uv) 



(6.15) 



4.11 • 10- 11 

Note that ro depends only on the top contribution. The departure of go from unity measures 
the relative importance of the internal charm contributions. go ~ 1.35. 



The ellipse defined by ro, go an d cr given above intersects with the circle (1.32). This allows 
to determine g and fj with 

1 



oo - \Ve§ + C 1 " ^0 - ° 2 Rl) , n = ^Rl- Q 2 (6-16) 



= 1 9 
1 — <T 

and consequently 

R 2 t =l + R 2 b -2g, (6.17) 
where fj is assumed to be positive. Given g and fj one can determine Vtd'- 

V td = AX 3 (l-g- if)), \V td \ = AX 3 R t . (6.18) 

The determination of \Vtd\ and of the unitarity triangle requires the knowledge of (or A) and 
of |VubAc&|. Both values are subject to theoretical uncertainties present in the existing analyses 
of tree level decays. Whereas the dependence on |Vu&/^c&| is rather weak, the very strong 
dependence of Br(K + — * n + vv) on A or V c b makes a precise prediction for this branching ratio 
difficult at present. We will return to this below. The dependence of Br(K + — * ir + vT>) on m t is 
also strong. However mt is known already within ±4% and consequently the related uncertainty 
in Br{K + — > tt + vu) is substantialy smaller than the corresponding uncertainty due to V c b- 



6.2.4 Numerical Analysis of K + 



IT 1 VV 



The uncertainties in the prediction for Br(K + — > tt + uv) and in the determination of \Vtd\ 
related to the choice of the renormalization scales fit and [i c in the top part and the charm 
part, respectively have been investigated in M. To this end the scales \i c and fit entering 
m c (n c ) and mt(iM), respectively, have been varied in the ranges 1 GeV < \i c < 3GeV and 
100 GeV < fit < 300 GeV. It has been found that including the full next-to-leading corrections 
reduces the uncertainty in the determination of \Vtd\ from ±14% (LO) to ±4.6% (NLO). The 
main bulk of this theoretical error stems from the charm sector. In the case of Br(K + — > tt + ui>), 
the theoretical uncertainty due to fi c> t is reduced from ±22% (LO) to ±7% (NLO). 
Scanning the input parameters of table ^ we find 

Br{K + -> ir + vv) = (7.5 ± 2.9) • 10~ n (6.19) 

where the error comes dominantly from the uncertainties in the CKM parameters. 
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The predicted branching ratio decreased during the last years due to the improved lower 
bound on AM S . Indeed, it is possible to derive an upper bound on Br(K + — > ir + vv) 

2 



Br{K"< 



IT VV) 



K + 



P (X)+A 2 X(x t 



r_sd AM d 
A V AM, 



(6.20) 



where r ds = t,\J m B s / m B d with £ defined in ( |4.85| ) . This equation translates a lower bound on 
AM S into an upper bound on Br{K + — > ix^vv). This bound is very clean and does not involve 
theoretical hadronic uncertainties except for r S( i- Using 



AM, 



< 0.18 , A < 0.89 , P (X) < 0.48 , X(x t ) < 1.56 , r sd < 1.22 (6.21) 



we find 



11.5 • 10 



-ii 



(6.22) 



This limit could be further strengthened with improved input. However, this bound is strong 
enough to indicate a clear conflict with the SM if Br{K + —* ir + vv) should be measured at 
1.5 • 10~ 10 . 



6.2.5 \Vtd\ from 



Once Br(K+ 
uncertainties: 



a (Vf 



i: + vv) = Br(K + ) is measured, \V t d\ can be extracted subject to various 

a{Br{K+)) 



td\ 



±0.04 sca i e ± 



cb\ 



±0.7- 



a{m r 



±0.65- 



(6.23) 



\V td \ \Vc\ fh c Br(K+) 

Taking cr(|V^»|) = 0.002, a(fh c ) = 100 MeV and a(Br(K + )) = 10% and adding the errors in 
quadrature we find that \Vtd\ can be determined with an accuracy of ±10%. This number is 
increased to ±11% once the uncertainties due to m t , a s and |Ki&|/|^cb| are taken into account. 
Clearly this determination can be improved although a determination of \ Vtd\ with an accuracy 
better than ±5% seems rather unrealistic. 



6.2.6 Summary and Outlook 

The accuracy of the SM prediction for Br{K + — » ir + vi>) has improved considerably during the 
last decade. This progress can be traced back to the improved values of mt and \V c b\ and to 
the inclusion of NLO QCD corrections which considerably reduced the scale uncertainties in the 
charm sector. 

Now, what about the experimental status of this decay ? One of the high-lights of 1997 was 



the observation by AGS E787 collaboration at Brookhaven 238 ] of one event consistent with the 
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signature expected for this decay. Recently an analysis of a larger data set has been published 



[236]. No further events were seen, leading to the branching ratio: 

Br(K + -» n+vD) = (1.5±H) ■ 10~ 10 . (6.24) 

The central value is by a factor of 2 above the SM expectation but in view of large errors the 
result is compatible with the SM. The analysis of additional data on K + — ► k^vv present on 
tape at AGS E787 should improve this result in the near future considerably. In view of the 
clean character of this decay a measurement of its branching ratio at the level of 1.5 • 10 -10 with 
a small error would signal the presence of physics beyond the SM. 

The experimental outlook for this decay has been recently reviewed by Littenberg p40| ], 
A new experiment, AGS E949 [241] is scheduled to run in 2001-3. It is expected to reach a 
sensitivity of ~ 10 _11 /event. In the long term, the CKM experiment at Fermilab |242f| should be 
able to reach ~ 10~ 12 /event sensitivity. At this level an accurate measurement of the branching 
ratio should be possible. 

6.3 The Decay K L -> tt uu 
6.3.1 The effective Hamiltonian 

The effective Hamiltonian for — > ir°vv is given as follows: 

H ea = -| . 2 V*V td X(x t )(sd) V - A (vv)v-A + h.c. , (6.25) 
a/2 2ir sin Bw 

where the function X{xt), present already in K + — * ix^vv, includes NLO corrections and is 
given in ( |6.2| ). 

As we will demonstrate shortly, — > ifivv proceeds in the SM almost entirely through 
direct CP violation ]|243|| . The indirectly CP violating contribution and the CP conserving 
contribution analyzed in ] [233|] are fully negligible. Within the SM — ► tx^vv is completely 
dominated by short-distance loop diagrams with top quark exchanges. The charm contribution 
can be fully neglected and the theoretical uncertainties present in K + — * tt + vi> due to m c , 
fj, c and A.jfg are absent here. Consequently the rare decay — ► 7r°^P is even cleaner than 
K + — > tt + vi> and is very well suited for the determination of the Wolfenstein parameter rj and 
in particular ImAt- 

We have stated that the decay — » ifivv is dominated by direct CP violation. Now as 
discussed in Section 3 the standard definition of the direct CP violation requires the presence of 
strong phases which are absent in — > tt uu. Consequently the violation of CP symmetry in 

— ► ifivv can only arise through the interference between K° — K° mixing and the decay am- 
plitude. This type of CP violation has been discussed already in Section 3. However, as already 
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pointed out by Littenberg [243] and demonstrated explictly in a moment, the contribution of CP 
violation to — > ifivv via mixing alone is tiny. It gives Br(K^ — > tt^vv) ~ 2 • 10 15 

and its interference with the directly CP-violating contribution is 0(1O~ 13 ). Consequently, in 
this sence, CP violation in — > iftvv with Br(Ki, — ► tv°uv) = O(10 -11 ) is a manifestation of 
CP violation in the decay and as such deserves the name of direct CP violation. In other words 
the difference in the magnitude of CP violation in — ► tttt (e) and Ki, — > 7r°z/P is a signal of 
direct CP violation and measuring — > 7r°z^i/ at the expected level would be another signal of 



this phenomenon. More details on this issue can be found in [244, 245, 246]. 
6.3.2 Deriving the Branching Ratio 

Let us derive the basic formula for Br^K^ — ► ir®vv) in a manner analogous to the one for 
Br{K + — ► 7r + z^). To this end we consider one neutrino flavour and define the complex function 

F = ^ % V t * s V td X{x t ). (6.26) 

V2 27rsin 2 Gw 

Then the effective Hamiltonian in fl6.25j ) can be written as 

W eff = F(sd) y _ A (z^)y-A + F*(ds) v -A(vv)v-A ■ (6.27) 
Now, from fl3.8|) we have 

K L = -J=[(l + ^ + (l- e -)K°] (6.28) 

where we have neglected | e 1. Thus the amplitude for Ki —* 7r°z^P is given by 



A(K L -» vr'W) = F(l +e)(7r u |(sd)y_ A |^ U ) + - e)(7r u |(J S )y_ A |^ u ) (^)y-A 



1 



V2 

Recalling 
we have 



Ol/sA. .I^0\ , it*/-, =-\/„0|/J„N . . |&>0\ 



(6.29) 

CP|K°) = C|Er°) = (6.30) 

(7r°|(d S )y_A|^°) = -(Tr \(sd)v-A\K°), (6.31) 



where the minus sign is crucial for the subsequent steps. 
Thus we can write 

A{K L - vr ^) = i= [F(l + e)-F*(l-e)] (vr |(sd)y_A|^ )(Hy-A. (6.32) 

Now the terms e can be safely neglected in comparision with unity, which implies that the 
indirect CP violation (CP violation in the K° — K° mixing) is negligible in this decay. We have 
then 

F(l + e) -F*(l- s) = % a lra{V t * s V td ) ■ X(x t ). (6.33) 

V2 7T sm z 9w 
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Consequently using isospin relation 



(Aids] 



V-A 



\K 



(7r°\(su)v-A\K" 



(6.34) 



together with ( |6.5[ ) and taking into account the difference in the lifetimes of Kl and K + we 
have after summation over three neutrino flavours 

„2 



Br{K L ^^vv) r(K L 



a 



t{K+) \V us \ 2 2tt 2 sin 4 G w 



Br(K+ 
where X t = V t *V u . 

6.3.3 Master Formulae for Br(K L -» ir°vv) 
Using ( |6.35| ) we can write Br^K-^ — > -k^vv) simply as follows 

flmX t 



[ImXfX(x t )Y 



(6.35) 



Br(K L 



7T°^Z/) = Kl 



V A* 

= 1.80 • 10~ 10 



(6.36) 
(6.37) 



r K + t{K^ 

with k + given in fl6,9| ) and r/^ L = 0.944 summarizing isospin breaking corrections in relating 
K L -> vr ^ to if+ -> 7r°e+^ |p37 |- 

Using the Wolfenstein parametrization and (|6.2| ) we can rewrite ( |6.36 ) as 

4 

(6.38) 



Br(K^ -> vr ^) = 3.0 • 10~ n 



0.39 



m t (m t ) 



170 GeU 



eft I 

0.040 



The determination of rj using Br(K^ — ► ifivv) requires the knowledge of U c b and mt. The 
very strong dependence on or ^4 makes a precise prediction for this branching ratio difficult 
at present. 

On the other hand inverting ( |6.36 ) and using ( |6,2[ ) one finds [246]: 



ImA, = 1.36 • 10~ 4 



"170 GeV" 


1.15 


~Br(K L -» tt°vP) 


. m t (m t ) _ 




3 • 10- 11 



.1/2 



(6.39) 



without any uncertainty in \V c b\- ( 6.3S| ) offers the cleanest method to measure ImAt; even better 
than the CP asymmetries in B decays discussed in section 8 that require the knowledge of |V^,| 
to determine ImAt. Measuring Br(K^ — > -k^uu) with 10% accuraccy allows to determine ImAf 
with an error of 



IS 



69 



6.3.4 Numerical Analysis of — ► it vv 

The //^-uncertainties present in the function X(xt) have already been discussed in connection 
with K + — ► tt + vv. After the inclusion of NLO corrections they are so small that they can be 
neglected for all practical purposes. Scanning the input parameters of table ^ we find 

Br(K L n°v9) = (2.6 ± 1.2) • 10~ n (6.40) 

where the error comes dominantly from the uncertainties in the CKM parameters. 



6.3.5 Summary and Outlook 

The accuracy of the SM prediction for Br^K^ — > tt°uu) has improved considerably during the 
last decade. This progress can be traced back mainly to the improved values of m t and \V c b\ 
and to some extent to the inclusion of NLO QCD corrections. 



The present upper bound on Br(K^ — ► tt^vu) from the KTeV experiment at Fermilab [247] 
reads 

Br(K L -> n°uu) < 5.9 • 10 -7 . (6.41) 
This is about four orders of magnitude above the SM expectation (6.40). Moreover this bound 



is substantially weaker than the model independent bound [244] from isospin symmetry: 

Br(K L -» ir°i/v) < 4.4 • Br(K + -> tt + vu) (6.42) 

which through ( |6.24 ) gives 

Br(K L -» vr°zyP) < 2.6 • 10~ 9 (90%C.L.) (6.43) 



The experimental outlook for this decay has been recently reviewed by Littenberg [240]. The 
KEK E391a experiment [248] should reach sensitivity of ~ 10~ 10 /event which would give some 



events only in the presence of new physics contributions (see Section 9). KAMI at Fermilab [24£] 
should be able to reach a sensitivity of < 10~ 12 /event. Finally a very interesting new experiment 
KOPIO at Brookhaven (BNL E926) p50(] expects to reach the single event sensitivity of 2- 10~ 12 . 
Both KAMI and KOPIO should provide useful measurements of this gold-plated decay. 



6.4 Unitarity Triangle and sin 2(3 from K — > txvv 

The measurement of Br(K + — > -k + vv) and Br(K^ — > 7r°i^) can determine the unitarity trian- 
gle completely, (see fig. |l"0"|) , provided m t and V^, are known |235| . Using these two branching 
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(0,0) 



(1,0) 



Figure 10: Unitarity triangle from K — > ttuv. 



ratios simultaneously allows to eliminate |Vuft/V^| from the analysis which removes a consider- 
able uncertainty. Indeed it is evident from ( |6.8[) and ( 6.36| ) that, given Br(K + — > -k + vv) and 
Br(K^ — > ir°ui>), one can extract both ImAf and ReA(. One finds [ 235 , [§ 



ImAt = A L 



X(x t ) 



ReAj = -A 



where we have defined the "reduced" branching ratios 



Br{K + tt+vu) Br{K h iPvv) 

Bl = TT-t — -,»_n -°2 =" 



4.11 -10- 11 1.80 -10- 10 

Using next the expressions for ImAt, ReAt and ReA c given in ( 1.22j )- (]l,24| ) we find 

v 7 ^ 



Q = 1 + 



P (X) - vHBi - B 2 ) 



^4 2 X(x t ) ' 77 ^A(x ( ) 

with a defined in ( 6.12 ). An exact treatment of the CKM matrix shows that the formulae ( |6.46| ) 
are rather precise [235], 

Using ( |6.46|) one finds subsequently [ 235 1 

2r, 



(6.44) 



(6.45) 



(6.46) 



r s = r s (Bi, B 2 ) = _ g = cot (5 , sin 2/3 



1 + rl 



with 



r s {B u B 2 ) = 



s /a(B 1 - B 2 ) - P (X) 



(6.47) 



(6.48) 



Thus within the approximation of (|6.46D sin 2/3 is independent of V c b (or A) and mj. 
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Table 10: Illustrative example of the determination of CKM parameters from K — > irvv. 





a(\V cb \) = ±0.002 


a(\V cb \) = ±0.001 


*{\Vtd\) 


±10% 


±9% 


a{g) 


±0.16 


±0.11 


a(fj) 


±0.04 


±0.03 


a(sm 2(3) 


±0.05 


±0.05 


o-(ImAj) 


±5% 


±5% 



It should be stressed that sin 2(3 determined this way depends only on two measurable branch- 
ing ratios and on the function Pq(X) which is completely calculable in perturbation theory. 
Consequently this determination is free from any hadronic uncertainties and its accuracy can be 
estimated with a high degree of confidence. 



An extensive numerical analysis of the formulae above has been presented in [246]. Assuming 
that the branching ratios are known to within ±10% and mt within ±3 GeV one finds the results 
in table 10 [246]. We observe that respectable determinations of all considered quantities except 
for q can be obtained. Of particular interest are the accurate determinations of sin 2(3 and 
of ImAf. The latter quantity as seen in ( 6.39| ) can be obtained from — > ifivv alone and 
does not require knowledge of V cb - The importance of measuring accurately ImAj is evident. 
It plays a central role in the phenomenology of CP violation in K decays and is furthermore 
equivalent to the Jarlskog parameter Jcp the invariant measure of CP violation in the SM, 
J CP = A(l - A 2 /2)ImA t . 

The accuracy to which sin 2(3 can be obtained from K — > ixvv is, in the example discussed 
above, comparable to the one expected in determining sin 2(3 from CP asymmetries in B decays 
prior to LHCB and BTeV experiments. In this case sin 2(3 is determined best by measuring 
CP violation in Bj — ► J/tpK$. Using the formula for the corresponding time-integrated CP 
asymmetry one finds an interesting connection between rare K decays and B physics J235|] 



.2 



2rs{Bl > B2) -a^B^J/^Ks) 1 -^ (6.49) 



l + r*(B u B 2 ) — ~ a x d 

which must be satisfied in the SM. Here x d = AMd/T(B®) is a B® — B® mixing parameter. We 
stress that except for Pq(X) all quantities in ( 6.4S| ) can be directly measured in experiment and 
that this relationship is essentially independent of mt and V cb . Due to very small theoretical 



uncertainties in (6.49), this relation is particularly suited for tests of CP violation in the SM and 
offers a powerful tool to probe the physics beyond it. We will return to this topic in Section 9. 



72 



The relation ( |6.49| ) is one of several correlations between K — > ttvu and observables in 
B physics. Another example is the upper bound in (|6.20| ), A recent numerical analysis of 



such correlations can be found in [251 1. Finally analyses of K — > ixvv in non-super symmetric 
extensions of the SM can be found in 252 1. K — > -nvv in supersymmetric extensions is discussed 
in section 9. 

7 Express Review of Rare K and B Decays 

7.1 Rare K Decays 
7.1.1 K L -> ir°e + e- 

There are three contributions to this decay: CP-conserving, indirectly CP-violating and directly 
CP-violating. Unfortunately out of these three contributions only the directly CP-violating can 
be calculated reliably. This contribution is dominated by Z°-penguin diagrams. The enhance- 



ment of these diagrams for large mt, pointed out in the context of this decay in [253], gave in 



fact the motivation to study Z°-penguin effects in e'/e [163, 164]. Including NLO corrections 



31| and scanning the input parameters of table |2j we find 

Br(K L -► 7r°e + e-) dir = (4.3 ± 2.1) • 1(T 12 , (7.1) 
where the errors come dominantly from the uncertainties in the CKM parameters. The calcula- 



tions of indirectly CP-violating contribution are plagued by theoretical uncertainties [254]. On 
the other hand, this contribution taken alone is given simply by 

Br(K L - A+e-) indir = \e\ 2 ^^-Br(K s - vr e + e -) = 3.0 • l^Br(K s -» ^e+e") , (7.2) 

t{K S ) 

with t^Kl^s) denoting the K^s life-times and Br(K$ — ► 7r°e + e~) hopefully measured in the 
next years by KLOE at Frascati. The two CP violating contributions will in general interfer 
with each other. Given the present uncertainty on Br{Ks — > 7r°e + e _ ) the total CP-violating 
contribution could be as high as few x 10~ n [255] but taking into account the theoretical 
estimates of the indirectly CP-violating contribution, one should expect it below 10~ n within 
the SM. 

The upper bound on the CP-conserving contribution governed by Kl — ► 7r°77 — » 7r°e + e _ 
can be obtained with the help of chiral perturbation theory |]254|] and the data on Kl — > 7r°77. 
The recent results on the latter decay [)256| imply that this contribution is smaller than 2- 10 -12 . 
Consequently it is smaller than expected by some authors in the past. As this contribution 
does not interfere with the remaining larger CP-violating contributions, it does not present a 
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significant problem but in order to be able to extract CKM parameters from this decay its better 
estimate is clearly needed. 

The most recent experimental bound from KTeV f257f | reads 

Br{K L -> 7r°e+e-) < 5.1 • 1(T 10 (90%C.L.). (7.3) 

Considerable improvements are expected in the coming years. 

7.1.2 K L n+n- 

The K]_, — > branching ratio can be decomposed generally as follows: 

BR(K L -> n + fi~) = \ReA\ 2 + |W| 2 , (7.4) 

where Re^4 denotes the dispersive contribution and ImA the absorptive one. The latter contri- 
bution can be determined in a model independent way from the Kl — > 77 branching ratio. The 
resulting \ImA\ 2 = (7.07 ± 0.18) • 10~ 9 is very close to the experimental branching ratio p5q] 

Br(K L -» n+fT) = (7.18 ± 0.17) • 10~ 9 , (7.5) 

so that |Re^4| 2 is substantially smaller and extracted to be p5Sf | 

|ReA exp | 2 < 3.7- 10~ 10 (90% C.L.). (7.6) 

Now ReA can be decomposed as 

Re.4 = Re^LD + ReA S D , (7.7) 

with 

|ReylsD| 2 = Br(K h m+zOsd (7.8) 

representing the short-distance contribution which can be calculated reliably. An improved 
estimate of the long-distance contribution Re^LD has been presented in ||259[] 

\ReA LD \ < 2.9 • 10" 5 (90% C.L.). (7.9) 

Together with ( [7.6|) this gives 

Br(K L -» m + /Osd < 2.3 • 10~ 9 . (7.10) 



This result is close to the one presented by Gomez Dumm and Pich 260] . More pesimistic view 
on the extraction of the short distance part from Br(K\^ — ► /i + /i~) can be found in [261]. 

The bound in ( 7.10| ) should be compared with the short distance contribution within the SM 
for which we find 

Br(K L -> ^+/Osd = (8.7 ± 3.7) • lO" 10 . (7.11) 

This implies that there is a considerable room for new physics contributions. Reviews of rare K 
decays are given in |], £| [240j [262 ]. 
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7.2 Rare B Decays 
7.2.1 Preliminaries 



The most interesting here are the decays B — * X s ^, B — > X s ^l + l~ , B — * X s ^vis, B s d — > 
and the exclusive counterparts of the first three decays with X s d replaced by K*(g) or K(tt). 
Within the SM the dominant operators are 

Qi = s L (T^m b b R F^, Q 9 = SL^bJ^l , (7.12) 

Qio = SLl^hh^l, Qv = s L l^b L v^^(l - . (7.13) 

B — > X s ^ is governed by Q-j and B — > X s ^l + l~ by Q7, Q9 an d Qio- -6 — > X s ^vv and 
-Bs.rf - ► involve only Q„ and Qio respectively. 

There is a vast literature on these decays within the SM and its extensions. I can recommend 
in particular the most recent analyses and reviews in p6. 
where further references can be found. Theoretically cleanest are the decays B — » X s ^vv and 
B s ,d ~^ as they involve only Z°-penguin and box diagram contributions. B — > Xgd'y and 
13 — > X S) dl + l~ involving magnetic 7-penguins and the standard photon-penguins are subject to 
long distance contributions which have to be taken into account. In the case of exclusive channels 
there are additional uncertainties in the relevant hadronic formfactors. On the experimental side 
the exclusive channels are easier to measure. Moreover, knowing them experimentally will help 
to distinguish between B — * X s l + l~ and B — * X^l + l^ . 



264, 



268, 26£, W! 



7.2.2 B 



A lot of efforts have been put into predicting the branching ratio for the inclusive radiative 
decay B — > X s ^f including NLO-QCD corrections and higher order electroweak corrections. The 
relevant references can be found in table |] and [jl], 263, 264 1, where theoretical details are given. 
The final result of these efforts can be summarized by [91] 



Br(B -> X,7) th = (3.29 ± 0.21 ± 0.21) • 10" 



(7.14) 



The two errors correspond to parametric and scale uncertainties. The main theoretical achieve- 
ment was the reduction of the scale dependence through NLO calculations, in particular those 



given in [67| and [47]. In the leading order this uncertainty alone would induce an error of ±0.6 



W3- 



The theoretical result in ( [7.14 ) should be compared with experimental data from CLEO 

whose combined branching ratio reads 



| 2731 , ALEPH j27| and BELLE []27 

Br(B X s7 ) cxp 



(3.21 ± 0.40) • 10" 



(7.15) 
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Clearly, the SM result agrees well with the data. In order to see whether any new physics can 
be seen in this decay, the theoretical and in particular experimental errors should be reduced. 
This is certainly a difficult task. On the other hand already the available data put powerful 
constraints on the parameter space of the supersymmetric extensions of the SM. Most recent 



analyses can be found in | 10£ , 276 , 277]. 

It is easier to measure the exclusive branching ratios Br(B^ — > i^* ± 7) and Br(B^ — > K*°j). 
The most recent data from CLEO, BaBar and Belle can be found in |278f |. They are found in 
the ball park of 3 • 10~ 5 and 5 • 10 -5 respectively. Theoretical calculations of these branching 
ratios are plagued by uncertainties in the relevant formfactors. Most recent theoretical reviews 



are given in 279 1. 



7.2.3 B 



XJ+l~ 



and B -» K*l + l 



These rare decays have been the subject of many theoretical studies. The NLO-QCD corrections 
have been calculated in [B3, Recently dominant NNLO corrections to the relevant Wilson 



coefficients have been calculated in [81]. The remaining scale dependence in the branching ratios 
is roughly 13%. In order to remove this dependence two-loop matrix elements of the relevant 
four-quark operators have to be calculated. 

The main interest in these decays is their large sensitivity to new physics contributions 
that are encoded in the Wilson coefficients CV, C9 and C\q of the operators in ( 7.12] ) and 
( 7.13| ). These contributions can be studied through the dilepton mass distribution, the leptonic 
forward-backward asymmetry and CP-asymmetries. Most recent analyses can be found in [265, 
|266| , 267 , 268 , 269 , 270] . While B — > X s ^ is only sensitive to the absolute value of C7, the 
decays B — > X s l + l~ and B — > K*l + l~ depend also on the sign of this coefficient that in more 
complicated models can be opposite to the one in the SM. In addition these decays involve the 
coefficients Cg and C10 that in the extensions of the SM can take rather different values than 
found in the SM. Of particular interest is the leptonic forward-backward asymmetry (FBA) 
in B — > K*fi + fi~ and B — * K*e + e~ which vanishes for a certain value of the dilepton mass. 
The position of this zero is a sensitive function of the ratio C7 /Re(Cg) with small uncertainties 



from the formfactors [280, 281, |266| [. Moreover FBA is proportional to C±q and consequently 
sensitive to the magnitude and the sign of this coefficient. It is clear that this asymmetry will 
offer useful tests of the SM and of its extensions. Similar comments apply to CP asymmetries 
that are very small in the SM but can be substantial in general supersymmetric models. The 
calculations of and 1/m^ can be found in 282 , |272| [ and 283 , 284 ], respectively. In the 

process of including non-perturbative corrections induced by intermediate cc states one has to 
avoid double counting. The most efficient method here at present is the dispersion approach of 



76 



Kriiger and Sehgal [285], which makes use of experimental data. 



Within the SM we have §6§ HI 

Br(B -> X s // + /0 = (5.7 ± 1.2) • 10~ 6 , Br(B -> #V>~) = (2.0 ± 0.7) • 10~ 6 , (7.16) 
to be compared with the experimental upper bounds 

£r(£ ->X„n + ir) < 4.2- 10" 5 [286], Br{B -> K*n + fi~) < 4.0 • 10~ 6 p87fl . (7.17) 
The first events for these decays are expected from B factories and Tevatron already this year. 



7.2.4 B -► 

The decays 5 — ► X s ^yv are the theoretically cleanest decays in the field of rare B-decays. They 
are dominated by the same -penguin and box diagrams involving top quark exchanges which 
we encountered already in the case of K + — * k^vv and — > tt^vv, except for the appropriate 
change of the external quark flavours. Since the change of external quark flavours has no impact 
on the m t dependence, the latter is fully described by the function X(xt) in Q6.2p which includes 
the NLO corrections. The charm contribution is fully negligible here and the resulting effective 
Hamiltonian is very similar to the one for — ► 7T°z>t> given in Q6.25| ). For the decay B — ► X s vv 
it reads 



n 



a 



cff 



-Vt b Vt s X(xt){bs)v-A{vv)v-A + h.c. 



(7.18) 



^2 2vr sin 2 9 W 
with s replaced by d in the case of B — ► X^vv. 

The theoretical uncertainties related to the renormalization scale dependence are as in 

Kl — > 7r°^P and can be essentially neglected. The same applies to long distance contributions 

considered in [283]. The calculation of the branching fractions for B — > X S( ivv can be done in 

the spectator model corrected for short distance QCD effects. Normalizing to Br(B — > X c ev) 

and summing over three neutrino flavours one finds 



Br(B 



X s vu) 



3a 2 



\V ts \ 2 X 2 (x t ) «(0) 



Br{B^X c ev) 4vr 2 sin 4 9 W I V cb \ 2 f(z) 



K(Z 



(7.19) 



Here f(z) is the phase-space factor for B — » X c eP with z = m^/m 2 and k(z) = 0.88 |288] , [289|| is 
the corresponding QCD correction. The factor k(0) = 0.83 represents the QCD correction to the 
matrix element of the b — > svv transition due to virtual and bremsstrahlung contributions. In 
the case of B — > X^uv one has to replace Vt s by Vtd which results in a decrease of the branching 
ratio by roughly an order of magnitude. 

Setting Br{B -» X c ev) = 10.4%, f(z) = 0.54, k(z) = 0.88 and using the values in ( |Q)D we 
have 



Br(B 



X s vv) 



3.7-10" 



V, 



ts 



|2 r 



\v*\< 



m t lm t 



170 GeV 



2.30 



(7.20) 
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Taking next, f(z) = 0.54 ± 0.04 and Br(B — ► X c ev) = (10.4 ± 0.4)% and scanning the input 
parameters of table |2| we find 

Br(B -» X s vv) = (3.5 ± 0.7) • 10~ 5 (7.21) 
to be compared with the experimental upper bound: 

Br(B -> X s vv) < 7.7 • 10~ 4 (90% C.L.) (7.22) 



obtained for the first time by ALEPH |290[| . This is only a factor of 20 above the SM expectation. 
Even if the actual measurement of this decay is very difficult, all efforts should be made to reach 
this goal. One should also make attempts to measure Br{B — > Xdyv). Indeed 

Br{B^X d vv) _ \V t tf 

Br{B^X s vv) |Vto| 2 1 ' ' 

offers the cleanest direct determination of |Vtd|/|Vf S | as all uncertainties related to mt, f(z) and 
Br{B — > X c ev) cancel out. 

It might be much easier to measure the exclusive mode B — > K*vv. Most recent discussions 



can be found in [291] 



7.2.5 B 8td -> l+r 

The decays B s ^ — > are after B — ► X s ^vv the theoretically cleanest decays in the field of 
rare S-decays. They are dominated by the Z°-penguin and box diagrams involving top quark 
exchanges which we encountered already in the case of B — ► X s ^vv except that due to charged 
leptons in the final state the charge flow in the internal lepton line present in the box diagram is 
reversed. This results in a different mt dependence summarized by the function Y(xt), the NLO 
generalization |56| , 57, 5S] of the function Yo(xt) given in ( |2.22|) . The charm contributions are 



fully negligible here and the resulting effective Hamiltonian is given for B s — * l + l as follows: 

H ea = -- 1 - , 2 V t tV ta Y(x t )(bs) V - A (ll)v-A + h.c (7.24) 
V2 27rsm 9w 

with s replaced by d in the case of B d — > . 
The function Y{xt) is given by 

Y(x t ) = Y {x t ) + -^Yi(xi) = rfyYoixt), r] Y = 1.012 (7.25) 

47T 



where Y\(xt) can be found in |56|, |5^, ^8[. The leftover /^-dependence in Y(xt) is tiny and 
amounts to an uncertainty of ±1% at the level of the branching ratio. With mt = m t (mt) the 
QCD factor r\y depends only very weakly on m%- The dependence on A-p^ can be neglected. 
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The branching ratio for B s — > is given by |p6| 



Br(B, 



i + i-) = t(b s ) 



a 



it V47rsin 0w 



F^ s mfm Bs 



1 



mr 



m. 



V t tV ts \'Y 2 (x t ) (7.26) 



where B s denotes the flavour eigenstate (6s) and Fg s is the corresponding decay constant. Using 
( 6.10| ) and ( |7.25| ) we find in the case of B s — > fi + fj,~ 



Br(B s 



3.3 • 10 -9 



r(B s 



1.5 ps 



F, 



\v t 



0.040 



m t (m t ) 
170 GeV 



3.12 



(7.27) 



_210MeV. 

The main uncertainty in this branching ratio results from the uncertainty in Fb s ■ Scanning 
the input parameters of table |2| together with t(B s ) = 1.5 ps and Fb„ = (210 ± 30) MeV we find 



Br(B s 



(3.2 ± 1.5) • 10 



-9 



(7.28) 



For B d — ► fM + fi~ a similar formula holds with obvious replacements of labels (s — > d). 
Provided the decay constants Fb 3 and F Bd will have been calculated reliably by non-perturbative 
methods or measured in leading leptonic decays one day, the rare processes B s — > fi + ur and 
Bd — * should offer clean determinations of \Vt s \ and \Vtd\- In particular the ratio 

Br(B d -+ n+jT) _ r{B d ) m Bd B% \V td \ 2 

Br(B s ^ /j+fi-) T(B s )m Bs Fl s \V ts \ 2 K ' ' 

having smaller theoretical uncertainties than the separate branching ratios should offer a useful 
measurement of |T4d|/|^fcs|- Since Br{B d — > u, + u,~) = O(10 -10 ) this is, however, a very difficult 
task. For B s — > t + t~ and B s — > e + e~ one expects branching ratios 0(1O -6 ) and 0(1O -13 ), 
respectively, with the corresponding branching ratios for i^-decays by one order of magnitude 
smaller. 

The bounds on B s ^ d — > ^ip, are still many orders of magnitude away from SM expectations. 
The 90%C.L. bounds from CDF read 



Br(B s -» ix + n~~) < 2.0- 10" 



(90%C.L.) [292] 



(7.30) 



and Br(B d -> < 6.8- 10" 7 . CLEO |9f] provided the bound Br(B d fj+ n~) < 6.1- 10~ 7 . 

CDF should reach in Run II the sensitivity of 1 • 10 -8 and 4 • 10 -8 for B d — ► u.jl and B s — > fj,£L, 
respectively. Thus if the SM is the whole story one will have to wait until LHC-B and BTeV to 
see any events. On the other hand in a Two-Higgs-Doublet-Model and in particular in the MSSM 



one can find substantially larger branching ratios provided tan (3 is large [294]. This means that 
either this decay will be measured in Run II at Fermilab or the allowed parameter space in these 
models will be considerably reduced. The usefulness of this decay and of B — * t + t~ in tests of 



the physics beyond the SM is discussed in these papers and in [295]. 
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8 CP Violation in B Decays 



8.1 Preliminaries 

CP violation in B decays is certainly one of the most important targets of B-factories and of 
dedicated B-experiments at hadron facilities. It is well known that CP-violating effects are 
expected to occur in a large number of channels at a level attainable experimentally in the near 
future. Moreover there exist channels which offer the determination of CKM phases essentially 
without any hadronic uncertainties. Therefore, it is expected that during the coming years very 
important progress in our understanding of CP violation will come through the measurements 
of CP asymmetries in B decays as well as through various strategies for extracting the angles 
a, (3 and 7 from two-body B decays. They are extensively discussed in the books |], [7|], in the 
working group reports in (8|, |j| and in |||, 

Here after recalling some useful formulae, I will confine the discussion to three topics: 

• An express review of the classic strategies for the determination of the angles of the 
unitarity triangle, 

• Short description of recent strategies for the determination of the angle 7 from the decays 
B -> ttK, 

• Strategies involving U-spin symmetry 



131, 296, E97 



8.2 A Few Useful Formulae and Examples 

Let us begin our discussion with neutral B decays to CP Eigenstates. A time dependent asym- 
metry acp(t, f) in the decay B° — > / with / being a CP eigenstate is given in ( 3.69| ) where 



we have separated the direct CP-violating contributions from those describing CP violation in 
the interference of mixing and decay. As demonstrated in Section 3.7, an interesting case arises 
when a single mechanism dominates the decay amplitude or the contributing mechanisms have 
the same weak phases. Then acp(t, f) is given simply by 

a C p{t, f) = V f sm(2<t> D - 2cj> M ) sin(AMt) (8.1) 

where (pp is the weak phase in the decay amplitude and 4>m the weak phase in the — B® s 
mixing, rjt = ±1 is the CP parity of the final state. In this particular case the hadronic matrix 
elements drop out, the direct CP-violating contribution vanishes and the CP asymmetry is given 
entirely in terms of the weak phases (po an d 4>m- 
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If a single tree diagram dominates, the factor sin(2</>£> — 2(f) m) can be calculated by using 



7 b^u f -0 B Q d 



b^c -p a B° 



i.2) 



where we have indicated the basic transition of the b-quark into a lighter quark. (3 S = 0(10 2 ) 
On the other hand if the penguin diagram with internal top exchange dominates one has 




".w= { ^ !l (8-3) 



-Ps B, 



Let us practice these formulae. Assuming that B^ — > i/jKs and B^ — > ir + n are dominated by 
tree diagrams with b — > c and b — > it transitions respectively we readly find 

a C p(t, ^5) = - sin(2/3) sin(AM d t), (8.4) 

acp(vr + 7r") = - sin(2a) sin(AM d i). (8.5) 

Now in the case of B^ — » ^Jfg the penguin diagrams have to a very good approximation the 
same phase {(po = 0) as the tree contribution and are moreover Zweig suppressed. Consequently 



.4) is very accurate. This is not the case for B^ — > 7r + 7r _ where the penguin contribution could 
be substantial. Heaving weak phase 4>D = which differs from the tree phase (f>u = 7, this 
penguin contribution changes effectively sin(2a) to sin(2a + #p) where Op is a function of [3 and 
hadronic parameters. Strategies to determine Op and consequently a are discussed below. 

Similarly the pure penguin dominated decay B& — * 4>Ks is governed by the b — > s penguin 
with internal top exchange which implies that in this decay the angle (3 is measured. The 
accuracy of this measurement is a bit lower than using B^ — > ipKs as penguins with internal u 
and c exchanges may introduce a small pollution. 

Finally we can consider the asymmetry in B s — > t/jcf), an analog of Bj — > ipK s . In the leading 
order of the Wolfenstein parametrization the asymmetry acp(t,i/j(f)) vanishes. Including higher 
order terms in A one finds [|29£ 



a C p(t,ip((>) = 2X 2 rjsm(AM s t) (8.6) 
where A and r\ are the Wolfenstein parameters. 

8.3 Classic Strategies 
8.3.1 The Angle a 

The classic determination of a by means of the time dependent CP asymmetry in the decay 
B® — > 7r + 7r~ as given by ( |S.5D is affected by the "QCD penguin pollution" which has to be taken 
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care of in order to extract a. We have just mentioned this problem. The CLEO results for 



penguin dominated B — * ttK decays indicated 1 300 1 that this pollution could be substantial as 



stressed in particular in [pOlfl . The well known strategy to deal with this "penguin problem" 
is the isospin analysis of Gronau and London [302]. It requires however the measurement of 



Br(B° —* 7r°7r°) which is expected to be below 10 6 : a very difficult experimental task. For 
this reason several, rather involved, strategies |303| have been proposed which avoid the use of 



Bd — ► 7r°7r° in conjunction with acp(7r + 7r~, t). They are reviewed in [Q, ||, It is to be seen 
which of these methods will eventually allow us to measure a with a respectable precision. 

The most promising at present appears the method of Quinn and Snyder |303] ]. It uses the 
Dalitz plot for B® — > qtt — > tt + tt^tt°. This method allows to determine both sin 2a and cos 2a, 
reducing possible discrete ambiguities. The remaining discrete ambiguity can be removed with 
the help of CP asymmetry in B® — ► 7r + 7r _ and a theoretical assumption 304 1. 



While the determination of a remains a challenge for both theorists and experimentalists, 
the situation looks more promising now than in 1999. First the new measurements of Br(B® — ► 
7r + 7r~) from BaBar and Belle 

f (9.3 ±2.4 ±1.3)- 10 -6 (BaBar) [3051 , , 

Br(B° d -> 7T+7T- = I V ' K J 8.7 

^ (6.3 ±4.0) • 10- 6 (Belle) goj , 



show that this decay mode may be less suppressed relative to B® — > K + ir~ than was suggested 
by CLEO, implying smaller penguin pollution. In addition the QCD factorization approach 
Wu, |307| could help in calculating the penguin pollution from first principles. 



8.3.2 The Angle (3 

The CP-asymmetry in the decay B^ — ► i/jKs allows in the SM a direct measurement of the angle 



(3 in the unitarity triangle without any theoretical uncertainties [305]. The relevant formula is 



given in fl3.4| ). Of considerable interest |298| , |3Qg| 1 is also the pure penguin decay Bd —>■ 4>Ks, 
which is expected to be sensitive to physics beyond the SM. Comparision of (3 extracted from 
Bd — ► <pK$ with the one from B^ — > ipK$ should be important in this respect. An analogue 
of Bd — > ipKs in i? s -decays is B s — > ilxp. As shown in (p.6|), the corresponding CP asymmetry 



measures here r\ [299] in the Wolfenstein parametrization. It is very small, however, and this 
fact makes it a good place to look for the physics beyond the SM. In particular the CP violation 
in B® — B® mixing from new sources beyond the Standard Model should be probed in this decay. 
Another useful channel for (5 is B& — > D + D~ . 
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8.3.3 The Angle 7 



The two theoretically cleanest methods for the determination of 7 are: i) the full time dependent 
analysis of B s — ► DfK~ and B s — ► D~K + |310| 1 and ii) the well known triangle construction due 
to Gronau and Wyler Jttl] which uses six decay rates B ± -> Z^pif*, 5+ -> 1)°^+ 
and -B — > D°K~, D°K~. Both methods are unaffected by penguin contributions. The first 
method is experimentally very challenging because of the expected large B® — B® mixing. The 
second method is problematic because of the small branching ratios of the colour supressed 
channel B + — > D°K + and its charge conjugate, giving a rather squashed triangle and thereby 
making the extraction of 7 very difficult. Variants of the latter method which could be more 
promising have been proposed in [312, 313]. Very recently the usefulness of B c — * DD S for the 
extraction of 7 was stressed in [314 1. It appears that these methods will give useful results at 
later stages of CP-B investigations. In particular the first and the last method will be feasible 
only at LHC-B. Other recent strategies for 7 will be discussed below. 



8.4 Constraints for 7 from B — > irK 

The most recent developments are related to the extraction of the angle 7 from the decays B — * 
PP (P=pseudoscalar). Several of these modes have been observed by the CLEO collaboration 
POQ J. In the summer of 2000 BaBar [305] and Belle [315] announced their first results for 
B — > Ktt branching ratios. In the future they should allow us to obtain direct information on 7. 
At present, there are only experimental results available for the combined branching ratios of 
these modes, i.e. averaged over decay and its charge conjugate, suffering from large uncertainties. 



They are collected in table 11. 



Table 11: Branching ratios for B — > Ktt Values in units of 10 



Decay 


CLEO 


BaBar 


Belle 


B° -> tt^± 
B ± 

B° d - n°K 


17.2±|;| ± 1.2 

11 fi +3.0+1.4 
11.0_2.7-1.3 

18.2i£o ± 1.6 

14 ft+5.9+2.4 


19 r+3.0+1.3 
iz ' d -2.6-1.7 


17.4+1-1 ± 3.4 
18.8l|| ± 2.3 

91 n +9.3+2.5 
zi - u -7.8-2.3 



There has been a large activity in this field during the last three years. The main issues here 
are the final state interactions (FSI) R316| , SU(3) symmetry breaking effects and the importance 
of electroweak penguin contributions. Several interesting ideas have been put forward to extract 



the angle 7 in spite of large hadronic uncertainties in B — > irK decays [317, |318], 31S, 32C, 321 
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322]. Reviews can be found in g2j| |1| gg§| . 

Three strategies for bounding and determining 7 have been proposed. The "mixed" strat- 
vr ^ and B ± -> 71^ if. The "charged" strategy |S2l involves 5 ± -> 



egy [fl7| uses B l d 



it® K^~, ir^K and the "neutral" strategy J320| ] the modes B® — > tt^K^, tt°K . General parametriza- 
tions for the study of the FSI, SU(3) symmetry breaking effects and of the electroweak penguin 
contributions in these channels have been presented in [ |319| , 320, 321]. 

As demonstrated in [317, 319, 322, 320 , |32l| , already CP-averaged B — ► irK branching ratios 
may imply interesting bounds on 7 that may remove a large portion of the allowed range from 
the analysis of the unitarity triangle. In particular combining the neutral and charged strategies 
|320| ] one finds that the most recent CLEO data favour 7 in the second quadrant, which is 
in conflict with the standard analysis of the unitarity triangle as we have seen in section 4. 
Other arguments for cos 7 < using B — * PP, PV and VV decays were given in [300, |325|j. 



Simultaneously to 7, the CLEO data provide some information on strong phases. The present 
pattern of these phases indicates either new-physics contributions to the electroweak penguin 



sector, or a manifestation of large non-factorizable S'C/(3)-breaking effects | 320 |. There is no 
doubt that these strategies will be useful in the future. 

Finally ratios of CP-averaged B — ► ttK and B — ► tttt rates as functions of 7 have been 
studied within the QCD factorization approach in |326[| . Interestingly, the ratio R w of the 
B d — * 7r + 7r~ and B® — * tt^K^ decay rates is in disagreement with the CLEO experimental 
value 0.25 ± 0.10, unless the weak phase 7 were significantly larger than 90°. Similar results 
using QCD factorization approach have been found in |327]]. This is in accordance with the 
findings in [32C, |300| , |325| ], On the other hand the most recent value R n = 0.74 ± 0.29 from 
BaBar |305f| does not necessarily require 7 > 90° and is in accordance with the expectations. 
With improved data, also from Belle, the situation should become clearer already this year. 



8.5 Employing U-Spin Symmetry 



New strategies for 7 using the U-spin symmetry have been proposed in [325, 32£]. The first 
strategy involves the decays B® s — ► ipKs and B^ s — > D^ S DJ S J32g[ |. The second strategy 
involves B® 



K+K- and B° d 



7r + 7r~ [329]. These strategies are unaffected by FSI and 
are only limited by U-spin breaking effects. They are promising for Run II at FNAL and in 
particular for LHC-B. 

A method of determining 7, using B + — > K°it + and the U-spin related processes B d — ► K + tt~ 
and Bg — > ir + K~~ , was presented in [330]. A general discussion of U-spin symmetry in charmless 



B decays and more references to this topic can be found in [331] 



84 



8.6 Summary 

This review together with the general discussion of CP asymmetries in Section 3, was only a short 
excursion into the reach field of CP violation in B decays. Clearly this field being dominated 
by the interplay of flavour and QCD dynamics will give us an insight into sofar poorly explored 
sector of particle physics. It is exciting that in the coming years the new experimental data will 
uncover various patterns of CP asymmetries and branching ratios bringing some order into the 
structure of non-leptonic B decays and hopefully giving some clear signals of new physics. 

9 Looking Beyond the Standard Model 

9.1 General Remarks 

We begin the discussion of the physics beyond the Standard Model with a few general remarks. 
As the new particles in the extensions of the SM are generally substantally heavier than W^, 
the impact of new physics on charged current tree level decays should be marginal. On the other 
hand these new contributions could have an important impact on loop induced decays. From 
these two observations we conclude: 

• New physics should have only marginal impact on the determination of \V US \, \V c b\ and 
\V ub \. 

• There is no impact on the calculations of the low energy non-perturbative parameters B{ 
except that new physics can bring new local operators implying new parameters B{. 

• New physics could have substantial impact on rare and CP violating decays and conse- 
quently on the determination of the unitarity triangle. 

9.2 CP Violation Beyond the SM 

The pattern of CP violation in the extensions of the SM deviates generally from the KM picture 
of CP violation. Let us discuss it briefly by comparing the special features of CP violation in 
the SM with those which appear in its extentions. A more detailed discussion can be found in 

HI- 

• In the SM CP is explicitly broken through complex phases in the Yukawa couplings. In 
the extensions of the SM CP can also be spontaneously broken. The latter case takes place 
if the scalar vacuum expectation values contain phases which cannot be removed. 
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• In the SM there is a single complex phase (S or if). In the extensions of the SM new 
complex phases can be present. Examples are models with extended Higgs sector and 
super symmetric models. 

• In the SM CP violation occurs only in charged current weak interactions of quarks, which 
are necessarly flavour violating. As a consequence CP violation is strongly suppressed 
in neutral current transitions (Z°, 7, H°) as it can appear there only as a loop effect. 
Similarly CP violation is very strongly suppressed in flavour diagonal transitions, which 
implies for instance unmeasurable electric dipole moments. 

• In the extensions of the SM CP violation may occur not only in charged current transi- 
tions but also in neutral current transitions at the tree level. Moreover it can occur in 
flavour diagonal interactions. Such interactions can be found in supersymmetric models 
and generally in models with an extended Higgs sector. Moreover flavour violating pro- 
cesses mediated by gluinos (0(a s )) can be CP-violating. Consequently in contrast to 
the SM large CP-violating effects in neutral and flavour diagonal transitions are possible. 
Good examples are substantial, possibly in the near feature measurable, electric dipole 
moments in some supersymmetric models and models with an extended Higgs sector. 

9.3 Classification of New Physics 

Classification of new physics contributions can be done in various ways. We find it useful to 
classify these contributions from the point of view of the operator structure of the effective weak 
Hamiltonian, the complex phases present in the Wilson coefficients of the relevant operators and 
the distinction whether the flavour changing transitions are governed by the CKM matrix or by 
new sources of flavour violation. 

Let us then group the extensions of the SM in five classes. For the first four classes we 
assume that there are only three generations of quarks and leptons. The last class allows for 
more generations. 

Class A 

• There are no new complex phases and flavour changing transitions are governed by the 
CKM matrix. 

• There are no new operators beyond those present in the SM. 

• The Wilson coefficients of the SM operators receive new contributions through diagrams 
involving new internal particles. 
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These new contributions will modify the SM expressions for rare and CP violating decays but 
these modifications can be formulated in a very transparent manner as we will see below. The 
presence of new contributions will have generally impact on the determination of a, f3, 7, \Vtd\ 
and ImAt- In an analysis of the unitarity triangle that uses the SM formulae of Section 4 these 
new contributions will be signaled by 

• Inconsistencies in the determination of (g, fj) through e, B® d — B® d mixing and rare decays. 

• Disagreement of (g, fj) extracted from loop induced decays with (g, ff) extracted using CP 
asymmetries. 

Examples are the Two Higgs Doublet Model II and the constrained MSSM if tan (3 is not too 
large. This class of models, to be named MFV-models 1332] , |333| ] (MFV= Minimal Flavour 
Violation) will be discussed in some detail below. 
Class B 

This class of models differs from class A through the contributions of new operators not 
present in the SM. It is assumed, however, that no new complex phases beyond the CKM phase 
are present. Examples are again the two Higgs doublet model II and the constrained MSSM 
with large tan/3 and all new phases set to zero. 

Class C 

This class of models differs from class A through the presence of new complex phases in the 
Wilson coefficients of the usual SM operators. Contributions of new operators can be, however, 
neglected. In these models new flavour changing transitions appear that are not governed by 
the CKM matrix. An example is the MSSM with not too a large tan/3 and with non-diagonal 
elements in the squark mass matrices. 

This kind of new physics will also be signaled by inconsistencies in the (g, fj) plane. However, 
new complications arise. Because of new phases CP violating asymmetries measure generally 
different quantities than a, (3 and 7. For instance the CP asymmetry in B — > ipK$ will no longer 
measure but (3 + 9np where 6np is a new phase. Strategies for dealing with such situations 
have been developed. See for instance [244, 334, 335, 336 1 and references therein. 



Class D 

Here we group models with new complex phases, new operators and new flavour changing 
contributions which are not governed by the CKM matrix. The phenomenology in this class of 



models is more involved than in the classes B and C [335, 337]. 

Examples of models in classes C and D are multi-Higgs models with complex phases in 
the Higgs sector, general SUSY models, models with spontaneous CP violation and left-right 
symmetric models. 
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Class E 

Here we group the models in which the unitarity of the three generation CKM matrix does 
not hold. Examples are four generation models and models with tree level FCNC transitions. 
If this type of physics is present, the unitarity triangle does not close or some inconsistencies in 
the (g, fj) plane take place. 

Clearly in order to sort out which type of new physics is responsible for deviations from the 
SM expectations one has to study many loop induced decays and many CP asymmeteries. Some 



ideas in this direction can be found in [244, 334]. 



9.4 Models with Minimal Flavour Violation 
9.4.1 General Remarks 

We will assume in accordance with the experimental findings that all new particles have masses 
higher than M\y. As in this class of models the effective local operators are the same as in the 
SM, the impact of new contributions in a renormalization group analysis is then only felt in the 
initial conditions for the Wilson coefficients taken usually at /i = O(Mw). The renormalization 
group transformation from [i = 0(Mw) down to \x = 0(1 GeV) is on the other hand the same as 
in the SM. This simplifies the inclusion of QCD corrections in the new contributions considerably. 
The only thing one has to do is to calculate QCD corrections to the new contributions in the full 
theory. The remaining QCD corrections present in the effective theory are the same as in the 
SM. Similarly all hadronic matrix elements and the related non-perturbative parameters are as 
in the SM. 

In view of the special manner in which the new contributions affect the decay amplitudes, 
it is useful to use the penguin-box expansion discussed in section 2.6 which suppressing CKM 
parameters reads 

A = P Q +Y,PrF r • (9.1) 

r 

The coefficients Pi depend on non-perturbative parameters and include NLO-QCD and QED 
corrections related to effective theory. They are common to the SM and all MFV models 
discussed here. The functions F r resulting from various box and penguin diagrams contain both 
the SM and new physics contributions involving new exchanges such as charged Higgs particles, 
charginos, squarks etc. They depend, in addition to m t , on the masses of new particles such as 
charged Higgs particles, charginos, squarks and sleptons as well as on a number of new physics 
parameters. They include also QCD corrections calculated in the full theory. 

The explicit expressions for F r in the SM without QCD corrections have been listed in 
section 2.6. The corresponding generalization valid for the Two-Higgs-Doublet Model II and 



the MSSM can be found in [10S]. Clearly the calculations of QCD corrections in the full and 



effective theories must be compatible with each other. This is, however, well understood by now 
II- 

Before presenting some specific results in a restricted MSSM with minimal flavour violation, 
we would like to discuss two specific features of this class of models. 

9.4.2 Universal Unitarity Triangle 

Even if this class of extentions of the SM is very simple, the extraction of the CKM parameters in 
these models appears at first sight to be more complicated than in the SM because new physics 
contributions depend on unknown parameters, like the masses and couplings of new particles, 
that pollute the extraction of the CKM parameters. 

Yet as pointed out recently |333[| , in this class of extensions of the SM it is possible to 
construct measurable quantities that depend on the CKM parameters but are not polluted by 
new physics contributions. This means that these quantities allow a direct determination of the 
"true" values of the CKM parameters which are common to the SM and this particular class of 
its extensions. Correspondingly there exists a universal unitarity triangle common to all these 
models. 

The determination of this universal unitarity triangle and of the corresponding CKM pa- 
rameters has four virtues: 

• The CKM parameters can be determined without the knowledge of new unknown param- 
eters present in these particular extensions of the SM. 

• Because the extracted values of the CKM parameters are also valid in these models, the 
dependence of various quantities on the new parameters becomes more transparent. In 
short: the determination of the CKM parameters and of the new parameters can be 
separated from each other. 

• The comparison of the predictions for a given observable in the SM and in this kind of 
extensions can then be done keeping the CKM parameters fixed. 

• Interestingly, the extraction of the universal CKM parameters is essentially free from 
hadronic uncertainties. 

In what follows we will list the set of quantities which allow a determination of the universal 
unitarity triangle. Subsequently we will indicate how the models in this class can be distinguished 
from each other and from more complicated models which bring in new complex phases and new 
operators. We refer to |333fl for details. 
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The strategy is very simple. As the new physics contributions enter only through the func- 
tions F r , one has to look for certain ratios in which these functions cancel out. 

The determination of Rt, one of the sides of the universal unitarity triangle, can be best 
achieved by using the ratio AA^/AM S . Indeed as discussed in section 4.1 one has to a very 
good accuraccy 



V, 



til 



Is 



A" 



1 + A 2 (2£-1) 



A 2 i? 2 



Consequently using 



\Vtd\ 



\v t 



ts 



imBs 



AM ' 



FbaBb, 



(9.2) 



(9.3) 



FB d \JB Bd 

one can determine Rt independently of new parameters characteristic for a given model and of 
wit. If necessary the 0(A 2 ) corrections in ( |9.2j ) can be incorporated. This will be only required 
when the error on £ will be decreased below 2%, which is clearly a very difficult task. 

While the ratio AM^/ AM S will be the first one to serve our purposes, there are at least two 
other quantities which allow a clean measurement of Rt within the class of the extensions of the 
SM considered. These are the ratios 

2 



Br{B -> X d vv) 



Br{B - 
Br{B d -» M V 



X s vv) 



Vt 



ts 



r Bd m Bd F% d 



Vt 



td 



Vt 



ts 



(9.4) 



(9.5) 



Br{B s -> n+ir) t Bs m Ba F% a 
which similarly to AM<f/AM s measure | V t d\/\ Vt s \ As discussed in section 7, out of these three 
ratios the cleanest is (p.4|), which is essentially free of hadronic uncertainties. Next comes 
(|9.5| ), involving 577(3) breaking effects in the ratio of £?-meson decay constants. Finally, SU(3) 
breaking in the ratio of bag parameters BB d /BB s enters in addition in Q9.3j). These SU(3) 
breaking effects should eventually be calculable with high precision from lattice QCD. 

Note that the branching ratio for the rare decay K + — > n + vi> provides a clean measurement of 
Vtd and consequently of Rt as discussed in section 6. However, this branching ratio alone cannot 
serve our purposes because it is sensitive to new physics contributions through the function X. 

In order to complete the determination of g and fj in the universal unitarity triangle one can 
use sin 2(3 extracted either from the CP asymmetry in Bj — * ipKs or from K — > nvv decays. 
We have discussed the relevant formulae in sections 8 and 6 respectively. They remain valid in 
the models in question as there are no new complex phases involved. As these formulae do not 
depend on the functions F r , both extractions of sin 2/3 are independent of the new parameters 
characteristic for a given model. 
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Concerning the determination of the angle 7, any method for the determination of 7 in 
which new physics of the type considered here can be eliminated or neglected can be used here, 



in particular the methods in [310, [311| , 312, 313] that involve only tree diagrams. 



Once Rt and sin 2(3 have been determined as discussed above, g and rj can be found through 



H |9|] 



R 1 

fj = a—=Jsm2(3 ■ r_b(sin2/3) , g = 1 — fjrb (sin 2/3) (9.6) 
V2 v 



where 



r b (z) = (l + b\/l- z 2 )/z, a,b = ±. (9.7) 
Thus for given values of (Rt, sin 2(5) there are four solutions for (g,fj) corresponding to (a, b) = 



(+, +), (+, — ), (— , +), (— , — ). As described in | 338 | three of these solutions can be eliminated 
by using further information, for instance coming from | V^t,/ V^f, | and e, so that eventually the 
solution corresponding to (a, b) = (+, +) is singled out, 

R 1 

fj = — I Jsin2/3- r_(sin2/3) , q = 1 - fjr+ (sin 2/3) . (9.8) 
y2 v 



A numerical analysis can be found in [332] 



On the other hand g and ij following from Rt and 7 are simply given by 

?7 = i?fesin7, £ = i4cos7 (9-9) 



with 



R b = cos 7 ± dR 2 t - sin 2 7. (9.10) 



Comparing the resulting R b with the one extracted from [F^b/V^l (see ( 1.39Q ) one of the two 
solutions can be eliminated. Similarly using [T/^/V^I and 7 one can construct the universal 
unitarity triangle by means of ( p.9| ). 

Alternatively one could use the measurement of R b by means of [V^/V^l together with Rt 
from AM d /AM s to find 



fj = yjR* - e 2 , q = 1(i + rI-rI). (9-ii) 

This determination suffers from hadronic uncertainties in the extraction of [Ktb/V^l but could 
turn out to be useful once the uncertainties in |14fe/^cfe| have been reduced. We will give 
numerical examples below. 

We observe that all these different methods determine the "true" values of fj and g indepen- 
dently of new physics contributions in the class of models considered. Since A and |V^,| = A}? 
are determined from tree level K and B decays they are insensitive to new physics as well. Thus 
the full CKM matrix can be determined in this manner. The corresponding universal unitarity 
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triangle common to all the models considered can be found directly from formulae like (|K 
(fTp and ( plU) . 



Having determined the CKM parameters by one of these methods one can calculate e, e' /e, 
AAfrf, AM, and branching ratios for rare decays. As these quantities depend on the parameters 
characteristic for a given model the results for the SM, the MSSM and other models of this class 
will generally differ from each other. Consequently by comparing these predictions with the data 
one will be able to find out which of these models is singled out by experiment. Equivalently, e, 
e' /e, AMj, AM S and branching ratios for rare decays allow to determine non-universal unitarity 
triangles that depend on the model considered. Only those unitarity triangles which are the same 
as the universal triangle survive the test. 

It is of course possible that new physics is more complicated than discussed here and that 
new complex phases and new operators beyond those present in the SM have to be taken into 
account. These types of effects would be signaled by: 

• Inconsistencies between different constructions of the universal triangle, 

• Disagreements of the data with the AMj s and the branching ratios for rare K and B 
decays predicted on the basis of the universal unitarity triangle for all models of the class 
considered here. 

In our opinion the universal unitarity triangle provides a transparent strategy to distinguish 
between the MFV models and to search for physics beyond the SM. Presently we do not know 
this triangle as all the available measurements used in section 4 for the construction of the 
unitarity triangle are sensitive to physics beyond the SM. It is exciting, however, that in the 
coming years this triangle will be known once AM S has been measured and sin 2(5 extracted 
from the CP asymmetry in B^ — ► ipKg. At later stages IV^/V^I, K — ► ttvP, B — > X^^vv^ 
Bd,s an d future determinations of 7 through CP asymmetries in B decays will also be 

very useful in this respect. 

9.4.3 sin 2/3 from AM S and \V ub /V cb \ 

In view of the forthcoming precise measurements of (sin 2/3) ^ s it is of interest to investigate 



its dependence on AM S and |^b/V^b| in the MFV models. We show this in table 12. To this 
end we have set £ = 1.16. The interesting feature is a very weak dependence on AM S for high 
values of |T^b/V^,|. For low values of (V^/V^l this dependence is rather strong. We observe 
that sin 2(5 has to be above 0.5 when AM S < 18/ps. On the other hand for AM, = 26/ps it can 
be as low as 0.17. The question then arises whether such low values are consistent with other 
known measurements. 
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Table 12: Values of sin 2(3 in the MFV models for specific values of AM S and \ V u b/V c b 



AM s /\V ub /V cb \ 


0.065 


0.075 


0.085 


0.095 


0.105 


17.0/ps 


0.528 


0.608 


0.682 


0.749 


0.809 


18.0/ps 


0.514 


0.599 


0.677 


0.747 


0.809 


19.0/ps 


0.496 


0.587 


0.669 


0.742 


0.807 


20.0/ps 


0.473 


0.572 


0.659 


0.736 


0.803 


22.0/ps 


0.412 


0.531 


0.632 


0.718 


0.793 


24.0/ps 


0.322 


0.476 


0.595 


0.694 


0.777 


26.0/ps 


0.165 


0.403 


0.549 


0.664 


0.757 



9.4.4 An Absolut Lower Bound on sin 2/3 



It turns out that in the MFV models there exists an absolute lower bound on sin 2(5 that follows 
from the interplay of AM, and e and depends only on \V c b\, IK^/V^I and the nonperturbative 



parameters Bx, FB d y and £ [147]- The derivation of this bound is straightforward. We first 
generalize the SM expressions for the e-hyperbola in ( |4.81 ) and Rt resulting from AMj in ( |4.83| ) 



and (4.84) to the MFV models simply as follows 



fj (1 - Q)A 2 rj 2 F tt + P c (e) A 2 B K = 0.226 



Rt = 1.26 



Ro 1 

A VWt 



Ro = 1.03 



/ (AM), 
0.50/ps 



200 MeV 
\[B~ d F Bd 



'0.55 



(9.12) 



(9.13) 



That is So(xt) in the SM formulae and given in ( 2.2C| ) is just replaced by Ftt, that now in 
addition to the contributions of the box diagrams with top quark exchanges includes all possible 
new physics contributions within a given MFV model. 

The most important feature of the formulae ( p. 12 ) and (9.13) relevant for the discussion 
below is that in the context of the standard analysis of the unitarity triangle the different MFV 
models can be characterized by the value of the function Ftt- This has been stressed recently 



in particular by Ali and London [137 1 . Moreover, as shown in [147], the new physics effects 
cancel in the ratio 772 jf]B ■ In view of this it is convenient in what follows to use the SM values 
772 = 0.57, t\b = 0.55 and absorb all QCD corrections related to new physics contributions into 
Ftt- 

Noting that 

sin2/3 = MLzl) (9 . 14 ) 



Rt 
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and combining ( |9.12| ) and ( p. 13 ) we find [^7], 147 1 

1.26 



sin 2(3 



Rim 



0-226 



K 



(9.15) 



wherby the first term in the parenthesis is typically by a factor of 2-3 larger than the second 
term. This dominant term is independent of F tt and involves the QCD corrections only in the 
ratio 772/ f]B- Consequently it is independent of mt and the new parameters in the extensions 
of the SM. The dependence on new physics is only present in fj entering the second term that 
would be absent if charm contribution to e was negligible. In particular for g > 0, the value of 
fj decreases with increasing Ftt- 

In spite of the sensitivity of the second term in ( |9.15 ) to new physics contributions, there 



exists an absolut lower bound on sin 2(3 in the MFV models, simply because for Bk > the 
unitarity of the CKM matrix implies 

< fj < R b - (9.16) 

At first sight one would think that the lower bound for sin 2/3 is attained for fj = R^, but this is 
clearly not the case as fj depends on the values of the parameters in ( 4.98; ). Consequently there 
is a correlation between the values of the two terms in ( 9.15| ) . 



In fig. 11 we show (sin2/3) m ; n as a function of Ftt obtained in |147| 1 by means of the scanning 
method for slightly different set of input parameters than given in table |2[ To this end the 
standard analysis of section 4 has been used except that So(xt) has been replaced by F U - While 
with increasing F tt < 6.2, (sin 2(3) m \ n decreases, it increases for larger values of F tt , so that indeed 
an absolut minimum for sin 2(3 in MFV models is found. In the case of the future ranges of the 



input parameters considered in |147| , the minimum is found for Ftt ~ 5.2. For Ftt > 13.5 (7.8) 
in the case of the present (future) input parameters no solutions for sin 2(3 are found. 



A number of super symmetric MFV models has been reviewed by Ali and London [137] 



where references to the original literature can be found. Using the results of [137] we find 
as characteristic values Ftt = 3.0, Ftt = 3.4, Ftt = 4.3 for minimal SUGRA models, non- 
minimal SUGRA models and non-SUGRA models with EDM constraints respectively. Moreover, 
Ftt = 2.46 and Ftt = 5.2 for the SM and the MSSM version of [|109([ respectively. We observe a 
rather weak dependence of (sin2/3) m i n on Ftt- This is in agreement with the analysis of [137], 



where sin 2(3 has been studied in the range 2.5 < Ftt < 4.3. Consequently the measurement of 
(sin 20)^k s will not be able to distinguish easily different MFV models. 

On the other hand, the existence of an absolut bound on sin 2(3 in the MFV models allows 
in principle to rule out this class of models if (sin 2(3)^k s is found below (sin 2/?) m i n . The impli- 
cations of such a possibility would be rather profound. With the measurement of (sin 2/3)^ s 
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Figure 11: Lower bound for sin 2/3 as a function of F tt for present and future ranges of the input 



parameters [147|. 



alone one would be able to conclude that new CP violating phases and/or new local operators in 
the weak effective Hamiltonians for K° — K° and s — B^ s mixings are necessary to describe 
the data. 



Repeating the analysis of [147] we find using the input parameters of table || 

(sin 2/3) min = 0.42 (9.17) 



to be compared with 0.34 in [147]. We would like to emphasize that this bound should be con- 
sidered as conservative. Afterall it has been obtained by scanning independently all parameters 
in question. The main reason for the improved bound is the higher value of FB d \J~~B~d used in 
the new analysis. The plot of (sin2/3) m i n versus Fu in the case of the parameters of table § is 



given in fig. 12 where this time the "future" represents simply the parameters of table H with 



the constraint F Bd \l B d > 210 MeV. 



The anatomy of the bound is given in table 13. Here we show (sm2f3) m - m as a function of 
Bk and |V^>| with all the remaining parameters scanned within the ranges of table ^. In the 
fourth column we show the impact of the reduced uncertainty in FB d \f~B~d. In the fifth column 
the impact of the future measurement of AM S is shown. This measurement giving an upper 
bound on AM S will provide the lower bound on Rt by means of ( |4.87| ) in addition to the known 
upper bound. This in turn will exclude high values of Ftt, see ( |9.13|) , implying an improved 
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Figure 12: Updated lower bound for sin 2(5 as a function of Fu for present and future ranges of 
the input parameters. 



lower bound on sin 2(3. The chosen value £ e g > 1.05 with £ e g defined in (4.87) corresponds for 
instance to (AM) S < 18.0/ps for £ = 1.15. The numbers in the parentheses are explained below. 

On the basis of this table and additional numerical analysis we find the following features in 
accordance with ( |9.15| ): 

• (sin2/?) m j n increases with decreasing Bk and \V c b\ and increasing |Vufe/V c {,| and Fb cL \Jb~ ( i. 

• In the ranges considered, the dependence of (sin2/3) m i n on Bk, \V c b\ and Fs d \fB~d is 
stronger than on |Kt6/^cb|- This is evident from fl9.15| ) in which |Vu&Ac&| is n °t explicitly 
present but affects the bound only through the value of f) in the subleading term and 
indirectly through its impact on the allowed ranges of the remaining parameters. 

• On the other hand when the upper bound on AM S will be known, the lower bound on 
|Kife/^cfe| wih have an important on (sin2/3) m i n . This is seen in table ^ and in the fifth 
column of table 13 where the values in the parentheses show the impact of the lower bound 
l^ufe/V^I > 0.085. In this case the dependence of the bound on Bk and \V c i,\ essentially 
disappears in the ranges considered. 

• One can check that (sin2/3) m i n is roughly proportional to Fg B^ and its dependence on 
Bk and A can be approximately given by a single variable r = A 2 Bk- These features 
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are clear from ( p. 15 ), table [l3| and fig. 12. The observed small departures from these 
regularities are caused by the correlations between various parameters as discussed below 



( )9lq ). 

With the parameters of table the impact of the measurement of AM S is only felt for 



Bk > 0.85. It was larger in the case of parameters considered in [147]. With increasing 
£ e ff the impact of the measurement of AM S becomes larger. 



Table 13: Values of ( sin 2/?) m i n in the ]VIFV models for specific values of Bj^ and |V^| and. 
different scenarios as explained in the text. 



Bk 


\Vcb\ 


Present 


F Bd ^B d > 210 MeV 


&// > 1-05 




0.039 


0.71 


0.82 


0.71 (0.71) 


0.70 


0.041 


0.65 


0.75 


0.65 (0.70) 




0.043 


0.59 


0.69 


0.59 (0.70) 




0.039 


0.59 


0.68 


0.59 (0.70) 


0.85 


0.041 


0.54 


0.63 


0.54 (0.70) 




0.043 


0.50 


0.57 


0.50 (0.70) 




0.039 


0.51 


0.59 


0.56 (0.70) 


1.00 


0.041 


0.46 


0.54 


0.52 (0.70) 




0.043 


0.42 


0.49 


0.49 (0.70) 



At present the lower bound on sin 2(3 is fully consistent with the experimental findings. On 
the other hand this bound could become stronger when our knowledge of the input parameters 
in question improves and when the upper bound on AM S will be experimentally known. In 
particular if the upper bounds on Bk and \ V c b\ and lower bounds on IK^/V^I, FB d \jB~a and £ c g 
will be improved, (sin2/3) m i n could be shifted above 0.5. If simultaneously the future accurate 
measurements of a^K s wm confirm the low values reported by BaBar and Belle, all MFV models 
will be excluded. 

Clearly other measurements, in particular those of the rare decay branching ratios and various 
CP asymmetries, will have an additional impact on (sin 2/3) m j n , but this is a different story. 

Finally as pointed out in [147], the absolute lower bound on sin 2/3 implies within the MFV 
models an absolute lower bound on the angle 7. We find 

(sin 7 ) min = 0.31 (9.18) 

with 7 in the first quadrant. The second quadrant in the MFV models is excluded through the 
lower bound on AM,. 
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It will be exciting to watch the experimental progress in the values of a^x s and AM, and 
the theoretical progress on Bk, \V u b/V c b\, FB d \[~B~d and £. Possibly we will know already next 
summer that new CP violating phases and/or new operators in the effective weak Hamiltonians 
are mandatory. 

On the other hand if (sin2/3),^ s will be found above 0.5 the MFV models will remain 
as a vital alternative to more complicated models unless one can exclude them by different 
measurements as discussed in the section on the universal unitarity triangle. In this context 
and in the distinction between different MFV models the measured value of AM S will play an 
important role as \Vts\ ~ l^cfel and AM S is directly proportional to Ftt- To this end Fb s \Jb~ s has 
to be precisely known. 



9.4.5 Results in the MSSM 

There are many new contributions in MSSM such as charged Higgs, chargino, neutralino and 
gluino contributions. However, in the case of minimal flavour and CP violation it is a good 
approximation to keep only charged Higgs and chargino contributions. 

The most recent analysis of e'/e and of rare decays in this scenario can be found in |109|. In 



this analysis constraints on the supersymmetric parameters from e, AM^ S , B — > Xg'y, Ag in the 
electroweak precision studies and from the lower bound on the neutral Higgs mass have been 
imposed. Supersymmetric contributions affect both the loop functions F r present in the SM 
and the values of the extracted CKM parameters like |V^| and ImA(. As the supersymmetric 
contribution to the function Fu in ( 9.12| ) are always positive, the extracted values of \Vtd\ and 



ImAj are always smaller than in the SM. Consequently quantities sensitive to \ Vtd\ and ImAt are 
generally suppressed relative to the SM expectations. Only for special values of supersymmetric 
parameters, the supersymmetric contributions to the loop functions Xq, Yq and Zq can over- 
compensate the suppression of \Vtd\ and ImA( so that some enhancements of branching ratios 
are possible. 

Let us define by T{Q) the MSSM prediction for a given quantity Q normalized to the SM 
result. Setting |V^|, \V c b\ and the non-perturbative parameters Bi, all unaffected by SUSY 
contributions, at their central values one finds f|10£ 



0.53 < T(e'/e) < 1.07, 0.65 < T(K + tt + uu) < 1.02 (9.19) 

0.41 < T(K L n°vp) < 1.03, 0.73 < T(K L 7r°e + e _ ) < 1.10 . (9.20) 

We observe that suppressions by a factor of 2 relative to the SM expectations are still possible. 

As the CKM element \Vt s \ is practically fixed by the unitarity of the CKM matrix and 
consequently unaffected by supersymmetric contributions, the supersymmetric contributions 
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to B — ► XgH/j, and -B s — > /i// enter only through the loop functions. Consequently visible 
enhancements, but also suppressions, in the corresponding branching ratios are found 

0.73 < T(B -> X,np) < 1.34, 0.68 < T(B S -> ^u) < 1.53 . (9.21) 



Reference [109] provides a compendium of phenomenologically relevant formulae in the 
MSSM, that should turn out to be useful once the non-perturbative parameters Bi will be 
better know and the relevant branching ratios have been measured. The study of the unitarity 



triangle can be found in [137]. 



9.5 Going Beyond the MSSM 

In general supersymmetric models the effects of supersymmetric contributions to e'/e and to rare 
branching ratios can be much larger than discussed above. In these models new CP-violating 
phases and new operators are present. Moreover the structure of flavour violating interactions 
is much richer than in the MFV models. The new flavour violating interactions are present 
because generally the sfermion mass matrices can be non-diagonal in the basis in which all 
quark-squark-gaugino vertices and quark and lepton mass matrices are flavour diagonal. Instead 
of diagonalizing sfermion mass matrices it is convenient to consider their off-diagonal terms as 
new flavour violating interactions. Denoting by A the off-diagonal terms in the sfermion mass 
matrices, the sfermion propagators can be expanded as a series in terms of 8 = A/rh 2 , where m 
is the average sfermion mass. As long as A is significantly smaller than m 2 , we can just take the 
first few terms of this expansion and compute any given process in terms of these 8's. This is 
equivalent to an approximate diagonalization of the squark mass matrices around their diagonal 
part. 



The method just described is the so-called mass-insertion approximation [336]. It has been 



reviewed in the classic papers [337, |335|], where further references can be found. The basic 



parameters in this approach are the insertions 

(<%)lL, {5ij)LR, (Sij)RL, (8ij)RR (9.22) 

with i,j = 1,2,3 denoting flavour indices and L,R the helicities of the fermionic partners of 
sfermions. These parameters being generaly complex constitute new sources of CP violation. 
Their values can be constrained through the existing data on flavour violating (i ^ j) and 
flavour conserving (i = j) processes. For instance (812) xy can be constrained through AS = 1 
and AS = 2 transitions, whereas (8is)xy and {82z)xy through analogous transitions involving 
B® and B® mesons, respectively. 

There is a vast literature on the phenomenological applications of the mass insertion method 
to rare and CP-violating decays. Here I would like to review very briefly my own work on this 
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subject in a wonderful collaboration with my Italian friends Gilberto Colangelo, Gino Isidori, 



Andrea Romanino and in particular Luca Silvestrini [ 220 , 228 , 336 1, I n [336] a model inde- 
pendent analysis of K — > irvv with an application to general supersymmetric models has been 
presented. A similar analysis can also be found in pOR . In p§ rare kaon decays in models 



with an enhanced sdZ vertex, for instance supersymmetric models, have been considered. The 
generalization of this analysis to B decays has been presented in |267|] . Let me concentrate here 



on |220|1 . 

Despite the presence of a large number of parameters within the mass insertion approach 
to general supersymmetric models, only a few of them are allowed to contribute substantially 



to e'/e. Phenomenological constraints, coming mainly from AS" = 2 transitions [337], make the 
contribution of most of them to AS" = 1 amplitudes very small compared to the SM one. The only 
parameters which survive are the left-right mass insertions contributing to the Wilson coefficients 
of Z- and magnetic-penguin operators. The reason for this simplification is a dimensional one: 
these are the only two classes of operators of dimension less than six contributing to e'/e. 
Supposing that the enhancement of the Wilson coefficients of either of these two (or both) type 
of operators is responsible for the observed value of e'/e, a corresponding effect in the rare kaon 



decays should be observed. In [22C] the relation between these new contributions to e'/e and 
the corresponding contributions to the rare decays Kl — > tt^uP, K + —* 7r + z^P, Kl — » 7r°e + e _ 
and Kl — > [i + fi~ has been analysed in detail. 

Before presenting the results of this analysis let me recall the observation of Colangelo and 



Isidori [230] that the branching ratios of rare kaon decays could be considerably enhanced, in a 
generic supersymmetric model, by large contributions to the effective sdZ vertex due to a double 
left-right mass insertion. This double mass insertions had not been included in the analyses 



[336, 340|, where only single mass insertions were taken into account. Consequently only modest 
enhancements of rare decay branching ratios, up to factors 2-3 at most, can be found in these 
papers, as opposed to the possible enhancement of more than one order of magnitude found in 



[23C]. This interesting observation has been challenged by Silvestrini and myself in [ |228| ], where 



we have shown that the data on e'/e and Kl — > fip, may constrain considerably the double 
left-right mass insertion and the corresponding enhancement of the rare kaon branching ratios. 
Our model independent analysis which went beyond supersymmetry, but assumed that the only 
new effect is an enhanced sdZ vertex, resulted in the following bounds: 

Br(K + -> tt + uv) < 2.3 • 10~ 10 , Br(K L -> tt°uP) < 2.4 • 10" 10 (9.23) 

and Br(K^ —* 7r°e + e~) < 3.6 • 10 -11 , which are substantially stronger than the bounds found 
in ]p3C 
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The purpose of our joint collaboration [|22C| 1 was to use the strategy in p2q] specifically in 
the case of super symmetry and to include also the effects of chromomagnetic and 7-magnetic 
penguins to e'/e and — > 7r°e + e _ , respectively. The latter two contributions were assumed 
to be small in the models considered in |228|]. Moreover, we have investigated whether the 



large double mass insertions suggested in [230] could be further constrained within the specific 
framework of supersymmetry. 

As demonstrated in [|220|| , radiative effects relate the double left-right mass insertion to the 
single left-left one. The phenomenological constraints on the latter imply then a stringent bound 
on the supersymmetric contribution to the Z-penguin. Using this bound and those coming from 
the data on e'/e and Kl — > /ifi one finds 

Br(K + -> 7r + i/i/) < 1.7- 10~ 10 , Br(K L -» vr ^) < 1.2 • 10~ 10 (9.24) 



and Br(K\j — > 7r°e + e~) < 2.0-10~ n , that is slightly stronger bounds than given in (9.2S). Larger 
values are possible, in principle, but rather unlikely. Thus the most probable maximal enhance- 
ments over the SM expectations are roughly by factors 2, 4 and 4 for these three branching 



ratios, respectively. This analysis confirms the previous findings in [21£] that the most natural 
enhancement of e'/e, within supersymmetric models, comes from chromomagnetic penguins. In 
this case sizable enhancement of Br(Ki, — ► 7r°e + e _ ), as seen above, can be expected. 

There is of course a number of other analyses of general supersymmetric effects in rare and 
CP-violating decays. Some references are collected in [341]. In particular, the constraints on 



phases of supersymmetric flavour conserving couplings can be best obtained from the upper 



bounds on electric dipole moments [342| 



9.6 Spontaneous CP Violation 

Spontaneous CP violation (SCPV) is a very interesting topic that we cannot cover in these 
lectures. It requires extended Higgs sector. A very comprehensive presentation of SCPV can be 
found in M where left-right symmetric models and multi-Higgs models are discussed in detail. 



Some recent papers on SCPV are given in [342] 



10 Perspectives 

I hope I have convinced the students that the field of CP violation and rare decays plays an 
important role in the deeper understanding of the SM and particle physics in general. Indeed the 
field of weak decays and of CP violation is one of the least understood sectors of the SM. Even if 
the SM is still consistent with the existing data for weak decay processes, the near future could 
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change this picture considerably through the advances in experiment and theory. In particular 
the experimental work done in the next ten years at BNL, CERN, CORNELL, DA^NE, DESY, 
FNAL, KEK, SLAC and eventually LHC will certainly have considerable impact on this field. 

Let us then make a list of things we could expect in the next ten years. This list is certainly 
very biased by my own interests but could be useful anyway. Here we go: 

• The error on the CKM elements \V c b\ and \V u bl^d\ w hl be decreased below 0.0015 and 
0.01, respectively. This progress should come mainly from CLEO III, ^-factories and 
new theoretical efforts. It will have considerable impact on the unitarity triangle and 
will improve theoretical predictions for rare and CP-violating decays sensitive to these 
elements. 

• The error on mt should be decreased down to ±3 GeV at Tevatron in the Main Injector 
era and to ±1 GeV at LHC. 

• The measurement of a non- vanishing ratio e'/e by NA31, KTeV and NA48, excluding 
confidently the superweak models, has been an important achievement. Yet, as we have 
discussed in section 5, the experimental value for e'/e requires considerable improvements 
before we could be satisfied with it. In particular, the difference between KTeV and NA48 
results by roughly a factor of two is rather disturbing. Let us hope that some progress will 
be made in this direction in the coming years. It should be stressed that the determination 
of e'/e with the accuraccy of ±(1 — 2) • 10 -4 from NA48, KTeV and KLOE will give some 
insight into the physics of direct CP violation inspite of large theoretical uncertainties. 
In this respect measurements of CP-violating asymmetries in charged B and K decays 
will also play an outstanding role. The situation concerning hadronic uncertainties is 
quite similar to e'/e. Therefore one should hope that some definite progress in calculating 
relevant hadronic matrix elements will also be made. 

• One of the most exciting measurements in this year will be the measurement of AM S . 
LEP and SLD experiments have done already a fantastic progress by providing the lower 
bound AM S > 15/ps. The actual measurement of AM S should come first from Run II at 



FNAL. With the improved calculations of £ in ( 4.85 ) this will have an important impact 
on the determination of \Vtd\, on the unitarity triangle and as discussed in section 9 on 
models with minimal flavour violation. 

Clearly future precise studies of CP violation by BaBar, Belle, CDF, DO, CLEO III, LHC-B 
and BTeV providing direct measurements of a, [3 and 7 may totally revolutionize our field. 
The first results from CDF, BaBar and Belle are very encouraging. During the recent years 
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several, in some cases quite sophisticated and involved, strategies have been developed to 
extract these angles with small or even no hadronic uncertainties. Certainly the future 
will bring additional methods to determine a, (3 and 7. Obviously it is very desirable to 
have as many such strategies as possible available in order to overconstrain the unitarity 
triangle and to resolve certain discrete ambiguities which are a characteristic feature of 
these methods. A recent review on discrete ambiguities with the relevant literature can be 



found in [|344]. 



Improved data for K + — > n^vv should be reported by AGS E787 collaboration this year. 
In view of the theoretical cleanliness of this decay the measured branching ratio at the 
1.5 • 10~ 10 level would signal physics beyond the SM. A precise measurement of this very 
important decay requires, however, new experimental ideas and new efforts. The new 



efforts [241, p42[ | in this direction allow to hope that a measurement of Br(K + — > tt + uu) 
with an accuracy of ±10% should be possible before 2005. This will provide an important 
test for the SM and its extensions. 



The newly approved experiment at BNL |250| 1 and the planned experiments at KEK | 24£ 



and FNAL [249] to measure Br(K^ — > ifivv) at the ±10% level before 2005 may make a 



decisive impact on the field of CP violation. — > -k^vv allows the cleanest determination 
of ImAf and taken together with K + — > tt + uu a very clean determination of sin 2(3. 

The theoretical status of — > ir°e + e~ and of — > /j,fi, should be improved to confront 
future data. Experiments at DA<I>NE should be very helpful in this respect. 

The future improved inclusive measurements B — > X S( f) confronted with improved SM 
predictions could give the first signals of new physics. It appears that the errors on the 
input parameters could be lowered further and the theoretical error on Br{B — > X s j) 
could be decreased confidently down to ±8% in the next years. The same accuracy in 
the experimental branching ratio will hopefully come from BaBar and Belle. This may, 
however, be insufficient to disentangle new physics contributions although such an accuracy 
should put important constraints on the physics beyond the Standard Model. It would 
also be desirable to look for B — > Xj'y, but this is clearly a much harder task. 

Similar comments apply to transitions B — > X s l + l~ and B — » K*l + l~ which are much 
reacher and more sensitive to new physics contributions than B — » X s ^- Observations of 
B — > X s fj,p, and B — > K* \x]x are expected from Run II at FNAL and B factories already 
this year. The distributions of various kind when measured should be very useful in the 
tests of the SM and its extensions. 
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The measurement of B — > X s ^vv and -B^d — > //// will take most probably longer time but 
as stressed in these lectures all efforts should be made to measure these transitions. 

According to the SM, CP violation outside the K-meson and B-meson systems is expected 
to be essentially unobservable. On the other hand, new sources of CP violation present in 
multi-Higgs models or models with supersymmetry could give rise to measurable effects. 
Attempts to trace these effects include charm meson decays, D° — D° mixing, hyperon 
decays at Fermilab (E756) R345| and the searches for the electron and neutron electric 



dipole moments. The most interesting at present are the results from CLEO [346] and 



FOCUS (FNAL) [347] which possibly indicate D° — D° mixing at a much higher level 



than expected in the SM and in the ball park of some super symmetric expectations |34£ 



Top quarks and Higgs particles [34S, 350] may also be used to probe CP violation once 



they are produced in large numbers at the Tevatron, the LHC and future linear colliders. 



Finally CP violation in neutrino oscillations 351 ] and in the context of baryogenesis [ 352 ] 



as well as lepton flavour violation [352] belong to the class of phenomena outside the SM. 



Considerable progress in this area should be made in this decade. 

• On the theoretical side, one should hope that the non-perturbative methods will be consid- 
erably improved so that various B{ parameters will be calculated with sufficient precision. 
It is very important that simultaneously with advances in lattice QCD, further efforts are 
being made in finding efficient analytical tools for calculating QCD effects in the long dis- 
tance regime. This is, in particular very important in the field of non-leptonic decays, where 
the progress in lattice calculations is slow. The accumulation of data for non-leptonic B 
and D decays in the coming years should teach us more about the role of non-factorizable 
contributions and in particular about the final state interactions. In this context, in the 
case of K-decays, important lessons will come from DA^NE which is an excellent machine 
for testing chiral perturbation theory and other non-perturbative methods. 

This list of topics shows that flavour dynamics and the related CP violation and rare processes 
have a great future. They will surely play an important role in particle physics in this decade. 
Clearly the next ten years should be very exciting in this field. 
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